Influence of Generic Scale Invariance on Classical and Quantum Phase 
Transitions 



D.Belit£| 

Department of Physics and Materials Science Institute, University of Oregon, Eugene, OR 97403 
T.R.Kirkpatriclfl 

Institute for Physical Science and Technology, and Department of Physics, University of Maryland, College Park, 
MD 20742 

Thomas 

Dept. of Physics, University of Missouri-Rolla, Rolla, MO 65409 
(Dated: February 2, 2008) 

This review discusses a paradigm that has become of increasing importance in the theory of 
quantum phase transitions; namely, the coupling of the order-parameter fluctuations to other 
soft modes, and the resulting impossibility of constructing a simple Landau-Ginzburg- Wilson 
theory in terms of the order parameter only. The soft modes in question are manifestations 
of generic scale invariance, i.e., the appearance of long-range order in whole regions in the phase 
diagram. The concept of generic scale invariance, and its influence on critical behavior, is explained 
using various examples, both classical and quantum mechanical. The peculiarities of quantum 
phase transitions are discussed, with emphasis on the fact that they are more susceptible to the 
effects of generic scale invariance than their classical counterparts. Explicit examples include: 
the quantum ferromagnetic transition in metals, with or without quenched disorder; the metal- 
superconductor transition at zero temperature; and the quantum antiferromagnetic transition. 
Analogies with classical phase transitions in liquid crystals and classical fluids are pointed out, 
and a unifying conceptual framework is developed for all transitions that are influenced by generic 
scale invariance. 
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I. PHASE TRANSITIONS 

Phase transitions are among the most fascinating phe- 
nomena in nature. They also have far-reaching impli- 
cations. The liquid-gas and liquid-solid transitions in 
water, for instance, are common occurrences of obvious 
importance. The transition from a paramagnetic phase 
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6. Other quantum phase transitions where GSI might plav a tttefci lferromagnetic one in the elements iron, nickel, and 

cobalt made possible the invention of the compass. The 
structural transition in tin from the /3-phase (white tin) 
to the a-phase (gray tin) is responsible for the degra- 
dation of tin artifacts, below a temperature of about 
286 K, known as the tin pest. One could continue with a 
long list, but these three examples may suffice to demon- 
strate that phase transitions come in a wide variety of 
phenomenologies with no entirely obvious unifying fea- 
tures. Accordingly, early attempts at theoretically un- 
derstan ding pha se transitions focused on particular ex- 
amples. Ivan der Waala (1873) developed the first exam- 
ple of what was later to become known as a "mean-field 
theory" , in this case, to describe the liquid-gas transi- 
tion. Weiss (1903) gave a mean-field theory of ferromag- 
netism which was based on the concept of the "mean 
field" seen by each 'elementary magnet' (spin had yet to 
be discovered at that time) that is produced by all other 
elementary magnets. A unification of al l mean-field the- 
ories was achieved by iLandaul l)l937albl l3ldl) (all of these 
pape rs, or their translations, are reprinted in iLandauL 
119651) . He introduced the general concept of the "order 
parameter" , a thermodynamic observable that vanishes 
in one of the phases separated by the transition (the 
"disordered phase" ^), and is nonzero in the other (the 
"ordered phase"). In the case of the ferromagnet, the 
order parameter is the magnetization; in the case of the 
liquid-gas transition, the density difference between the 
two phases. Landau theory underlies all later theories of 
phase transitions, and we therefore discuss it first. 



Preamble 

The theoretical understanding of classical, or thermal, 
phase transitions, which occur at a nonzero temperature, 
is very well developed. A characteristic feature of con- 
tinuous phase transitions, or critical points, is that the 
free energy, as well as time correlation functions, obey 
generalized homogeneity laws. This leads to scale invari- 
ance, i.e., power-law behavior of various thermodynamic 
derivatives and correlation functions as functions of the 
temperature, external fields, wave number, time, etc. 
Quantum phase transitions, which occur at zero temper- 
ature as a function of some nonthermal control parameter 
such as composition, or pressure, are not yet as well un- 
derstood. One feature that has slowed down progress has 
turned out to be a connection between quantum phase 
transitions and the phenomenon known as generic scale 
invariance, which refers to power-law decay of correlation 
functions in entire phases, not just at an isolated criti- 
cal point. It turns out that the critical behavior at some 
classical phase transitions is also heavily influenced by 
a coupling between generic and critical scale invariance, 
but these phenomena have usually not been cast in this 
language. The goal of this review article is to give a uni- 
fying discussion of the coupling between critical behavior 
and generic scale invariance, for both classical and quan- 
tum phase transitions, and for both static and dynamic 
critical behavior. Accordingly, we limit our discussion to 
phase transitions where this coupling is known or sus- 
pected to be important. 

The structure of this paper is as follows. In Sec. HI 
we give a brief introduction to phase transitions, both 
classical and quantum. In Sec. |n]we discuss the concept 
of generic scale invariance, and illustrate it by means of 
various examples. These concepts are less well known 
than those pertaining to phase transitions, so our expo- 
sition is more elaborate. In Sees. !III! and llV! we return 
to phase transition physics and discuss the influence of 



A. Landau Theory 

Landau theory is based on one crucial assumption, 
namely, that the free energy F is an analytic function 
of the order parameter m,^ and hence can be expanded 
in a power series. 



+ 0{m^) . (1.1) 



^ We will use the term "disordered phase" in the sense of "phase 
without long-range order". This is not to be confused with the 
presence of quenched disorder, which we will also discuss. 

^ In this section we consider a scalar order parameter for simplic- 
ity, but later we will encounter more general cases. 
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FIG. 1 Phase diagram of UGe2. After ISaxena et al\ \200\ 



Here r, u, etc., are parameters of the Landau theory 
that depend on all of the degrees of freedom other than 
m. Fi^ is sometimes referred to as the Landau functional, 
although it actually is just a function of the variable to. 
The physical value of to is the one that minimizes Fj^ . 

Landau theory is remarkably versatile. For sufficiently 
large r, the minimum of F is always located at m — 0, 
while for sufficiently small r it is located at some to ^ 0. 
If t; 7^ 0, the transition from m = to m 7^ oc- 
curs discontinuously, and the theory describes a "first- 
order transition" , with the liquid-gas transition, except 
at the liquid-gas critical point, being the prime exam- 
ple. If w = 0, either accidentally or for symmetry rea- 
sons, it describes a second-order transition, or critical 
point, at r = 0, provided u > 0. Prime examples are the 
ferromagnetic transition in zero magnetic field, and the 
liquid-gas transition at the critical point. Furthermore, 
Landau theory applies to both classical and quantum sys- 
tems, including systems at zero temperature (T = 0). In 
the latter case, the free energy F = U — TS reduces to 
the internal energy U.^ This is not an academic case. 
Consider a ferromagnet with a low Curie temperature 
Tc, e.g., UGc2. By varying a non-thermal control pa- 
rameter, e.g., hydrostatic pressure, one can suppress the 
Curie temperature to zero, see Fig. ^ The phase tran- 
sition can then be triggered in particular by varying the 
pressure at T = 0. Obviously, Landau theory predicts 
the same critical behavior in either case: The magneti- 
zation TO vanishes as m{t — > 0) = ^—r/2u, irrespective 
of how r is driven to zero. This is an example of the 
'super- universality' inherent in Landau theory: It pre- 
dicts the critical exponents to be the same for all critical 
points. For instance, the critical exponent /3, defined by 



TO cx \r\^ , is predicted to have the value /? = 1/2. 

A complete description of a statistical mechanics sys- 
tem requires, in addition to the thermodynamic proper- 
ties encoded in the free energy, knowledge of time cor- 
relation functions. Of particular interest is the order 
parameter susceptibility Xm. As a critical point is ap- 
proached from the disordered phase, the harbingers of 
the latter's instability are diverging order-parameter fluc- 
tuations, which lead to a divergence of Xm at zero fre- 
quency and wave number. Within Landau theory, these 
fluctuations are descr i bed i n a Gaussian approximation 
ijLandau and Lifshitd . Il980|) . For the wave number de- 
pendent static order parameter susceptibility Xm{k) this 
yields the familiar Ornstein-Zernike form 



(1.2a) 



with c a constant. In real space, this corresponds to 
exponential decay, 



Xm{x) OC \x\ ^ e 



(1.2b) 



where the correlation length ^ diverges, for r 0, ac- 
cording to ^ cx r~^/^. 

The order-parameter fluctuations at criticality, r = 0, 
are an example of a "soft mode",'* i.e., fluctuations that 
diverge in the limit of small frequencies and wave num- 
bers as illustrated by Eqs. H1.2|l . The soft-mode concept 
will play a crucial role in the remainder of this article. 
In addition to modes that are soft only at the critical 
point, we will discuss, in Sec. ^ modes that are soft in 
an entire phase. The coupling between such "generic soft 
modes" and the critical order-parameter fluctuations can 
have profound consequences for the critical behavior, as 
we will see in Sees. IIIll and II VI These consequences are 
the main topic of this review, but before we can discuss 
them we need to recall some additional aspects of critical 
behavior. 

Equations 1)1. 2|l define three additional critical expo- 
nents, namely, the correlation length exponent i^, deflncd 
by ^ cx |?'|~'', and the susceptibility exponents 7 and 77, 
defined by Xm{k = 0) cx r"''', and Xm{r = 0) cx |fe|^^+''. 
Landau theory universally predicts v ~ 1/2, 7=1, and 
?7 = 0. 

Universality is actually observed in experiments, but 
it is weaker than the superuniversality predicted by Lan- 
dau theory, and the observed values of the exponents are 
in general different from what Landau theory predicts. 
For instance, all bulk Heisenber g ferromagnets have a 
(3 « 0.35 ('e.g.. lZinn-JustinLll996fl . which is significantly 
smaller than the Landau value. Similarly, all bulk Ising 
ferromagnets have a common value of /3, but that value. 



^ At T > 0, a stable phase can have a higher internal energy U 
than an unstable one, as long as its entropy S is also higher. 
At T = 0, the stable phase must represent a minimum of the 
internal energy. 



* The terms "soft", "gapless", and "massless" are used inter- 
changeably in this context, with the first and second most pop- 
ular in classical and quantum statistical mechanics, respectively, 
and the third one borrowed from particle physics. 
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/3 « 0.32, is different from the one in Heisenberg systems. 
Also, the critical exponents turn out to be different for 
systems of different dimensionality, again in contrast to 
the prediction of Landau theory. For insta nce, in two- 
dimensiona l Ising ferromagnets, /3 = 1/8 (jOnsageilll948t 
lYaneilTflfil . 

B. Landau-Ginzburg-Wilson Theory, the Renormalization 
Group, and Scaling 

The reason for the failure of Landau theory to cor- 
rectly describe the details of the critical behavior is that 
it does not adequately treat the fluctuations of the or- 
der parameter about its mean value. The deviations of 
these fluctuations from a Gaussian character in general 
are stronger for lower dimensionalities, and for order pa- 
rameters with fewer components. This explains why the 
critical behavior of Ising magnets deviates more strongly 
from the Landau value than that of Heisenberg magnets, 
and why the exponents of bulk systems are closer to the 
mean-field values than those of thin films. This obser- 
vation suggests that Landau theory might actually yield 
the correct critical behavior in systems with a sufficiently 
high dimensionality d. Indeed it turns out that in general 
there is an upper critical dimensionality, dj*", such that 
for d > d^ fluctuations are unimportant for the leading 
critical behavior and Landau theory gives the correct an- 
swer (this fol lows from the Ginzburg criterion; see, e.g., 
ICardvl Il996l Sec. 2.4). For the ferromagnetic transi- 
tion at T > 0, d^ — 4. For d < d^, fluctuations need 
to be taken into account beyond the Gaussian approxi- 
mation in order to obtain the corre ct critical behavior . 
This problem was solved by Wilson ( Wilson and Koeul 
Il974t see also, e.g.. iFishcr. 1983; Ma, 1976J. He gener- 
alized the Landau functional, Eq. Hl.l|l . by writing the 
partition function Z = as a functional integral,^ 

Z = e-^/^ = y e-^[*l , (1.3a) 

where 

SW\ ^ ^/^^ [FL(0(a;)) +c(V0(a;))2] (1.3b) 

is usually referred to as the "action". Here V is the 
system volume, and is a fluctuating field whose aver- 
age value with respect to the statistical weight exp(— S*) 
is equal to m. This Landau-Ginzburg-Wilson (LGW) 
functional is then analyzed by means of the renor- 
malization group (RG);^ the lowest order, i.e., saddle- 
point/Gaussian, approximation recovers Landau theory. 



^ We use units such that Boltzmann's constant, Planck's constant, 

and the Bohr magneton are equal to unity. 
^ Ped agogical expo s itions o f the Wilsonian RG have b e en g i ven by , 

e.g.. ICardvl JlflflfiHilFish^ 119831) i lOoldenfeldl Jl992ll : lMal l|l97fii) : 
IWilson and Kogud il974D . 



Wilson's RG takes advantage of the fact that at a critical 
point there is a diverging length scale, namely, the cor- 
relation length ^, which dominates the long-wavelength 
physics. By integrating out all fluctuations on smaller 
length scales, one can derive a succession of effective the- 
ories that describe the behavior near criticality. The RG 
made possible the derivation and proof of the behavior 
near criticality known as "scaling" , which previously had 
been observed empirically and summarized in the "scal- 
ing hypothesis" ? Due to scaling, all static critical phe- 
nomena are characterized by two independent critical ex- 
ponents, one for the reduced temperature r oc |T — Td,^ 
and one for the fleld h that is conjugate to the order pa- 
rameter. Let /i — (r, K) define the space spanned by these 
two parameters. Under RG iterations, the system moves 
away from criticality according to 

^i{h)^(rh^l'' .hhy-^ , (1.4a) 

where 6 > 1 is the RG rescaling factor, v is the correla- 
tion length exponent defined above, and j/^i is related to 
the susceptibility exponent r\ by 

2/^ = (d + 2-ry)/2. (1.4b) 

The exponents Xjv and yh illustrate the more general 
concept of scale dimensions, which determine how pa- 
rameters change under RG transformations. If p is some 
parameter with scale dimension [p ] , then p (&) = p o ^ ' > 
with po — P (^ = !)■ Accordingly, [r] = 1/;/, and [/i] — 
Uh- Parameter values fi* with the property fi*{b) = /i* 
constitute a RG fixed point. Parameters with positive, 
zero, and negative scale dimensions with respect to a par- 
ticular fixed point are called relevant, marginal, and irrel- 
evant, respectively. A critical fixed point is characterized 
by the presence of two and only two relevant parame- 
ters, r and h. The critical behavior of all thermody- 
namic quantities follows from a generalized homogeneity 
law satisfied by the free energy density / = — {T/V) In Z, 

f{r,h)^h-''f{rh^'\hhy-). (1.5) 

The derivation of this relation was one of the major tri- 
umphs of the RG. In conjunction with the Wilson-Fisher 
e-expansion llWilson and Fisheilll972l) . it allows for the 
computation of the critical exponents in an asymptotic 
series about d^ .'^ 



The state of affairs just be fore the invention of th e RG was 
reviewed by K adan o ff et al . (. 19671) and lStanlevI <197lD . 

® One needs to distinguish between the exact or fully renormal- 
ized value of r, which appears in formal scaling arguments, the 
bare value of r, which appears in Landau theory or in the LGW 
functional, and any partially renormalized values. We will ex- 
plicitly make this distinction when doing so is essential, but will 
suppress it otherwise. 

® It is a popular misconception that the RG is useful only for 
dealing with critical phenomena. In fact, the technique is much 
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In addition to the diverging length scale set by ^, there 
is a diverging time scale t^. Its divergence is governed by 
the dynamical critical exponent z, according to oc 
and it r esults in the phenomenon of "critical s l owing 
down" (ivan Hovel 11954 1 Landau and Khalatmkovlll954i 
for a translation of the latter paper, see. iLandaiA Il965l 



p. 626), that is, the very slow relaxation towards equilib- 
rium of systems near a critical point. The critical behav- 
ior of time correlation functions C{k, fl; r, h) is given by 

C{k, n; r, h) = b'^^'Cikb, Q 6^; r b^'" , hby^). (1.6) 

Here fl denotes the frequency, which derives from a 
Fourier transform of the time dependence, xq is an 
exponent that is characteristic of the correlation func- 
tion C. This relation was fir st postulated as the "dy- 
namical s caling hvpothes is" i j Ferrell et gj] . 119671 Il968t 
iHalperiiT and Hohenbc^ Il967tl. and later also derived 
by dy namical RG techniques l)Hohenberg and Haloerinl 
Il977i) . 



C. Classical versus Quantum Phase Transitions 

In classical statistical mechanics, the dynamical critical 
exponent z is independent of the static cri tical behavior 
Pohenberg and Halper3 Il977t iMal Il976j) . The reason 
is as follows. Classically, the canonical partition function 
for a system consisting of N particles with / degrees of 
freedom is given as an integral over phase space of the 
Boltzmann factor, 

Z = ^( dpdq e-f^^^P^"^ 



1 

m 



dp e-Z'^-nW J dq e-'3^-'(«) 
= const. / dq e-'3^p°'('?). (1.7) 



Here P — 1/T is the inverse temperature, and p and q 
represent the generalized momenta and coordinates, re- 
spectively, H is the Hamiltonian, and i/kin and Hpot are 



more powerful and versatile, and can describe entire phases as 
well as the transitions between them, see Sec. [H] below. Al- 
though th i s was realiz ed by the foun ding fathers of the RG (see, 
lAndersont 11984 ch. 5. iFisheii [19981) . it has been exploited only 
recently (e.g., IShankaiJ. 11993). It is interesting to note that the 
systematic application of the Wilsonian RG to condensed matter 
systems implements a program that in a well-defined way is the 
opposite of that in high-energy physics. In condensed matter 
physics, the microscopic theory is known for all practical pur- 
poses, viz., the many-body Schrodinger equation. By applying 
the RG one derives effective theories valid at lower and lower en- 
ergy scales. In high-energy physics, some effective theory valid at 
relatively low energies is known (say, the standard model), and 
the goal is to deduce a 'more microscopic' theory valid at higher 
energies. 



the kinetic and potential energy, respectively.^" Due to 
the factorization of the phase-space integral indicated in 
Eq. (|1.7|l . one can integrate over the momenta and solve 
for the thermodynamic critical behavior without refer- 
ence to the dynamics. 

In quantum statistical mechanics, the situation is dif- 
ferent. The Hamiltonian H, and its constituents i/kin 
and -ffpot, are now operators, and i/kin and -ffpot do not 
commute. Consequently, the grand canonical partition 
function 



(1.8) 



does not factorize, and one must solve for the dynami- 
cal critical behavior together with the thermodynamics. 
This becomes even more obvious if o ne rewrites the parti- 
tion function as a function al integral ijCasher et a^J . ll968t 
iNeeele and OrlandL Il988|) . We will consider fermionic 
systems, in which case the latter is taken with respect 
to Grassmann- valued (i.e., anti-commuting) fields tp and 



Z= / D[^,^] 



The action S is determined by the Hamiltonian, 



(1.9a) 



J dr J dx'^ipcr{x,T) {-dr 



dr / dx H {ipcrix,T),'tp„{x,T)) , (1.9b) 



Here x denotes the position in real space, a is the spin 
index, and the fields i^aix, r) and ipa{x, t) are, for each 
value of r, in one-to-one correspondence with the creation 
and annihilation operators of second quantization, ijj^{x) 
and ipa-ix), respectively. They obey antiperiodic bound- 
ary conditions, ^paix, r = 0) = —ipcrix, t — j3). The func- 
tion H is defined such that H = J dx H{'^'^{x),'i(ja{x)). 

jjL is the chemical potential, and N in Eq. (|1.8|l is the 
particle number operator. A quantum mechanical gen- 
eralization of the LGW functional can be derived from 



Eq. I|1.9bl) by constraining appropriate linear combina- 
tions of products of fermion fields to a bosonic order 
parameter field 0, and integrating out all other degrees 
of freedom. Often one can also just write down a LGW 
functional based on symmetries and other general con- 
siderations. 

Due to the coupling of statics and dynamics, the scal- 
ing rela tion (11.51) fo r th e free energy must be generalized 
to (e.g.. [Sachdevtll999|) 

f{r,h,T) = b-^'^+''\f{rb^/'',hby'',Tb'). (1.10) 



We assume that there are no velocity dependent potentials. 
For a thorough treatment of Grassmannian algebra and calculus. 



see. iBerezinI 119661) . 

a may also comprise other quantum numbers, e.j 
depending on the model considered. 



a band index, 
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This relation reflects the fact that the temperature is 
necessarily a relevant operator at a T = critical point, 
and that temperature and frequency are expected to scale 
the same way.^"^ 

Equation ljl.9b|) displays another remarkable property 
of quantum statistical mechanics. The auxiliary vari- 
able T, usually referred to as imaginary time, acts effec- 
tively as an extra dimension. For nonzero temperature, 
(3 = 1/T < oo, this extra dimension extends only over 
a finite interval. If one is sufficiently close to critical- 
ity such that the condition > 1/Tc is fulfilled, the 
extra dimension will not affect the leading critical be- 
havior. Rather, it will only lead to corrections t o scaling 
that are dete rmined by finite-size scaling effects l|Barbeil 
ll983HCardvL ri9961. Wc thus conclude that the asymp- 
totic critical behavior at any transition with a nonzero 
critical temperature is purely classical. However, a tran- 
sition at T = is described by a theory in an effectively 
different dimension, and will therefore in general be in a 
different universality class. This raises the question how 
continuity is ensured as one moves to T = along the 
phase separation line. The resolution is that, at low tem- 
peratures, the critical region, i.e. the region around the 
phase separation line where critical behavior can be ob- 
served, is divided into several regimes. Asymptotically 
close to the transition the critical behavior is classical 
at any nonzero temperature,^"* but since this asymptotic 
classical regime is bounded by a crossover at l/r^ w T, 
it shrinks to zero as T ^ 0. In the vicinity of the quan- 
tum critical point, quantum critical behavior is observed 
except in the immediate vicinity of the phase boundary. 
This quantum critical regime is in turn divided into a 
region characterized by T < r"^, where one sees static 
quantum critical behavior approximately independent of 
the temperature, and a region characterized hy T > r"^, 
where one observes dynamic or temperature scaling. The 
latter is often called "quantum critical regime" in the lit- 
erature. Finally, there is a regime inside the critical re- 
gion, but outside both the asymptotic classical and quan- 
tum scaling regimes. This is characterized by crossover 
scaling governed by both the classical and quantum fixed 
points. These various regimes are shown schematically in 

Fig.m 



Equation (|1.9b|l suggests that imaginary time or in- 
verse temperature always scale like a length, i.e., that the 
dynamical critical exponent z — 1 at any quantum phase 
transition. Indeed, early work generalizing the Wilsonian 




QCP 



FIG. 2 Schematic phase diagram in the T-p plane, as in Fig. 
The solid line is the phase separation line, separating the 
ferromagnetic (FM) phase from the paramagnetic (PM) one, 
which ends in the quantum critical point (QCP). The critical 
region in the vicinity of the phase transition line is bounded 
by the dashed lines, which are defined as the loci of points 
where the correlation length is a certain multiple of the mi- 
croscopic length scale. The critical region is separated into 
four regimes that show different scaling behaviors: Region I 
(light shaded) displays classical scaling governed by the clas- 
sical fixed point. Regions Ila and lib (dark shaded) display 
static or temperature-independent quantum scaling governed 
by the zero-temperature fixed point, and region III (medium 
shaded) displays dynamic or temperature quantum scaling. 
In regions IVa and IVb one has crossover scaling governed by 
both fixed points. The edges of the critical region, as well 
as the boundaries between the various scaling regimes, are 
not sharp. Notice that the dashed lines are not parallel to 
the phase separation line, since both the exponents v and the 
critical amplitudes associated with the classical and quantum 
fixed points, respectively, in general have different values. 



RG to quant urn syst ems either dealt with systems where 
z = 1 (Suzuki, 1976), or assumed that z — 1 universally 
l)Beal-Monod. 1974 ). This is misleading, however. The 
point is that Eq. Ijl.9b|l represents the bare microscopic 
action. Under renormalization, or even when integrat- 
ing out degrees of freedom to derive a bare effective ac- 
tion, the relation between space and imaginary time can 
change. As a res ult, z can as sume any positive value, 
as was realized bv iHertj l)l976j) .^^ Hertz also noted the 



As we will see later, the latter expectation can be violated, due 
to, (1) the existence of multiple temperature and/or frequency 
scales , and (2) the existence of dangerous irrelevant variables 
jFisher. 1983). 

This includes the transitions at nonzero temperature in, say, 
superconductors or superfluids. Although the occurrence of the 
transition in these cases, and indeed the very existence of the or- 
der parameter, depend on quantum mechanics, the critical prop- 
erties are described by classical physics. 



One might object that the underlying microscopic theory must 
be Lorentz invariant, so ^ = 1 after all. Again, this argument 
ignores the fact that the RG derives an effective low-energy the- 
ory valid at small wave numbers and frequencies. In the critical 
theory, the relevant frequency, l/rj, is zero. Another way to say 
this is that the speed of light has been renormalized to c = oo. 
The critical theory is therefore rigorously nonrelativistic. As we 
have seen above, at any nonzero temperature h renormalizes to 
zero, and the critical theory is also rigorously non-quantum me- 
chanical. 
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following interesting consequence of the space-time na- 
ture of quantum statistical mechanics, which holds in- 
dependently of the exact value of z. Since the quan- 
tum mechanical theory is effectively higher dimensional 
than the corresponding classical one, fluctuation effects 
are weaker, and mean-field theory will be more stable. In- 
deed, since Landau theory is valid classically for d > d^, 
quantum mechanically it should be valid for d > d^ — z. 
Since d^ — 4 for many phase transitions, and usually 
z > 1, it seems to follow that most quantum phase tran- 
sitions in the most interesti ng dirnension , d = 3, show 
mean-field critical behavior. iHertj l)l97(f ) demonstrated 
this conclusion by means of a detailed LGW theory of the 
quantum phase transition in itinerant Heisenberg ferro- 
magnets, where d^ = 4 and z = 3 in Hertz's theory. 



D. The Soft-Mode Paradigm 

The above chain of arguments implies that a quantum 
phase transition is closely related to the corresponding 
classical phase transition in a higher dimension. Since 
this higher dimension will usually be above the upper 
critical dimension, the conclusion seems to be that quan- 
tum phase transitions generically show mean-field critical 
behavior, and are thus uninteresting from a critical phe- 
nomena point of view. In this review we discuss one of the 
reasons why this is in general not correct. More generally, 
we will discuss how the entire concept of a LGW theory 
in terms of a single order-parameter field can break down 
for either classical or quantum transitions, and why this 
breakdown is more common for the latter than for the 
former. The salient point is that a correct description 
of a phase transition, i.e., of phenomena at long length 
and time scales, must take into account all soft modes. If 
there are soft modes in addition to the order-parameter 
fluctuations, and if the former couple sufhciently strongly 
to the latter, then the behavior cannot be described in 
terms of a local field theory for the order-parameter fluc- 
tuations only. Away from the quantum phase transition, 
these extra soft modes, via interactions or nonlinear cou- 
plings, lead to power-law correlations in various physical 
correlation functions. This phenomenon is called generic 
scale invariance and we will discuss it in Sec. ^ 

In a derivation of a LGW functional from a microscopic 
theory this manifests itself as follows. As mentioned 
above, such a derivation involves the integrating-out of all 
degrees of freedom other than the order-parameter field. 
The parameters in the LGW functional are given in terms 
of integrals over the correlation functions of these degrees 
of freedom. If there are soft modes other than the order- 
parameter fluctuations, then these integrals may not ex- 
ist, which can lead to diverging coefficients in the LGW 
theory. At the level of Landau or mean-field theory this 
manifests itself as a nonanalytic dependence of the free 
energy on the mean-field order parameter. We will see 
later, in Secs. lIII.A.3l and lIV.A.3l respectively, that both 
in liquid crystals and in quantum ferromagnets, extra soft 



modes lead to an effective Landau free energy 

F = rni^ + vm'^lnm'^ + um'^ + . . . (l-H) 

as opposed to Eq. Hl.l|) . The situation gets worse if order 
parameter fluctuations are taken into account. Writing 
down a LGW theory based on symmetry considerations, 
or even deriving one in some crude approximation, will 
in general lead to erroneous conclusions if additional soft 
modes are present. A correct derivation, on the other 
hand, will lead to coefficients that are singular functions 
of space and imaginary time, which renders the theory 
useless. A more general approach, which keeps all of the 
soft modes on equal footing, is thus called for. Getting 
in trouble by integrating out, explicitly or implicitly, ad- 
ditional excitations is not unprecedented in physics. For 
instance, the Fermi theory of weak interactions can be 
considered as resulting from integrating out the W gauge 
bosons in the standard model. Thi s results in a th eory 
that is not renormalizable (see, e.g., lLeBellacLll99l . 

Before we enter into details, let us elaborate somewhat 
on the general aspects of this breakdown of LGW the- 
ory. Imagine a transition from a phase with an already- 
broken continuous symmetry to one with an additional 
broken symmetry. Then the Goldstone modes of the for- 
mer will in general influence the transition. A good classi- 
cal example for this mechani sm is the nematic-smectic-A 
transition in liquid crystals ('PcGcnn es and Prostl Il993|) 
which we will discuss in Sec. lIII..^ However, for the most 
obvious classical phase transitions, e.g., for the liquid- 
gas critical point, the effects of generic soft modes affect 
only the dy namics, i.e., the critical behavior of transport 
coefficients (jHohenberg and HalperinL Il977|) . Quantum 
mechanically, the concept is more important, for two rea- 
sons: (1) There are more soft modes at T = than at 
T > 0. An example that will be important for our discus- 
sion are the particle-hole excitations in an electron fluid 
that are soft at T = and acquire a mass proportional to 
T at T > 0. They couple strongly to the magnetization, 
with consequences for the quantum ferromagnetic tran- 
sitions that will be discussed in Sec. If VI below. They also 
provide a good example of the phenomenon of generic 
scale invariance, which we will discuss in detail in Sec.lTII 
(2) Quantum mechanically, the statics and the dynamics 
are coupled, as we have discussed above, and therefore 
effects that classically would affect the dynamics only in- 
fluence the static critical behavior as well. As a result 
of these two points, the original concept of mapping a 
quantum phase transition onto the corresponding clas- 
sical transition in a higher dimensionality turns out to 
be mistaken. Rather, the mapping is in general onto a 
classical transition with additional soft modes. 

In systems with generic scale invariance the correlation 
length, which diverges at the critical point, must be de- 
fined differently from Eq. Ijl.2b|) . The reason is that the 
order parameter correlations decay as power laws both 
at and away from the critical point, albeit with different 
powers, the decay being slower at criticality. In this case, 
^ is the length scale that separates the generic power-law 
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correlations from the critical ones. Specifically, in the 
scaling region, and for |a;| < ^, the correlations obey the 
critical power law, while for |a;| > ^ the correlations show 
the generic power law. 

In addition to this mechanism, there are other rea- 
sons why quantum phase transitions can be more com- 
plicated than one might naively expect. In low- 
dimensional systems, topological order can lead to 
very nontrivial effects. O ne exam p le is t he int ricate 
behavior of spin chains ijHaldaneL Il982l Il98,'^ . an- 
other, the recent proposals of exotic quantum crit- 
ical beha vior in two-d i mensiona l quantum antiferro- 
magnets llSenthil et all l2003albD . Also, the quan- 
tum phase transition may not have any classical coun- 
terpart, and hence no classical upper critical dimen- 
sion. Examples include, metal-insulator transitions 
ijEelit^^and ^KiA patrickL 119941: iKramer and MacKinnonL 
ll993';'MottVl99(T), and transitions between quantum Hall 
states ( Sondh i et a l., 1990). These phenomena are out- 
side the scope of the present review. 

Finally, we mention that there are also examples of 
quantum phase transitions where non e of the ab ove 
complications arise, and Hertz theory ijHertd . Il976fl is 
valid. This class of transitions include, quantum meta- 
magnetic transitions (Millis et ai, 2002l|, see Fig. IT^ 
for an example), and the liquid-crystaHike transitions 
in high-T^^ superconductors and quantum Hall s y stems 
llDu et all 119991: [Emerv et all Il999t iLillv et all Il999t 
lOganesvanerTI] . 2001*). One would also expect it to in- 
clude clean antiferromagnetic systems without local mo- 
ments, although no convincing experimental examples 
have been found so far. A general classification of which 
quantum phase transition should b e affected by generi c 
soft modes has been attempted bv iBelitz" et all l|2002|) . 
Although they are very interesting, we will not further 
consider these 'simple' quantum phase transitions in this 
review. 



II. GENERIC SCALE INVARIANCE 

Correlation functions characterize how fluctuations at 
one space-time point ar e correlated wi th fluctuations at 
another point (see, e.g.. iForsteil Il975|) . The default ex- 
pectation is that they decay exponentially for large dis- 
tances in space or time; fluctuations that are far apart we 
expect to be weakly correlated. If they do, then there is 
a characteristic length or time scale associated with the 
decay, and the correlations are said to be of short range. 
It is well known that at special points in the phase dia- 
gram, at critical points in particular, certain correlations 
may decay only as power laws in space or/and time (e.g., 
|Ma, 1976). In this case, the correlations are called long- 
ranged (in space) or long-lived (in time); we will call 
them long-ranged irrespective of whether the reference 
is to space or to time. Such correlations functions ex- 
hibit scale invariance: A scale change in space and/or 
time can be compensated by multiplying the correlation 



function by a simple scale factor. This property is con- 
veniently expre ssed in terms of generali z ed ho mogeneous 
functions (e.g.. IChaikin and Lubenskvl Il995l) . see, e.g., 
Eqs. H1.5II1.6|I . In addition, systems can exhibit scale in- 
variance in whole regions of the phase diagram, rather 
than only at special points. In this case one speaks 
of "generic scale invariance" (GSI)fDorfman. et al ]. 11994 
Eaw and Nieuw oudt. 1989; Nagel, 1992). In this section 
we discuss various mechanisms for GSI in both classi- 
cal and quantum systems, highlighting their similarities 
and differences. Before we go into details, we list four 
mechanisms that can lead to this phenomenon. 

(1) Spontaneous breaking of a global continuous sym- 

metry leads to Goldstone modes that are massless 
everywhere in the broken-symm etry phase (e.g., 
IForsteil IT973 IZinn- JustinL Il996j) . This mechanism 
is operative in both classical and quantum systems. 

(2) Gauge symmetries lead to soft modes since a 

mass would be incompatib le with ga uge invari- 
ance (iRvdeillTgsllWeinbera. .1996a.b.: , ■Zinn-Justinl 
Il996j) . The only case we will be interested in is the 
U(l) gauge symmetry that underlies the massless- 
ness of the photon. 

Both of these mechanisms, if operative, lead to soft 
modes that directly result in obvious ways from the sym- 
metry in question, and hence to GSI that we will refer 
to as 'direct' GSI effects. What is usually less obvious is 
that these soft modes, via mode-mode coupling effects, 
can lead to other modes becoming soft as well. We will 
see various examples for these 'indirect' GSI effects. An- 
other source of direct GSI are 

(3) conservation laws, which lead to power-law t emporal 

decay of local time correlation functions ijForsteil 
Il975fl . There are two ways in which conservation 
laws can lead to indirect GSI, namely: 

(3a) The conservation laws can appear in con- 
junction with mode-mode coupling effects, 
or nonlinearitie s in the equations of motion 
jBoon a nd Yhl Il991t |P omeau and Resiboia . 
[l975). In classical systems, this leads to long- 
ranged time correlation functions; in quantum 
systems, it leads to long-ranged time correla- 
tion functions and thermodynamic quantities. 

(3b) Conservation laws in co njunction with a 
nonequilibrium situation 



Kirkoatrick et all 
ISchmittmann and Zia» .199511 can 
leads to long-ranged time correlation functions 
and thermodynamic quantities even in classi- 
cal systems. 

In what follows we discuss specific examples, starting 
with classical systems. 
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A. Classical Systems 

We now illustrate the four mechanisms listed above by 
means of four examples that show how they lead to long- 
ranged correlations in classical systems. For each case we 
give a general discussion, followed by the specification of 
a suitable model, and a demonstration of how explicit 
calculations and, where applicable, RG arguments, yield 
GSI. 



1. Goldstone modes 

To illustrate the Goldstone mechanism for long-ranged 
spatial correlations in equilibrium, we choose as an ex- 
ample a generalized Heisenberg ferromagnet with an TV- 
dimensional order parameter. A'^ = 3 represents the 
usual Heisenberg magnet, which will be relevant later in 
this review. Let (p^x) be the fluctuating magnetization, 
with m = {(Pix)) its average value, and let h be an ex- 
ternal magnetic field conjugate to the order parameter . 
Goldstone's theorem l)Goldstonel Il96ll: iGoldstone et all 
^962) states that whenever there is a spontaneously bro- 
ken global continuous symmetry, there will be massless 
modes. (This is not true for local, or gauge, symmetries, 
see Sec. III.A"2l below.l To understand this concept, sup- 
pose that the action is invariant under rotations of the 
vector field </> provided that h — 0. This implies that the 
order parameter in zero field vanishes, Tn{h = 0) = 0, 
due to rotational invariance. However, in the limit — s- 
there are two possible behaviors, depending on the tem- 
perature, which determines whether the system is in the 
disordered phase, or in the ordered phase. Namely, 

Tn{h ^ 0) = (disordered phase), (2.1a) 

m{h ^ 0) 7^ (ordered phase). (2.1b) 

The situation in the ordered phase is sometimes also char- 
acterized by saying that the action obeys the symmetry 
while the state does not, and it is referred to as a spon- 
taneously broken continuous symmetry. In the ordered 
phase, one still has invariance with respect to rotations 
that leave the vector m fixed, i.e., with respect to the 
subgroup 0{N — 1) that is the little group of m. Gold- 
stone's theorem says that this results in as many soft 
modes as the quotient space 0{N)/0{N — 1) has di- 
mensions, i.e., there are dim {0{N)/0{N — 1)) = TV — 1 
Goldstone modes. These soft modes are 'perpendicular' 
or 'transverse' to the direction m of the spontaneous or- 
dering. For example, if </> = ((f>i, . . . ,(t>N) and there is 
spontaneous ordering in the TV-direction, then one has. 



in the limit h ^ and for asymptotically small wave 
numbers, 

(0,(fc)0,(-fc)> =%/Cfe2 , (z,j = l,...,7V-l). 

(2.2a) 

C is called the stiffness coefficient of the Goldstone modes. 
In real space in, say, three-dimensions, this implies that 
for large distances, \xi — X2 \ ^ oo. 



i{xi) (f>j{x2)) oc l/|a;i - X2\. 



(2.2b) 



We see that Goldstone's theorem gives rise to power-law 
correlations, or GSI, in the entire ordered phase. In the 
terminology explained above, this is an example of direct 
GSI. 



a. Nonlinear a model To illustrate the Goldstone 
mechanism more explicitly, and for later reference, 
let us consider the deri va tion o f a nonlinear a mode l 
llBrezin and Zinn- JustinL Il976l: ICell - Mann and Lew 



1 196(1 iNelson and PelcovitsL 119771: iPolvakovl 1197. 
lzinmJustinr ^9^ for our Ol'iWsvmmetric Heisenberg 
ferromagnet. We specify the action by 

S[ct>] = J dx [rcl}^{x)+c {Vct>{x)f + u(l}^{x) 

-h<j)N{x)]. (2.3) 

Here </)2 = (P,^, and (Vc/))^ = da(j),d°'(f)\ with i = 
1, . . . N, and a = 1, ... c? in d-dimensions. Summation 
over repeated indices is implied, and we have assumed 
an external field in A^-direction. S determines the par- 
tition function via Eq. H1.3a|l . For h = 0, it is obviously 
invariant under 0(A^)-rotations of the vector field 0. In 
the low-temperature phase, where the 0(A^)-symmetry 
is spontaneously broken, it is convenient to decompose 
<p into its modulus p and a unit vector field 0, 



4>{x) ^ p{x) 4){x) , 4)^{x) = l. 



(2.4) 



4> parameterizes the (N — l)-sphere which is isomorphic 
to the coset space 0{N) / 0{N — 1). In terms of p and </>, 
the action is^* 

S[p,(i)\ ^ c(i) j dx p^{x)[V^{x)]'^ 

dx rp'^{x)+c^'^\Vp{x)f +up'^{x) 

(2.5) 



h I dx p{x) 4)m{x). 



In this case, the continuous symmetry is characterized by the 
group 0{N), but the concept is valid for arbitrary Lie g roups. A 
useful reference for Lie groups, or continuous groups, is lGilmora 
IT97I . 



Provide d that d > 2; no sponta neous symmetry breaking occurs 
in d < 2 tMermin and Wagneii llQeet. 

The change of variables from </> to (p, </>) also changes the inte- 
gration measure in Eq. lOal . If one elimi nates a in terms of tt, 
Eqs. 12.61 . the measure changes again, see IZinn- JustinI ^^9^. 
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The bare values of c^^^ and c*^^) are equal to c. Notice that 
the field 4> appears only in conjunction with two gradient 
operators. This implies that the </>- fluctuations represent 
the N — 1 Goldstone modes of the problem, while the p- 
fluctuations represent the massive mode. Assuming the 
system is in the ordered phase, and taking the order to 
be in the TV-direction, we parameterize cj) &s follows, 



0(a;) = {7v{x),a{x)) , 



where the vector tt represents the N — 1 transverse di- 
rections, and 



a{x) = [l - n^{x)] 



1/2 



(2.6b) 



Next we split off the expectation value of the massive 
field by writing p{x) = M + dp{x), with M ~ {p{x)). 
Absorbing appropriate powers of M into the coupling 
constants, we can write the action 



S[p, tt] ^ S'nlctmM -I- dS[p, tt]. (2.7a) 



Here 



'- J dx [iV7rix)r + i^aix))'] 
J dx a{x), (2.7b) 



is the action of the 0(A^)-nonlinear a model, with the 
bare value of ( equal to 2cM^. 6S is the part of the 
action that contains the 5yO -fluctuations, as well as the 
coupling between Sp and 



6S 



dx iSp{x)y + O Sp-', {VSpy,aSp, Sp{V<jyy 



(2.7c) 

This parameterization of the model accomplishes an ex- 
plicit separation into modes that are soft in the broken- 
symmetry phase (i.e., the Tr-fields), modes that are mas- 
sive (i.e., the -field), and couplings between the two. 



b. GSI from explicit calculations At the Gaussian level, 
the action reads 

So ^ ^Jdx {Vn{x)f + \jdx {^{x)f 

+ r j dx (Jp(a;))2 +c(2) j dx {V5p{x)f. (2.8) 

Here we see explicitly that the Gaussian 7r-propagators 
are soft as /i ^ 0, 



(^,;(fc)7r,(-fc)) =5,,/(Cfe2 + /l). 



(2.9) 



That is, they have the form given by Eq. H2.2a|l . while the 
(5/0 -pr opagator is massive. Explicit perturbative calcula- 
tions ijNelson and Pelcovitsl Il977(l show that this does 
not change within the framework of a loop expansion. It 



is also interesting to consider the explicit perturbation 
theory for the magnetization m = M{a{x)). To one-loop 
order, Eq. Ij2.6b(l yields 

m/M = ^-\J {Mp)^ji-P)) 

= m{h = 0)/M + const, x /i^^-^V^. (2.10) 



(2.6a) Here / = Jdp/{2nY, and we show only the leading 



nonanalytic dependence of m on h. 



c. Generic scale invariance from RG arguments The next 
question is whether these perturbative results are 
generic, and valid independe nt of perturbation the- 
ory. Withi n a pertur bative RG f Brezin and Zinn-Justinl 
il976: Nelson and Pelcovits . 1977n ) this can be checked or- 
der by order in a loop expansion. The fact that Eqs. H2.9I 
I2.10|l are valid to all orders in perturbation theory, and 
are indeed exact properties independent of any perturba- 
tive scheme, hinges on th e proof of the renorm alizability 
of the nonlinear a model (jBrezin et all Il976j) . 

Here we present a much simple r, if less complete, argu- 
ment based on power counting (Belitz a nd Kirkpatrickl 
Il997j) . We employ the concept of .Ma (,,1976^ ■ wherebv 
one postulates a fixed point, and then self-consistently 
checks its stability. Accordingly, we assign a scale di- 
mension [L] = —1 to lengths L, and look for a fixed 
point where the fields have scale dimensions 



[7r,{x)] = (d-2)/2, 
[Spix]] = d/2, 



(2.11a) 
(2.11b) 



in d-dimensions. This ansatz is motivated by the ex- 
pectation that the Gaussian approximation for the tt- 
correlation, Eq. H2.9I) . is indeed exact, and that 5p is mas- 
sive. Power counting shows that C, and r are marginal, 
while all other coupling constants are irrelevant except 
for h, which is relevant with [h] = 2. In particular, all 
terms in 5S that do not depend on h are irrelevant with 
respect to the putative fixed point, and so is the term 
(Vcr)^ in S'nlctM- The fixed point is thus stable and de- 
scribes the ordered phase. The fixed-point action plus 
the most relevant external-field term is Gaussian, 

5fp - ^ y dx {Vn{x)f + \jdx {TT{x)f 

+ r [ dx {Sp{x) f. (2.12) 



We see that the tt-tt correlation function for asymptot- 
ically small wave numbers is indeed given by Eq. H2.9|l . 
i.e., one has soft Goldstone modes, or GSI, everywhere 
in the ordered phase. Notice that the above RG argu- 
ments prove this statement, albeit the proof is not rig- 
orous in a mathematical sense. Furthermore, the fixed 
point value of the magnetization, m = M{a{x)), is given 
by M, and the leading correction is given by the tt-tt 
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correlation function. The scale dimensions of tt and h 
then yield 



and a new action 



i{h) = m{h 0) + const, x /i('*-2)/2 



(2.13a) 



in agreement with the perturbative result, Eq. H2.10|l . 
This in turn implies that the longitudinal susceptibility 



XL = dm/dh cx h^'^-^^'^ 



(2.13b) 



diverges in the lim it /i — > for all d < 4 
l)Brezin and Wallacel fl973). Alternatively, the zero-field 
susceptibility diverges in the homogeneous limit as 



xL(fc^o)cx |fcr(4-<i). 



(2.13c) 



In real space this corresponds to a decay as 1 / \x\'^^'^~'^\ 
The above arguments assume that no structurally new 
terms are generated under renormalization which might 
turn out to be marginal or relevant; this is one of the 
reasons why the proof is not rigorous. With this caveat, 
they show that several properties of isotropic Heisen- 
berg ferromagnets that are readily obtained by means 
of perturbation theory are indeed exact. They also il- 
lustrate how much more information can be extracted 
from simple power counting after performing a symme- 
try analysis, and separating the soft and massive modes, 
as opposed to doing power counting on the original ac- 
tion, Eq. t,Ma, 1976). More importantly for our 
present purposes, they illustrate how soft modes in the 
presence of nonlinearities lead to slow decay of generic 
correlation functions. In the present case, the transverse 
Goldstone modes couple to the longitudinal fluctuations, 
which leads to the long-range correlations expressed in 
Eq. H2.13c|l . This is an example of indirect GSI due to 
mode-mode coupling effects. In Sec. III.A.3I we will see 
that a very similar mechanism leads to long-ranged time 
correlation functions in classical fluids. 



2. (7(1) gauge symmetry 

The second mechanism on our list is GSI caused by a 
gauge symmetry. Let us consider 0^-theory again, Eq. 
(|2.3|) . with = 2 or, equivalently, with a complex scalar 
field (j). Let us further demand that the theory is invariant 
under local U{1) gauge transformations, 



JA{x) 



(2.14a) 



with an arbi trary real field A. It is well known (e.g., 
iRvdeil 119851) that this demand forces the introduction 
of a gauge field A with components (a = 1,2,3 in 
three-dimensions) that transforms as 



A{x) 



A\x)=A{x) + -VA{x), (2.14b) 



A detailed analysis shows that SS does not contribute to the 
leading corrections to scaling at the stable fixed point. 



S[<j>,A] 



dx 



'u\^ix)\' 



r\(f>{x)f + c\[V -iqAix)](f>{x)f 
F^f,ix)F"('ix)\, (2.15a 



where 



Fa[3 (x) = da A/3{x) - dp Aa{x). 



(2.15b) 



Here fi and q are coupling constants that characterize the 
gauge field A and its coupling to cj), respectively. The 
usual interpretation of this action is that of a charged 
particle, or excitation, with charge g, described by ip, 
that couples to electromagnetic fluctuations, or photons, 
described by the vector potential A. We will see in Sec. 
nil. A. II however, that it can describe other systems as 
well, at least in certain limits. 



a. Symmetric phase Let us flrst discuss the symmetric 
phase, where both fields have zero expectation values. 
In this phase, we expect the Gaussian A-propagator to 
be massless, as the A-field appears only in conjunction 
with derivatives. This is a consequence of the local gauge 
invariance; a term of the form m?A'^{x) would violate the 
latter. However, the action S, Eq. (|2.15a|l . comes with 
the usual problems related to gauge fields. That is, the 
A-propagator does not exist since the A-vertex has a 
zero eigenvalue. We deal with this problem by "gauge 
fixing", i.e., we work in Coulomb gauge, V • A{x) = 0, 
which we enforce by ad ding to the action a gauge fixing 
term Ce.g.. lRvdeilll985i ch. 7.1) 

5gf = - I dx {V-A{x)f, (2.16) 
V J 

with 77 0. One finds from Eqs. (I^T^ITIB 



{Aa{k)Ap{^k))=A^I^l 
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fc2 



(2.17a) 



so the photon is indeed soft. Local gauge invariance thus 
leads to GSI in an obvious way; this is another exam- 
ple of direct GSI. Notice that this mechanism is distinct 
from the one related to conservation laws to be discussed 
in Sec. III. A"3l below: While local gauge invariance is suf- 
ficient for the conservation of the charge q, it is not nec- 
essary; a global U{1) suffices to make q a conserved quan- 
tity. 

The 0-ficld, on the other hand, has two massive com- 
ponents with equal masses. Writing (/> = ((/>! -I- i02)/\/2, 
with (f)i and 02 real, we have 

(</),(fc)0,(-fc)) =5y/(r + cfc2) , (*,j = l,2). 

(2.17b) 

We thus have two massive scalar fields, and one mass- 
less vector field with two degrees of freedom. 



The vector field, or photon, has only two degrees of freedom, 
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b. The broken-symmetry phase We now turn to the phase 
where the local t/(l)-symmetry is spontaneously broken. 
Suppose the spontaneous expectation value of 4> is in the 
01-direction, so we have 

(t){x) =v+ (Mx) + iMx)) /V2, 



with V real and equal to y^—r/2u at tree level. If we 
expand about the saddle-point solution (j) ~ v, A ^ 0, 
we find to Gaussian order 



S 



dx 



+ cq\^A{x)f + -—F^p{x)F'^P{x) 

lOTT/i 

+ 1 {V(j)2{x)f - V2cqvA{x) ■ VMx) 



(2.18) 



There arc two interesting aspects of this action. First, 
the vector field has acquired a mass that is proportional 
to v"^. Second, the field (j)2, which is often called the 
Higgs field in this context, and which we would expect 
to form the Goldstone mode associated with the spon- 
taneously broken U{1) or 0(2)-symmetry, couples to the 
now-massive vector field. Indeed, (j)2 can be eliminated 
from the Gaussian action by shifting A, A{x) A{x) — 
V (j)2{x) / \/2qv . Notice that this shift just amounts to a 
change of gauge, and hence does not change the physical 
nature of the vector field. It does, however, give A a 
longitudinal component and thus increases the number 
of photon degrees of freedom from two to three. With 
this shift, the Gaussian action consists of only the first 
four terms in Eq. (|2.18|) . which leads to the following 
propagators. 



4.r /"^+,'"'f/"^' , (2.19a) 



1 



(2.19b) 



with = An^cq^v^. 

We see that now there are two massive fields, namely, 
a massive scalar field with one degree of freedom, and 
a massive vector field with three degrees of freedom. In 



rather than three. In the gauge we have chosen this is obvious, 
since the propagator, Eq. I2.17al . is purely transverse. 
^ By writing </> in terms of an amplitude and a phase, one can 
choose a gauge such that the entire action, not just the Gaus- 
sian part, is independent of 4>2- This is known as the "physical" 
or "unitary" gauge (e.g.. iRvdeil IToS^ ch. 8.3). Other choices, 
which retain </>2 and lead to a different A-propagator, are pos- 
sible. It has been shown that all these formulations are indeed 
physically equivalent; the contributions from any nonzero ip2- 
propagator cancel against pieces of the A-propaga tor and thus d o 
not contribute to any observab l e prop erties. See, iRvdeJ J1985D : 
IWeinber3 ll996bD ; IZinn-.Tustinl Jl99fili for detailed discussions of 
this point. 



particular, there is no Goldstone mode, despite the spon- 
taneously broken 0(2)-symmetry. Pictorially speaking, 
the gauge field has eaten the Goldstone mode and has 
become massive in the process. This phenome non is com- 
monly refer r ed to as the Higgs mechanism l)AndersonL 
1196,'^ iHiggsL ll964alH ) . Notice that the situation is com- 
plementary, in a well-defined sense, to the case without 
a gauge field in Sec. III. A. II Without local gauge invari- 
ance, one has two massive modes in the disordered phase, 
and one massive mode and one soft Goldstone mode in 
the ordered phase. In the gauge theory, there is a mass- 
less mode (the photon) in the disordered phase, and no 
massless modes in the ordered phase. While there are 
no obvious examples of indirect GSI involving photons, 
we will see in Sec. IIII.A.ll that the nature of the mass 
acquired by (t>i in the ordered phase can have a drastic 
influence on the nature of the phase transition. 

We close this subsection with one additional remark. 
In particle physics, local gauge invariance is an indis- 
pensable requirement, since it is necessary for Lorentz 
invariance. In the theory of phase transitions there is no 
such requirement, since the critical theory is rigorously 
nonrelativistic, see footnote Nevertheless, it is im- 
portant if an order parameter field couples to photons, 
and also in some other cases. We will elaborate on this 
in SeclTTIl 



3. Long-time tails 

We now turn to the third mechanism on our list; 
namely, long-ranged correlations in time correlation func- 
tions due to conservation laws. Let us consider time 
correlation functions involving the local currents of con- 
served quantities, e.g., mass, momentum, or energy. 
Time integrals over the spatial averages of these corre- 
lation fimctwnBj_j;efe as G reen-Kubo expressions 
ijGreenL 11954 iKubol Il957l Il959() . determine the trans- 
port coefficients of fluids. A simple example is the ve- 
locity autocorrelation function of a tagged particle. It is 
deflned as 



CD{t) = {v{t)-vm 



oqj 



(2.20a) 



where v{Q) is the initial velocity of the tagged particle, 
v{t) its velocity at a later time t, and (. . .)cq denotes 
an equilibrium ensemble average. The coefficient of self 
diffusion in a d-dimen sional system is then given by (e.g., 
iBoon and Yi^ . ll99l|) 



D=^ I dt Coit). 



(2.20b) 



Analogous expressions determine the coefficients of shear 
viscosity, r], bulk viscosity, C, and heat conductivity. A, 
in a classical fluid (and, when appropriate, in a classical 
solid). In general we denote these current correlation 
functions by C^i(i), where fi can stand for D, r], C, or A. 
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FIG. 3 Normalized velocity autocorrelation function poit) — 
C oit) / {v"^ (Q)) as a function of the dimensionless time t* — 
t/to, where to is the mean- free time. The crosses indicate com- 
puter results obtained bv lWood and Er pcnbcck (1975) for a 
system of 4000 hard spheres at a reduced density correspond- 
ing to V/Vo — 3, where V is the actual volume and Vo is the 
close-packing volume. The dashed curve represents the theo- 
retical curve poit) ~ cto {t*)~^^^. The solid curve represents 
a more complete evaluation of the mode-coupling formula 
with contributions from all possible hydrodynamic modes and 
with finite-size correc tions included (iDorfman. .1981') . From 
iDorfman eTah Jl99l). 



In traditional many-body theories (e.g., the Boltzmann 
equation) the were always found to decay exponen- 
tially in time lICercignaiiilllQSSHChapnian and Cowlind . 



I1952I: iDorfman and vaii Beiierenl Il977|l . with a char- 
acteristic decay time on the order of the mean-free 
time between collisions, and until the mid-1960s this 
was believed to be generally true. It thus came 
as a great surpris e when both numeric a l mol e cular- 
dyna mics studies ijAlder and WainwrightL Il967l Il968i 
|197C ), and, shortly thereafter, more sophisticate d theo- 
ries Dorfman and Cohenl 1197(1 Il972 . 1975; Erns t et all 
1197(1 11971I ll976albD . showed that all of these correla- 
tions decay for asymptotically long times as l/f^/"^. This 
power-law decay in time is a type of GSI that is referred 
to as long-time tails (LTT). In Fig. Owe show results of 
both computer simulations and theoretical calculations 
for Coit) for a hard-sphere fluid in three-dimensions. 
The LTT is clearly visible. 



a. Fluctuating hydrodynamics The simplest and most 
general way to understand the LTT mechanism in a clas- 
sical flu id is via the equa tions of fluctuating hydrody- 
namics ( Ernst et a/.lll97l[l . i.e., the Navier-Stokes equa- 
tions with appropriate Langevin forces added, although 



explicit many-body calculations, within the framework 
of a generalized kinetic theory, are also possible and 
give identic al results for the LTTs ()Dorfman and Cohenl 
I1972L Il975|) .^^ The hydrodynamic approach gives the 
exact long-time behavior because the slowest-decaying 
fluctuations in a classical fluid are the fluctuations of the 
conserved variables, which are described by the hydrody- 
namic equations. The conserved variables are, the mass 
density p, the momentum density g = pu, with u the 
fluid velocity, and the energy density e or, alternatively, 
the entropy density s. In the long- wavelength and low- 
frequency limit the fluctuations of these variables are de- 
scrib ed exactly by the equations (jLandau and Lifshit j . 
I1987D 



dtga + dpgau^ = -daP 



?7 {dau'' 



(2.21a) 

(2.21b) 
(2.21c) 



Here p denotes the pressure, and summation over re- 
peated indices is implied. In Eq. (|2.21c() we have ne- 
glected a viscous dissipation term that represents entropy 
production since it is irrelevant to the leading LTTs. The 
Langevin forces Paf3 and qa are uncorrelated with the ini- 
tial hydrodynamic variables, and satisfy 



jd-l) \,' 

1 V Oap Opiy 



6{x-x')5{t-t'), (2.22a) 



{qc.{x,t)qp{x\t')) = 2\TH^p5{x-x')5{t-t'), 

(2.22b) 

{P^p{x,t)q^{x\t')) = 0. (2.22c) 

The above equations can be derived in a number o f ways 
(IFox and UhlenbeckL Il97(l iLandau and LifshitzL Il987^ 

and are known to exactly describe the long-wavelength 
and low-frequency fluctuations in a fluid. 



b. GSI bom explicit calculations To illustrate how the 
LTTs arise we choose a slightly different example than 
the self-diffusion coefficient discussed above, namely, the 
shear viscosity r]. The appropriate time correlation func- 
tion that enters the Green-Kubo formula in this case is 
the autocorrelatio n of the tran sverse velocity current, or 
the stress tensor (|Ernst et a,l\ . Il97i|) . The basic idea is 



Notice that the term "kinetic theory", which is often used syn- 
onymously with "Boltmann theory" , is used much more generally 
by theses authors. 
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that r] in the above equations is really a bare transport co- 
efficient that gets renormalized by the nonlinearities and 
fluctuations. We \yill pre sent a simplified calculation of 77 
l|KirkT)atrick et all EOO^ and then discuss more general 
results. 

For simplicity, we assume an incompressible fluid. The 
mass conservation law then reduces to the condition 

V-u{x,t)^0, (2.23) 

and the momentum conservation law is a closed partial 
diiTerential equation for u, 

(2.24) 

Here v ~ 77/p is the kinematic viscosity, which we assume 
to be constant. The pressure term serves only the enforce 
the condition of incompressibility. In fact, it can be elim- 
inated by taking the curl of Eq. 1)2.24(1 , which turns it into 
an equation for the transverse velocity. The cause of the 
LTTs is the coupling of slow hydrodynamic modes due 
to the nonlinear term in Eq. H2.24|l : this is another ex- 
ample of the mode-mode coupling effects mentioned in 
the introduction to the current section. Since we will 
treat this term as a perturbation, we have formally mul- 
tiplied it by a coupling constant 7 whose physical value 
is unity. We will take the nonlinearity into account to 
lowest nontrivial order in 7. 

Consider the velocity autocorrelation function tensor, 

Cap{k, t) = (w„(fc, t) up{-k, 0)). (2.25) 

An equation for C can be obtained by Fourier trans- 
forming Eq. (|2.24|) . multiplying by u/3(— fc,0), and av- 
eraging over the noise while keeping in mind that the 
noise is uncorrelated with the initial fluid velocity. In the 
case of an incompressible fluid we need to consider only 
the transverse- velocity correlation function, C_\_. This 
is easily done by multiplying with unit vectors, k!"^ 
[i — l,2,...,d — 1), that are perpendicular to fc, which 
eliminates the pressure term. We obtain 

{dt + vk^)C^{k,t) = -tlk^k^llk%Yl 

q 

X (M^(fc - q, t) u'^iq, t) u^{-k, 0)). (2.26) 

Here we have used the incompressibility condition, Eq. 
(|2.23l) . to write all gradients as external ones. To ze- 
roth order in 7 we find, with the help of the f-sum rule 
llForsterllTflTl C^{k,Q) ^T/p, 

{k, t) = (T/p) e-"^^'* + 0(7). (2.27a) 

This is the st andard result obtained from li nearized hy- 
drodynamics ijChapman and Cowhnel [7952(1 . which pre- 
dicts exponential decay for k ^ 0.^^ Notice that it 



The reader might find it curious that this derivation did not 
make explicit use of the Langevin force correlations, Eqs. I2.22i . 
However, it needs various equal-time correlation functions as in- 
put, which contain the same information as Eqs. 12.221 . 



amounts to GST in time space for the local time correla- 
tion function, 

Ci(a; = 0,<) cx l/i'^/2. (2.27b) 

This is an immediate consequence of the conservation law 
for the transverse momentum, and hence an example of 
direct GSI. 

To calculate corrections due to the nonlinearity, we 
need an equation for the three-point correlation function 
in Eq. 12.26|l . The simplest way to obtain this is to use 
the time translational invariance properties of this cor- 
relation function to put the time dependence in the last 
velocity, and then use Eq. I(2.24|l again. The result is 
an equation for the three-point function in terms of a 
four-point one that is analogous to Eq. H2.26|l . To solve 
this equation, we note that, due to the velocity being 
odd under time reversal, the equal-time three-point cor- 
relation vanishes. By means of a Laplace transform, one 
can therefore express the three-point function as a prod- 
uct (in frequency space) of the zeroth-order result for 
C^, Eq. ((2.27a|l . and the four-point function. To leading 
(i.e., zeroth) order in 7 the latter factorizes into products 
of velocity autocorrelation functions. Upon transforming 
back into time space, and to quadratic order in the cou- 
pling constant 7, we obtain 

{dt + vk^)C^{k,t)+ f dT^{k,t-T)C_L{k,T) =0, 
Jo 

(2.28a) 

with 

S(fc,i) = 72|A:^fc,^[C«^(q,t)C'^''(fc-q,t) 

+ C^^q, t) C^^ik - q, t)] k^^Mp + 0(7'). (2.28b) 

The self-energy S is proportional to fc^, and thus pro- 
vides a renormalization of the bare viscosity z/, or the 
time correlation function that determines the shear vis- 
cosity, Cri {t) . This correction to v is time and wave num- 
ber dependent. 

The source of the LTTs is now evident. In our model 
of an incompressible fluid, only the transverse component 
of Cap is nonzero, 

Cap{q,t) = {Sa(3-qaqp)C^iq,t). (2.29) 

Putting 7 = 1, and defining SCri{t) = 
plimfe^o S(fc, i)/A;^, one obtains for asymptotically 
long times 

^^"(*)-^^(;|tI))(8^- ^'-''^^ 

The correction to the static, or zero-frequency, shear vis- 
cosity is given by Sr] oc dt6Crj{t), cf. Eq. (|2.20bl) . We 
thus have 

C^{t ^ 00) ^dCriit-^ 00) (xt-'^/^. (2.30b) 
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This is the well-known contribution of the t ransverse ve- 
locity m odes to the LTT of the viscosity l)Ernst et all 
^2Sil3)- Notice that Cr, describes correlations of shear 
stress, which is not conserved. The LTT therefore is a 
result of mode-mode coupling effects, and hence an ex- 
ample of indirect GSI. In a compressible fluid, a similar 
process coupling two longitudinal modes also contributes 
to the leading LTT. The other transport coefhcients, e.g., 
A and ^, also have LTTs proportional to l/t'^^'^, and all of 
them have less-leadi ng LTTs proportional to 1 /^*^'^"'"^-'/^ 
or weaker ( Dorfma^. [l98lt lErnst et aLLIl976aibl) . 

For the frequency-dependent kinematic viscosity, the 
nonexponential decay of Svit) implies a nonanalyticity 
at zero frequency. More generally, the algebraic l/f^^^ 
LTTs in the long-time limit imply, for the frequency or 
wave number-dependent transport coefficients /i, a non- 
analyticity at zero frequency fl, or wave number |fc|,^^ 



and ST. If we neglect the nonlinearity in Eq. Ij2.21b|l . 
Eqs. (|2.21|l can be written 

dtu^^o. = i^dpd^u^,^ + — (d^P^p)^, (2.32a) 
Po 



dt ST + Uc^d^'T = Dt dc^d'^ST 



PoToCp 



daq", (2.32b) 



where po is the average mass density, Cp is the specific 
heat per mass at constant pressure, and Dt ~ A/po Cp is 
the thermal diffusivity. These bilinear equations can be 
solved by means of Fourier and Laplace transformations. 
Focusing on static, or equal-time, correlations, one finds, 
for example. 



{\5p{k)\')^pT 



(dp\ 



^^("-"^ ■ , (2.33a) 

DT{iy + DT)k^ 



(2.31a) 
(2.31b) 



where the prefactors and 6^ are positive, and only the 
leading nonanalyticities are shown. For the implications 
of these results in d < 2, see footnote below. 



4. Spatial correlations in nonequilibrium steady states 

We lastly consider a fluid in a nonequilibrium steady 
state; the equations of fluctuating hydrody namics can be 
extended to this case (iRonis et alV Il980j) . It has been 
known for some time that these systems in general ex- 
hibit GSI in both time correl ation functions and th ermo- 
dynamic quantities (see, e.g.. lDorfman et all . ll994|) . The 
spatial correlations responsible for the GSI in thermody- 
namic susceptibilities are closely related to LTTs of the 
equilibrium time correlation functions. We will consider 
a fluid in a steady, spatially uniform, temperature gradi- 
ent VT, but far from any convective instability. Further, 
we use a number of approximations that enable us to fo- 
cus on the most interesting effects of such a gradient. For 
a justification of this procedure, as well as the underly- 
ing details, we refer the reader to the original literature 
iKirkpatrick et ai. 1982a,bl . 

We write the temperature T — Tq + ST as fluctuations 
ST about an average value Tq, and focus on the coupling 
between fluctuations of the transverse fluid velocity u± 



To avoid misunderstandings, we note that a static transport co- 
efficient, fi{k,Q = 0), is a time integral over a time correlation 
function, see the Grecn-Kubo formula, Eq. 12.20bl . As such, 
it is long-ranged in space, just as fi{k = 0, Q) is in time. A 
static susceptibility, on the other hand, is related to an equal- 
time c orrelation function via the fluctuation-dissipation theorem 
iFors tcr. 1975) and is not long-ranged in space in classical equi- 
librium systems. 



([fei • g(fc)] Sp{-k)) = pTaT 



k^-VT 

J^TTDtW 



where 



p \dT 



(2.33b) 



(2.33c) 



is the thermal expansion coefficient at constant pressure. 

There are several remarkable aspects of these results. 
First, Eq. H2.33a|l for the density correlations implies that 
the first term, which also exists in equilibrium, is delta- 
correlated in real space, while the secon d term decays 
as const. — \x\ in three-dime nsions l)Schmi tz and CohenL 
Il985t Ide Zarate et all l200l[l . Equation H2.33bp shows 
that the transverse-momentum-density correlation func- 
tion decays as l/|a;| in three-dimensions. Both of 
these results show that spatial correlations in a nonequi- 
librium steady state exhibit GSI. Second, the right- 
hand side of Eq. (|2.33b|) is essentially the integrand 
of a LTT con tribution to th e heat conductivity A 
l|Hohenberg and Halperiil ll977|) . This demonstrates the 
close connection between the LTTs in equilibrium time 
correlation functions and the spatial GSI of equal-time 
correlation functions in nonequilibrium situations. As in 
the case of the equilibrium LTTs, this is an example of 
indirect GSI. 

These fluctuations can be directly measured by small- 
angle light-scattering experiments. Specifically, the dy- 
namic structure factor Spp{k,t) — {Sp{k,t)Sp{—k,0)) , 
which is proportional to the scattering cross sec- 
tion, in a nonequilib rium steady state has the form 
l|Dorfman et all \1994 



Spp{k,t) = So {1 + At 



e 



-A, 



(2.34a) 



where 5o is the structure factor in equilibrium, and 
Cp V {k± ■ VT f 



An 



TDT{iy^ - Dl) 



fc4 



(2.34b) 



16 




' (10 K cm ) 



FIG. 4 Amplitudes At and of the nonequilibrium tem- 
perature and transverse-momentum (viscous) fluctuations in 
liquid hexane at 25° C as a function of (Vr)^/k''. The sym- 
bols indicate experimental data. The solid lines represent the 
values predicted by Eqs. II2.34|I . Notice th e excellent agre e- 
ment with no adjustable parameters. From lLi et all (^9J). 



A, = 



{ki_ ■ VT) 

fc4 



(2.34c) 



For t = Q one recovers the equal-time density correla- 
tion function, Eq. (|2.33a|l . The amplitudes At and A^j 
are proportional to (VT)^/fc^, which has been verified 
by experiments, see Fig. 0] Notice that the amplitude of 
the temperature fluctuations is enhanced by a factor of 
a hundred compared to the scattering by an equilibrium 
fluid. In real space, this strongly singular wave num- 
ber dependence corresponds to a decay as const. — in 
three-dimensions. 



B. Quantum Systems 

The same mechanisms discussed above for classical sys- 
tems also lead to long-range correlations in space and 
time in quantum systems. The chief distinctions between 
quantum systems (most interesting at T = 0) and clas- 
sical systems are as follows. 

(1) There are more soft, or gapless, modes in quan- 
tum systems. These extra modes are due to a Gold- 
stone mechanism that is absent at T 7^ 0. For example, 
particle-hole excitations across a Fermi surface are soft 
at T = 0, but acquire a mass at T ^ 0. 

(2) As discussed in Sec. II. CI there is a coupling between 
the statics and the dynamics in quantum statistical me- 
chanics that is absent in the classical theory. In general, 
this means that long-ranged correlations in time correla- 
tion functions imply long-ranged spatial correlations in 



equal-time correlation functions, and thus in thermody- 
namic susceptibilities. 

In this subsection we explain these points in some de- 
tail, and lay the foundation for later sections where we 
discuss the importance of GSI for understanding many 
quantum phase transitions. Our focus will be on itiner- 
ant interacting electron systems. To motivate our discus- 
sions, and to make contact with the classical results dis- 
cussed in the preceding section, we first will give some of 
the results showing GSI in zero (or very low) temperature 
electron systems, both disordered and clean. We choose 
to discuss the disordered case first, because in this case 
the nonanalyticities that reflect the GSI are well known, 
although they usually have not been thought of as exam- 
ples of GSI. In noninteracting electron systems they are 
usually referred to as "weak-localization effects", and 
in interacting ones, as "interaction effects" or "Altshuler- 
Aronov effects" ; we will collectively refer to both classes 
as "GSI effects" .^^ The fact that there are analogous ef- 
fects in clean electron systems is much less well known. 
We then discuss the most important soft modes in elec- 
tron fluids, paying particular attention to those that exist 
only at T = 0. This is followed by an account of two dis- 
tinct approaches to GSI in quantum systems: The first 
approach is via explicit many-body calculations, which 
is analogous to a generalized kinetic-theory approach to 
the classical fluid. The second approach relies on the 
concept of clean and disordered Fermi-liquid fixed points, 
and uses very general RG arguments. This is analogous 
to the treatment of the classical Heisenberg ferromagnet 
discussed in Sec. III.A.l.cl 



1. Examples of generic scale invariance in itinerant electron 
systems 

a. Disordered systems in equilibrium In the entire metal- 
lic phase of disordered interacting electron systems, var- 
ious transport coefficients and thermodynamic quanti- 
ties show nonanalytic frequency and wave number de- 
pendencies that represent GSI. Since we are dealing 
with a quantum system, there also are corresponding 
nonanalytic temperature dependencies. For instance, 
for 2 < c? < 4 the electrical conductivity a, the 



We will refer to systems with and without quenched disorder as 
'disordered' and 'clean', respectively. 

The term "weak localization" is ill-defined, and difli'erently ap- 
plied by different authors. The most restrictive meaning refers 
to the weakly nonmetallic weak-disorder regime in d = 2, but we 
use it to denote the nonanalytic dependence of observables on 
the frequency or the wave number in noninteracting disordered 
electron systems in any dimension. 

Occasionally both classes have been collectively c alled "weak- 
localization effects" (e.g.. in lKirkpatrick et aILl2nn2fl . Since most 
people associate "weak localization" with quenched disorder at 
least to some degree, this use of the term tends to obscure the 
fact that precisely analogous effects are found in clean systems. 
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FIG. 5 Static low-T conductivity of ni ne Si:B samples in a 
magnetic field plotted versus Vt. After iDai efai] J1992D . 



specific heat coefficient 7 = Cy/T, and the static 
spin susceptibihty Xs) as functions of the wave vec- 
tor fe, the frequency f2, and tlic temper ature T, are 
given by (for reviews, see. lAltshuler and Aronovl Il984t 
iLee and Ramakrishnanl Il985|i 

a{n^O,T = 0)/<T{Q,0) = 1 + (2.35a) 

(7{n = 0,T ^0)/a{0,0) = l + S'^T^'^-'^^l'^, (2.35b) 

7(T^0)/7(0) = 1 + c]T('^-2)/2^ (2.35c) 

X.(fe-->0,r = 0)/xs(0,0) = 1-cf Ifel'^-^. (2.35d) 

Here the 0^'"*^'^° are positive coefficients that depend on 
the disorder, the interaction strength,^* and the dimen- 
sionahty. In d = 2 and d = 4 the fractional powers in 
these equations are replaced by integer powers times log- 
arithms. In the time domain, Eq. (|2.35al) implies that 
current correlations decay as 1/t'^/^, while Eq. (|2.35d|) 
implies that spatial spin correlations decay as l/r^^"*^^'. 

All of the above effects are examples of indirect GSI, 
as will become clear from the derivations given below. 
There are other examples of long-ranged correlations in 
addition to the ones given above; the latter were chosen 
because of their experimental relevance. For instance, an 
experiment showing the y/T dependence of the conduc- 
tivity expressed by Eq. (|2.35bp is shown in Fig. O It is 



of interest to note, however, that not all obvious corre- 
lation functions have nonanalyticities and exhibit GSI. 
For example, explicit calculations show that the density 
susceptibility, or compressibility, has no |fc|''~^ nonana- 
lyticity, in contrast to the s pin susceptibi l ity. E g. Ij2.35d|l . 
This has been discussed bv lBelitz et al\ l)2002(l . 

b. Clean systems in equilibrium For clean electronic sys- 
tems, two interesting nonanalyticities that reflect GSI 
are, for 1 < d < 3 (for a recent overview, see, 
IChubukov and MaslovLl2003a|) . 

7(T^0)/7(0) = l + 6;jr'^-i, (2.36a) 
Xs(fe ^ 0, T = 0)/x.(0, 0) = 1 + bf \kf-\ (2.36b) 

with the b'^'^' positive, interaction dependent coefficients. 
In d = 1 and d = 3 the integer powers in these equations 
become multiplied by logarithms. 

Let us explain the interesting sign difference between 
the nonanalytic corrections to the spin susceptibility in 
the clean and dirty cases, respectively, which will become 
important in Sec. lIVI Quenched disorder is known to in- 
crease the effective spin-triplet interaction streng th be- 
tween the electrons l|Altshuler and Aronovl Il984|) . This 
effect is strongest at zero wave number. Xs(0,0) is there- 
fore enhanced, but the enhancement effect decreases 
with increasing wave number, which leads to the negative 
correction in Eq. (|2.35dp . In the clean case, on the other 
hand, the generic soft modes represent fluctuations that 
weaken the tendency towards ferromagnetism. This ef- 
fect decreases Xs(0, 0) compared to the Pauli value (more 
precisely, the Pauli value with mean-field, or Hartree- 
Fock, corrections), and the wave number dependent cor- 
rection is positive. 

2. Soft modes in itinerant electron systems 

Let us now consider the soft modes that are respon- 
sible for the results listed above. Apart from the soft 
modes due to the standard conservation laws, which are 
the same as in the classical case, the most interesting 
soft modes in itinerant electron systems are the Gold- 
stone modes of a broken symmetry that is characteristic 
of quantum systems. To understand these new modes, we 
consider the formal action given by Eq. Ijl.9b|) in more 
detail. Specifying the Hamiltonian, we consider a model 
action 

S^So + Spot + 5int, (2.37a) 



Some of these nonanalyticities, e.g., tlie one in Eq. I2.35al . 
are presen t even in noni nterac t ing electron systems (see, e.g ., 
lKramer~a nd MacKinnon! Il993t iLee and Ramakrishnanl llQSSTl . 
which are an alogous to the classical Lorentz model (see, e.g., 
iHaugeL Il97j^ . We will be mostly interested in interacting elec- 
tron systems, which are analogous to classical fluids. 



Note that this represents an effect of disorder that enhances the 
tendency towards ferromagnetism. It comes in addition to the 
more basic, and mor e obvious, oppos ite effect due to the dilution 
of the ferromagnet iGrinstehl Il985^ . The important point for 
our argument is that the former effect is wave number dependent, 
while the latter is not. 
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with 



So = dx^i^aix) [-dr + VV2mc + i'aix). 

(2.37b) 

S'pot = / dx u{x)t/ja{x)ipa{x). (2.37c) 



Here x = {x, r) and J dx ^ J dx J^^ dr, nic is the elec- 
tron mass, /i is the chemical potential, and u{x) is a 
one-body potential, u can represent an electron-lattice 
interaction, but for the physical applications we are inter- 
ested in this would not lead to qualitative deviations from 
the behavior of a simple 'jellium' model of a parabolic 
band, and we will therefore not pursue this possibil- 
ity. We will, however, want to treat the interaction of 
electrons with quenched (i.e., stat i c) ran dom impurities 
l)Abrikosov et alL lEdwardsl [ 1958|) . In this case 

u{x) is a random function. It is governed by a distri- 
bution that we take to be Gaussian with zero mean, and 
a second moment 



{u{x)u{y)}dis = Uo 5{x -y)= . (2.38) 



Here {. . .}dis denotes the disorder average, Ti-oI is the elas- 
tic relaxation time, and N-p is the density of states per 
spin at the Fermi surface. The factors in Eq. H2.38|l are 
chosen such that disordered average Greens function has 
a single-particle lifetime equal to 1/2t^cI- We will con- 
sider both clean systems, where C/o = 0, and disordered 
ones, where C/q > 0. 

S'int in Eq. (|2.37a() represents the Coulomb potential, 
but it is often advantageous to integrate out certain de- 
grees of freedom to arrive at an effective short-range in- 
teraction. For our present purposes it is not necessary to 
specify the precise form of S'int; it will suffice to postu- 
late that the ground state of the interacting system is a 
Fermi liquid for [/q = 0, or a disordered Fermi liquid for 
Uq ^ 0.'^*' Ho wever, for later referenc e we write down a 
popular model ijAbrikosov et aZ.l . ll963(l that describes the 
self interaction of the electron number density, n{x)^ and 
the electron spin density, ns{x)^ via point-like, instan- 
taneous interaction amplitudes Fg and Ft, respectively, 



f dx n{x) n{x) + j dx n s{x) ■ n s{x) . 



In terms of fermionic fields, n and rig are given by 

n{x) = ^V^^(x) V^(x), (2.39b) 

a 

ns{x) = '^ll)a{x)crc,^a' 1pa'{x). (2.39c) 

Here cr = {a^ ,a^) represents the Pauli matrices. 
An underlying repulsive Coulomb interaction leads to 
Fg,Ft > in Eq. fT^ . 

It is useful to perform a Fourier representation from 
imaginary time to fermionic Matsubara frequencies w„ — 
2TTT{n+l/2), 

iJnA=^) = Vt dr e'^^-^Mx), (2.40a) 

"'0 

^n^x) ^ Vt dr e'^-^^^{x). (2.40b) 
Jq 

The noninteracting part of the action can then be writ- 
ten. 

So = So + Spot = / dxy^^'ipn,a{x) \iujn + V^/2mc 

<T,n 

+fi + u{x)]ij„^„{x). (2.41) 

In the disordered case, it is further convenient to inte- 
grate out the quenched disord e r by m eans of the replica 
trick llEdwards and AndersonL Il975t for a pedagogical 
discussion of this technique, see. iGrinsteinL 11985). Ac- 
cordingly, one introduces N identical replicas of the sys- 
tem, labelled by an index a, and integrates out u(x). 
Spot then gets replaced by 



Sa 



AttNf 



1 ^ r 

J E / dxY,^nAx)^nA^) 



a, 13=1 



x^i.A^X.'i^)^ (2-42) 



(2.39a) 



and physical quantities are obtained by letting — > at 
the end of calculations. 

We now no te two crucial features 

l|Wegner and Schafeil Il980|) . First, for w„ = 0, the 
action is invariant under transformations of the 
fermionic fields that leave V'nV'n invariant. That is, 
is invariant under a continuous rotation in frequency 
space. '^^ Second, this symmetry is spontaneously broken 
symmetry whenever N-p 7^ 0. To see this, consider the 
order parameter 

Q ^ lim {■iijn,a{x)lpn.cr{x)) 

- lim {^n,^{x)il;„.aix)). (2.43) 



A disordered analog of La ndau's Fermi-liquid theory 
(e.g. , iBav m and Pethic kl ^9^1) h as been developed by 
ICastcUani and_Pi^X]aatro| ll9 8^ ; IC^ ggtellani and Di Castrol 
Iigsal ^ -Castellani et aM il987alll983.ll987bD . 



Because of the anticommuting nature of the fermion fields, the 
Lie group in question is symplectic; for a model with 2M Mat- 
subara frequencies it is Sp (2M). 
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Q is the difference between retarded and advanced Green 
functions, and is proportional to iVp. The causal Green 
function has a cut along the real axis for all frequencies 
where Q oc A^f 7^ 0, which breaks the symmetry between 
frequencies with positive and negative imaginary parts, 
or between retarded and advanced degrees of freedom. 
Since Si-at cannot explicitly break this symmetry,^^ this 
will result in soft modes according to Goldstone's theo- 
rem. Technically, these soft modes are most conveniently 
discussed in terms of fluctuations of 4 x 4 matrices that 
are isomorph i c to bilinear products of fermionic fields 
l|Efetov et adll980|) . 



/ 



V 



-i/'iT'02T 
-i/'itV'2T 



-■'/'iT'02i 

V:uV:2i 

-V'iT'02i 



-■01TV'2i 

-V:u'02i 

--01TV'2i 



'01tV'2T \ 

V;U'02T 
-'0UV'2T 

V'itV'2t y 

(2.44) 

Here all fields are understood to be at the position x, 
and 1 = (ni,ai), etc., comprises both frequency labels 
n and, for the disordered case, replica labels a. Since 
a nonzero frequency explicitly breaks the symmetry, one 
expects, for nin2 < 0, 



lim ((3„i„2(fe) 



- UJ„ 



(2.45) 



For nin2 > 0, on the other hand, we expect the corre- 
lation function to approach a finite constant in the limit 
of small wave numbers and frequencies. That is, Qnm is 
soft if the frequencies ojn and tOm have opposite signs, and 
massive if they have the same sign. For later reference, 
let us introduce a notation that distinguishes between the 
soft and massive components of Q. We write 



qnm{x) if n> 0,771 <0 
= { qlmix) if 77< 0, 771 > 
Pnm{x) if 77777 > 



(2.46) 



Explicit calculations confir m these expe ct ations 
llBehtz an d Kirkp atrickL Il997l: lEfetov et all Il98nt 
IWegner an d Schafed. Il98ri) . For technical reasons, it 
is convenient to expand the 4x4 matrix given by Eq. 
(|2.44(l in a spin-quaternion basis, 



Qi2{x) = ^ (r^ (g) Si) lQ^^{x), 



(2.47) 



For small wave numbers and low frequencies, and 
for disordered no ninteracting electrons, one finds 
\VM.o^ et r,iJ.IT98fl33 

G 

('r9l2(^)i934(P)) = -^5{k+P)5rs % ^ig (524 2?„i„2 (fc): 

(2.48a) 



with 



T^n\n2 (^) 



1 



(disordered) 



(2.48b) 

Here H = nNp/A, and G — S/ttcto oc I/ti-oi, with (Jq the 
conductivity in Boltzmann approximation, see Eq. (|2.5t)|) 
below. 1/GH = D is the electron diffusion coefficient. In 
the low-disorder limit, D — VpT^ci/d, with 7;f the Fermi 
velocity. Note that this structure is consistent with Eq. 
(|2.45l) . which was based on very the general arguments 
and Goldstone's theorem. For clean noninteracting sys- 
tems, Eq. (|2.48a() is still valid, but G = 2/NfVf, and V 
is given by'^'* 



T^n\n2 (^) 



1 



|fe| + GH\uJni ~ iOn 



(clean) . (2.48c) 



H is the same in either case. As expected, for disordered 
electrons V is diffusive, while for clean electrons it is 
ballistic. 

For interacting systems, clea n or disordered, the re- 
sults are structurally th e same l)Belitz and Kirkoatrickl 
Il997t iFinkelsteinL 1198,'^ . The correlation functions re- 
main diffusive and ballistic, respectively, but they are 
no longer the same in the spin-singlet [i = 0) and 
spin-triplet (i = 1,2,3) channels. Rather, instead of 
H in Eqs. Il2.48bl I2.48c|l . the combinations H + 
and H + Kt, respectively, appear in the two channels, 
with Ks^ = 7rAr|rs,t/2. We wiU not detailed corre- 
lation functions for the interacting case; we just men- 
tion the physical significance of these various coupling 
constants. H determines th e specific heat coefficient 
7 = Cv/T via 7 — 8ttH/3 (jCastellani and Di Castrol 
Il986j) .^^ a nd can be interpreted as a quasiparticle density 
of states l|Castellani et all Il987btl . H + = ndn/dfi 
and H + Kt = ttXs are proportional to the thermody- 
namic density susceptibil ity dn/dn and the spin suscep- 
tibility Xs, respectively (jCastellani et al 1 Il984bl Il98fit 



with To = So the 2x2 unit matrix, and Tj = —Sj = —icrj, 
{j = 1,2,3). In this basis, 7 = and i — 1,2,3 describe 
the spin-singlet and spin-triplet degrees of freedom, re- 
spectively, r = 0, 3 corresponds to the particle-hole chan- 
nel (i.e., products of V'V'); while r = 1,2 describes the 
particle-particle channel (i.e., products or ip^). 



This follows generally from time translational invariance, and 
can be checked explicitly for the interaction defined in Eq. 



The normalizations of these propagators used in the literature 
vary, depending on how the Q-fields have been scaled. The nor- 
malization used in Eqs. I2.48i corresponds to a dimensionless 
Q{x), in accord with the usual convention in the nonlinear a 
model (see Eqs. I2.62i below), while the Q{x) as defined by Eq. 
I2.44i is dimensionally a density of states. 

Equation I2.48ci is a schematic representation, which has the 
correct scaling properties, of a more complicated function. In 
contrast, Eq. I2.48bl represents the exact propagator in the 
limit of small frequencies and wave numbers. 

This reference established the relation between -yy and H for 
diso rdered systems only. A later proof by means of Ward identi- 
ties iCastellani et alin^^ can be generalized to apply to clean 
systems as well. 
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iFinkelsteirJ . 119844 - Notice that dn/d^i, Xs, and 3^/2^"^ 
all are equal to Np for noninteracting electrons, but differ 
for interacting systems. 

The above discussion assumed that the underlying con- 
tinuous symmetry (see the discussion above Eq. (|2.43|) ) 
is not broken explicitly. Such an explicit symme- 
try breaking occurs, e.g., in the presence of a mag- 
netic field, magnetic impurities, or spin-orbit scattering 
llCastellani et 'aIl. ll984atlEfetov et al. l llQSnHFinkelsteinL 
ll984albf L In the presence of any of these symmetry 
breakers, some of the channels classified by the spin and 
particle- hole indices i and r of the matrix become 
massive. For instance, magnetic impurities lead to a mass 
in both the particle-particle channel (r — 1,2) and the 
spin-triplet channel {i = 1,2,3). A summary of thes e 
effects has been given bv lBelitz and Kirkpatricl3 l)l994|) . 
The nonanalyticities in various observables quoted in Sec. 
III.B. II above are cut off accordingly, depending on which 
soft channel they rely on. 



3. Generic scale invariance via explicit calculations 



One way to obtain the results quoted in Sec. lII.B.Tl is by 
explicit many-body, or Feynman diagram, calculations to 
lowest order in a small par ameter l Abrikosov ^7^ . 11 9631 
iFetter and Waleckal Il97l[) . To illustrate this approach, 
we compute the electrical conductivity a in a noninteract- 
ing, disordered, electronic system, to lowest nontrivial or- 
der in the disorder about the Boltzmann value. T he cal- 
culat ion starts with the Kubo expression for (j{n) llKubrl 



I 



7r(q = 0,i^ln) 



ne 



(2.49a) 



with TT the longitudinal current correlation function, 

TT{q, in^) f dT e'""^ (Tr q jig, r) q -/(q, 0) ' 

(2.49b) 

Here Tr is the imaginary time ordering operator, and j 
is the current operator, 



-q/2 



(2.49c) 



with and a electron creation and annihilation opera- 
tors, respectively. Wick's theorem can be used to evalu- 
ate time ordered correlation functions such as the one in 
Eq. p.49b|l . 

For our present purposes, the small parameter for a 
perturbative treatment is the impurity density n^, or, al- 
ternatively, l/kpi, with kp the Fermi wave number, and £ 
the mean-free path between electron-impurity collisions. 
The averaging over the positions of the impurities can be 
performed using stan dard techniques ijAbrikosov et all . 
119631: lEdwardIll958(l . The building blocks of the theory 
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FIG. 6 Diagrammatic representation of the Green function 
G'"' and the impurity factor uq. 



FIG. 7 The electronic self energy and the Green function in 
Born approximation. 



are, the bare-electron Green function. 



G^"\q,iLO^) 



dr e-"- (r.aq(r)at(0))g 



(2.50) 

and the impurity factor Uq, Eq. (|2.38l) . The subscript in 
Eq. H2.5()|l indicates that both the average and the time 
dependence of a are taken with the free-electron part 
of the Hamiltonian only. Diagrammatically, we denote 
G^°^ by a directed straight line, and uq by two broken 
lines (one for each factor of the impurity potential) and 
a cross (for the factor ofn.i ~ l/kp£); see Fig.|Sl For our 
free-electron model one has 

G(")(q, tLJn) = {iuJn - qV2m, + (2.51) 

The exact disorder-averaged Green function can be writ- 
ten in terms of a self energy S as 

G{q, iun) = {iujn - q^/2mc + /i + S(q, iujn)) ^ ■ 

(2.52a) 

The Born approximation for E and for G is shown dia- 
grammatically in Fig. [7| Analytically, it is given by 



E(q, iw„) 



2t,c 



sgn (w„). 



(2.52b) 



In this approximation, single-particle excitations thus de- 
cay exponentially with a mean- free time r^ci- One com- 
monly defines a (longitudinal) current vertex function, F, 
by writing tt as 

7r{q,inn) = [ TY,^oip,q)Gip-q/2,iuJn -inn/2) 
X G{p + q/2, ibJn + inn/2) F(p, g; iw„, il7„). (2.53a) 



Here Jp = J dp/{27T)'^ in d-dimensions, and 
ro(p,g) = P ■ q/\q\ 



(2.53b) 



is the bare current vertex. It is often convenient to ex- 
press F in terms of another vertex function. A, via the 
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FIG. 10 Crossed-ladder contribution to A. 






+ 



FIG. 8 Diagrammatic representation of the current correla- 
tion function, and of the relation between the functions F and 
A. 



FIG. 9 The vertex function A in Boltzmann approximation. 

equation 

r(p,q;iw„,il7„) = ro(p,q') 

+ / ro(fc,q)G(fc-q'/2,it^„-if]„) 
Jk 

xG{k + q/2, iujn + if^,i/2) A{k,p, q; iojn, i^n)- 

(2.54) 

Diagrammatically, these equations are illustrated in Fig. 

El 

The Boltzmann approximation Ab for A (which is 
equivalent to solving the Boltzmann equation for a) is 
shown diagrammatically in Fig. |51 Analytically, we find, 
for small il„ and q, 



ABik,p,q;iLO„,iQn) 



= e (-cj„(cj„ + fi„)) 



Dq^ + \nn\ 



(2.55) 



Using this in Eq. H2.54|l . and the result in Eqs. H2.53a|) 
and H2.49a|l . the Boltzmann approximation for a is^^ 



(To = n e^Trcl/TOo 



(2.56) 



There are, of course, an infinite number of corrections to 
CTo of 0{l/kF£) and higher. Of particular interest to us 
are the contributions that couple to the electronic soft 
modes, Eqs. H2.48|l . As in the classical LTT calculation, 



Sec. III.A.3.bl it is the coupling of these modes to the cur- 
rent fluctuations that lead to the slow decay of current 
correlations, and to a singular dependence of the conduc- 
tivity on the frequency. The relevant diagrams for this 
soft-mode contribution to a are th e so-called maximall y 
crossed diagrams shown in Fig. ^ ^ Gorkov et cfII ■ ll979^ . 
For these contributions to A, which we denote by AmCj 
one finds 

^Mcik, P,q;ii^n,'i'^n) = 9 (-w„(w„ + f^„)) 

^i?(fc+p)2 + |a„| ■ ^2-57) 

Notice that this i s just the Boltzmann result , Eq. H2.55|l . 
with q^ k+p llVollhardt and Wolfd. ll980V For smaU 
real frequencies this contri bution leads to the following 
result for the conductivity l)Clorkov et al\ . \l97^ . 



1 



-in + Dq^ 



(2.58) 



In d = 2, this 'correction' to cto leads to a (negative) loga- 
rithmic divergence. This reflects the fact that, at least for 
n oninteracting electrons, there is no metallic state in d ~ 
2 ifAbrahams et a?J . ll979j) . This phenomenon (the "weak 
localization" proper, see footnote I26|l has generated a 
huge body of literature (for revie ws, sec, Ber gmann, 19^ 
Kramer and MacKinnoni 1199.4 Eee and Ramakrishnam 
1985|) . To the extent that it is caused by the LTT not 
being integrable, this is analogous to the breakdown of 
classical hydrodynamics in d < 2 (see Sec. III.A.3.b|l . al- 
though the physics behind it is quite different. 

In d > 2, Eq. H2.58|l predicts that a{n,) is a nonanalytic 
function of frequency, the small-frequency behavior for 
2 < d < 4 being 

a{n 0) /cr(0) = 1 + const, x rj(''-2)/2^ (2.59a) 

With const. > 0. This particular correction is one of sev- 
eral that contribute to Eq. (|2.35al) : in the nomenclature 
explained at the beginning of Sec. III.BI it is the weak- 
localization contribution. Anot her one is the interact ion 
contribution that was found bv lAltshuler et al\ lll98Cll) .^^ 
In the time domain this corresponds to a behavior of the 



Due to our point-like impurity potential, F in Boltzmann ap- 
proximation is equal to Fo, and the conductivity is given in 
terms of the single-particle relaxation time r^^i . For a more gen- 
eral scattering potential, the vertex corrections in Eq. 112.541 are 
nonzero, and the Boltzmann conductivity is given by Eq. I2.56i 
with Trpi replac e d by the transport relaxation time rtr- See, 
lAbrikosov et all <1963D Sec. 39.2. 



Historically, the interacti on effect preda tes weak localization. 
It was first discussed bv |Schmid| ] 1974f) for the ele ctron in- 
elastic lifetime, and by lAltshuler and Aronol 119791) for the 
density of states, and a closely related effect was found by 
iBrinkman and Engelsberd ll968f) . 
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current-current correlation function, Eq. (j2.49b|) . at long 
(real) times t, 



7r(q = 0,i ^ oo) l/f^^^ 



(2.59b) 



Since the current is not conserved in our model, this is an 
example of indirect GSI. Notice that the exponent is the 
same as for the corresponding LTT in a classical fluid, Eq. 
(|2.31a|l . but the sign is different. Physically, this means 
that, in the quantum case, the leading LTTs decrease the 
diffusion coefhcient, while in the classical fluid case they 
increase it. This is because the scattering in the electron- 
impurity model is due to static impurities, while in the 
classical fluid it is due to the other fluid particles. "^^ 



4. Generic scale invariance via renormalization-group arguments 



The results given in Secs. lII.B.ll and lil.B.3l can also be 
obtained independent of perturbation theory, by general 
RG and scaling arguments which establish these LTTs 
as the ex act leading behavior for l ong t imes or small fre- 
quencies (jBelitz and Kirkpatrickl 1199^ . These authors 
showed that in the disordered case, the GSI effects in 
d > 2 can be understood as corrections to scaling at 
a stable zero-temperature fixed point that describes a 
disordered Fermi liquid. This also makes more explicit 
the connection between weak-localization and Altshuler- 
Aronov effects on one hand, and the Goldstone modes of 
the 0{N) nonlinear cr model fSec. III. A. II) . or the hydro- 
dynamic theory of GSI fSec. lII.A.4|) . on the other. This 
is the first instance where we see a close, if not entirely 
obvious, analogy between classical and quantum effects. 



of the effective theory then becomes a two-step process. 
First, one needs to identify the massless modes of the 
system, cf. Sec. III. B. 21 Second, the microscopic theory, 
Eqs. (|2.37I2.39|I . must be transformed into an effective 
one that keeps only the soft modes, while all other de- 
grees of freedom are integrated out in some reasonable 
approximation. The result of this procedure is a general- 
ized nonlinear a model. Within this theory the partition 
function can be written as an integral over a matrix field 
Q that essentially contains the soft sectors of the field Q, 
Eq. 



D[Q] 



(2.60) 



The effective action A is given by 



A 



-1 
2G 



dx triVQ{x)f + 2H 



dx tr{nQ{x))+Aint- 
(2.61a) 

Here Aint ^ O(TrQ^) is the interaction part of the ac- 
tion, i.e., S'int from Eq. (|2.39a(l expressed in terms of Q- 
matrices. Schematically, leaving out the detailed struc- 
ture of the frequency, replica, and spin-quaternion labels, 
it is of the form""^ 



Aint^TT / dx Q{x)Q{x). 



(2.61b) 



For our present purposes, F can stand for either Fg or Ft. 
17 is a diagonal matrix with fermionic Matsubara frequen- 
cies as diagonal elements, and the coupling constants G 
and H were defined after Eq. (|2.48b|) . Q is subject to the 
constraints 



ix)^l , = Q , trQ{x) = 1. 



a. Disordered electron systems The long-wavelength and 
low-frequency excitations in a noninteracting disor- 
dered electron systems ca n be d escribed by a nonlin- 
ear a model (M cKanc and Stonel [1981; Weener, 1979; 
IWegner and Schafeit 1198(1 " This has been generalized 
to interacting svst ems bv IFinkelsteinI ll'983l): a deriva- 
tion in t he spirit oflWeg ner and Schafei lll98n|^ has been 
given by iBehtz and Kirkpatric^lfl997(l and lBelitz et all 
|)1998|) . As for the simpler a model derived and discussed 
in Sec. UFA. II the basic idea is to construct an action 
solely in terms of the massless modes. The derivation 



This makes the electron-impurity model more closely analogous 
to the classical Lorentz model, which shows a weaker 
LTT (with a negative sign) than the classical fluid I Haugel ll974D . 
The mode-mode coupling effects are stronger in the quantum 
system than in the classical one, which is why the electron model 
has a LTT with the same strength as the classical fluid, but with 
the sign of the classical Lorentz model. 

We note in passing that it is possibl e to cast classical hy drody- 
namics in the form of a field theory iMartin et all 119731) . This 
technique has not been applied to the LTTs so far; doing so 
would make the analogy even closer. 



(2.62a) 



A standard way to enforce these constraints, analogous 
to Eqs. H2.6|l . is to write Q as a block matrix in frequency 
space, 



Q = 



(1 



tU/2 



-(1 



t,)l/2 



(2.62b) 



where the matrix q has elements qnm with n ^ 0, m < 0. 

Insight into the GSI effects in a disordered metal 
can be gained by an RG analysis that focuses on a 
stable fixed point of the ab ove generalized cr model 
p5clitz and Kirkpat ricklll997j) . This fixed point provides 
an RG description of the disordered Fermi liquid ground 
state (see footnote I3U|) . in analogy to Shan kar's R G de- 
scription of a clean Fermi liquid ( ShankaJ Il994l) . The 
procedure in analogy to Sec. III.A.l.cl We again define 
the scale dimension of a length L to be [L] = — 1. The 



The complete expression IFinkelsteinI I1983D will not be needed 
here. 
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stable disordered Fermi-liquid fixed point is characterized 
by the choice 

[q{x)]^-{d~2)/2 (2.63) 

for the scale dimension of the field q, which corresponds 
to diffusive correlations of the q. This choice is consistent 
with what one expects for the soft modes in a disordered 
metal, see Eq. (|2.48b|) . In addition, the scale dimension of 
the frequency or temperature, i.e., the dynamical scaling 
exponent z = [f2] = [T], is needed. In order for the fixed 
point to be consistent with diffusion, the frequency must 
scale as the square of the wave number, so we choose 

z = 2. (2.64) 

Now we expand the action in powers of q. In a symbolic 
notation that leaves out everything not needed for power 
counting purposes, we have 

A ^ J dx {\Iqf + H j dx riq^+TT J dx q^ 

+ 0{V^q\nq\Tq''). (2.65) 

Power counting shows that all of the coupling constants 
in Eq. H2.65|l have vanishing scale dimensions with re- 
spect to our disordered Fermi- liquid fixed point, 

[G] = [H] ^ [F] = 0. (2.66) 

Now consider the leading corrections to the fixed-point 
action, as indicated by Eq. H2.55|l . Power counting shows 
that all of these terms are irrelevant with respect to the 
disordered Fermi-liquid fixed point as long as d > 2. Fur- 
thermore, all of the terms that were neglected in deriving 
the generalized nonlinear a model can be shown to be 
even more irrelevant than the ones considered here. The 
conclusion from these arguments is that the term given 
explicitly in Eq. (|2.65|) constitute a stable fixed-point ac- 
tion, and that the leading irrelevant operators (which we 
denote collectively by u) have scale dimensions 

M = -(d-2). (2.67) 

These results can be used to derive the GSI effects 
from scaling arguments. We first consider the dynamical 
conductivity a(i2). Its bare value is proportional to G, 
and according to Eq. H2.66|l its scale dimension is zero. 
We therefore have the homogeneity law (cf. Sec. II.B|I 

a{n, u) ^ a{n b',u b-^'^-^'>). (2.68a) 

By putting b = 1/0^/^, and using z = 2 as well as the 
fact that (t(1, a:;) is an analytic function of x, we find that 
the conductivity has a singularity at zero frequency, or 
LTT, of the form 

a{n) = a{n ^0) + const, x Q^'^-^^/'^. (2.68b) 

That is, we recover Eqs. (|2.35a|) and H2.59a|) . The present 
analysis proves (with the same caveats as in Sec. lII.A.l.cp 



that the fi^''"^)/^ is the exact leading nonanalytic behav- 
ior. 

The specific heat coefficient 7 = Gy /T is proportional 
to the quasiparticle density of states H, whose scale di- 
mension vanishes according to Eq. (|2.66|) . We thus have 

7(T,w) = 7(T6",u6-('^-2))^ (2.69) 
which leads to the low-temperature behavior given by Eq. 

Finally, we consider the wave vector dependent spin 
susceptibility Xsik). Xs is a two-particle correlation and 
thus can be expressed in terms of a Q-Q correlation func- 
tion. The leading correction to the finite Fermi-liquid 
value is obtained by replacing both of the Q by q. We 
thus have structurally Xs ^ T J dx {q^ q), with scale di- 
mension [xs] = 0. The corresponding homogeneity law is 

Xs{k,u)^Xs{kb,ub-^''-^y), (2.70) 

which leads to a nonanalytic dependence on the wave 
number as given by Eq. H2.35c|) . 

We conclude with an important caveat. As all scaling 
arguments, those presented above are exact, but very 
formal. There is no guarantee that the prefactors of the 
various nonanalyticities, i.e., the coefficients c^''*^''"''', are 
nonzero. Indeed, explicit calculations show that both 
and are nonzero on ly in the presence of electron - 
electron interactions fe.g.. 'Alts huler and Aronovlll984() . 
Further, the same arguments would lead to the conclu- 
sion that the density susceptibility has the same form as 
Eq. (|2.70|l for the spin susceptibility, while perturbation 
theory finds that the corresponding prefactor vanishes 
even in the presence of interactions, see the discussion 
at the end of Sec. III.B.l.al The situation is thus as fol- 
lows. If explicit calculation show that a nonanalyticity 
has a nonzero prefactor, then the scaling arguments can 
be used to show that the perturbative result is indeed the 
leading nonanalytic behavior, which one can never con- 
clude from perturbation theory by itself. However, the 
scaling arguments by themselves can never establish the 
existence of a particular nonanalyticity. 

b. Clean electron systems Structurally identical argu- 
ments can be made for clean electronic systems 
l)Belitz and Kirkpatrickl Il997(l . There are two motivat- 
ing factors. First, while the weak-localization effects are 
clearly caused by disorder, disorder plays a much less 
crucial role in the explicit calculations for the Altshuler- 
Aronov effects that lead to the same type of nonan- 
alyticities, and the scaling arguments presented above 
do not at all depend on disorder in any obvious way. 
Second, it is well known that in d = 1, a perturba- 
tive expansion in powers of the interaction breaks down 
due to logarithmic divergencies, and that this divergence 
is related to the breakdown of Fermi-liquid theory in 
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this d imension ijDzvaloshinski and Larkinl Il97lt ISchul j . 
Il995j) . A natural question is, what happens to these 
singularities for d > 1? Explicit perturbation theory 
shows that both clean and disordered interacting elec- 
tronic systems have related nonanalyticities that reflect 
the GSI in both systems, and that these singularities 



down of Fermi liquid theory in d = 1 ( Belitz et at. 


199(1 


Chitov and Minisll2nnil: 


Chubukov and Maslovll2003albl 


Galitski and Das Sarma 


120031). 



The formal scaling arguments for these effects proceeds 
as in Sec. III.B.4.al above. The only structural differ- 
ence is that the two-point vertex in the disordered case, 



p(2) 
disordered 



k^+GHQ, is replaced by 



(2) 
clean 



\k\+GHn 



in the clean case, and in the terms of higher order in q, 
the operator in the disordered case is effectively re- 
placed by a V operator. In terms of scale dimensions, 
the net result of these substitutions, using the notation 
of the previous subsection, is 



[q{x)] = -(d-l)/2, 
z = 1, 

[u] = 



(2.71a) 
(2.71b) 
(2.71c) 



with u denoting the least irrelevant variables. Note that 
technically the difference between the disordered and 
clean case is that (d — 2) in the disordered case is re- 
placed by (d — 1) in the clean case. 

Using the Eqs. H2.71|l and arguments identical to those 
used in the disordered case, one obtains the results given 
by Eqs. H2.36(l . As for the disordered case, this reasoning 
shows that the Eqs. H2.36|l are exact. 



5. Another example of generic scale invariance: Disordered 
systems in a nonequilibrium steady state 

The electron-impurity model in a nonequi- 
librium steady state has been studied 
l)Yoshimura and Kirkpatrickl Il996j) by means of 
many-body diagrammatic techniques in conjunction 
with Z ubarev's nonequilibrium statistical operator 
method l|Zubarevlll97^ . This work showed that effects 
analogous to those discussed in Sec. III. A. 41 for classical 
nonequilibrium fluids also exist in disordered electron 
systems at T = 0, and that quantum GSI effects due 
to nonequilibrium conditions are stronger, i.e., of longer 
range, than in classical systems. 

The particular model studied was that defined by Eqs. 
(|2.37|) and H2.38|l . without the electron-electron interac- 
tion term S'int, but in the presence of a chemical-potential 
gradient V/i. Consider the nonequilibrium part of the 
electronic structure factor 



(2.72) 



where (. . .) denotes a nonequilibrium thermal average. 
As one might expect, the calculation leads to a result 



that is analogous to the one in the classical fluid, Eqs. 
viz., 



C2{k) 



6d7r(L>fe2)2 



25(V/i)2 - 12(fc • Vn)^ . (2.73) 



As in the classical fluid, this corresponds to a decay of 
correlations in real space given by const. — |a;| in three- 
dimensions. 

In comparison, for a classical Lorentz gas teaTigeLll974l 
see also foonotes EHl and one flnds instead 



C2(fc) ^ (VA*)Vfc' 



(2.74) 



That is, in both quantum and classical cases C2 exhibits 
GSI, but the effect is stronger in the quantum case due 
to the existence of more soft modes. 



6. Quantum Griffiths phenomena: power laws from rare regions 
in disordered systems 

In this section we discuss another mechanism for long- 
range correlations in time in disordered quantum sys- 
tems, the so-called quantum Griffiths phenomena. In 
contrast to the examples above, the relevant soft modes 
are local in space. 

Griffiths phenomena are nonperturbative effects of rare 
strong-disorder fluctuations. They were flrst discovered 
in the context of classical pha se tra nsitions in quenched 
disordered systems (GriffithS. Il969l) . Griffiths phenom- 
ena can be understood as follows: In general, impurities 
will decrease the critical temperature Tc from its clean 
value T°. In the temperature region Tc < T < T° the 
system does not display global order, but in an infinite 
system one will find arbitrarily large regions that are de- 
void of impurities, and hence show local order, with a 
small but nonzero probability that usually decreases ex- 
ponentially with the size of the region. These static dis- 
order fluctuations are known as "rare regions" , and the 
order parameter fluctuations induced by them belong to 
a class of excitations known as "local moments" ; some- 
times they are also referred to as "instantons" . Since they 
are weakly coupled to one another, and flipping them 
requires changing the order parameter in a whole re- 
gion, the rare regions have very slow dynamics. iGriffith^ 
( 1969) was the first to show that they lead to a nonana- 
lytic free energy everywhere in the region Tc < T < Tc, 
which is known as the Griffiths phase, or, more appro- 
priately, the Grifflths region. In generic classical sys- 
tems, the contribution of the Griffiths singularities to 
thermodynamic (equilibrium) observables is very weak 
since the singularity in the free energy is only an essen- 
tial one. In contrast, the consequences for the dynamics 
are much more severe, with th e rare regions dominat- 
ing the behavior for long ti mes ijBravL 11985 iDhar et al\ . 
119881: iRanderia et alVmS^ . In the classical McCoy- Wu 
model, where the disorder along a certain direction is 
correlated, the Grif fiths singularities have drastic effects 
even for the statics (|McCov and WuLll968lll969|) . 
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At quantum phase transitions, the quenched disor- 
der is perfectly correlated in one of the relevant dimen- 
sions, viz., the imaginary time dimension. Therefore, the 
quantum Griffiths effects arc gcncri cally as strong as in 
the classical McCoy- W u model ( Rieeer and Youn3 . ll996t 
iThill and Hi'^1199,^ . One of the most obvious realiza- 
tions of quantum Griffiths phenomena can be found in 
random quantum Ising models. 



H 



(2.75) 



with Pauli matrices and cr^ 
over nearest neighbors on a 
couplings Jij and the transverse fields hi 

has 



and the summation runs 
lattice. Both the Ising 
are indepen- 
been studied 



model 
r 



dent random variables. This 

extensively in d = 1 jFisheij, '1995; Igloi et al, 199£ 
Igloi and Rieeeil Il998t lYoungl [l997: Young and Riegei 
19961) and d = 2 llMotrunich e< a/.l l200(]l: jpich et al 



199S : iRieger and Youn3 ~ 1996j) ! and the results are 
expected to qualitativ ely hold in d — 3 as well 
llMotrunich et "all. 1200(1) . Let us consider a system which 
is globally in the disordered phase, with a rare region of 
volume Vr that is locally in the ordered phase. The prob- 
ability density for such a region to occur is exponentially 
small, p oc exp(— oVr), with a a disorder dependent coef- 
ficient. However, since the volume de pendence of the lo- 
cal energy gap A is also exponential ijYoung and Riegeil 
IT99Q) . A ~ exp(-6VR), the probability density P(A) for 
findi ng a gap of si ze A varies as A to a power. One thus 
has l|Youngl ll997|) 



P{A)^piVR)\dVR/dA\o,A-/'>- 



(2.76) 



Thus, we obtain a power-law density of low-energy ex- 
citations. Notice that the exponent a/b (which in the 
literature is often called an inverse dynamical scaling ex- 
ponent, and denoted by 1/z) is a continuous function of 
the disorder strength. 

Many results follow from this. For instance, a region 
with a local energy gap A has a local spin susceptibility 
that decays exponentially in the limit of long imaginary 
times, XIoc{t — ^ oo) oc exp(— Ar). Averaging by means 
of P, Eq. P7^ . yields 



Xloc 



(r oo) (X T-°/^ 



(2.77a) 



The temperature dependence of the static average sus- 
ceptibility is then 



dT 



(2.77b) 



If a < 6, the local zero-temperature susceptibility di- 
verges, even though the system is globally still in the 
disordered phase. The typical correlation functi on falls 
off mu ch faster than the average one; for d = 1 . lYoun j 
lll997^ has found a stretched-exponential behavior 



xloci'^ ^ oo) (X exp 



(2.77c) 



with c another coefBcient. The large difference between 
the average and the typical values is characteristic of very 
broad probability distributions of the corresponding ob- 
servables. 



7. Quantum long-time tails from detailed balance 

Finally, we discuss a very simple source of long time 
tails, or temporal GSI, which is operative at zero tem- 
perature only. To be specific, consider a model with spin 
diffusion. The imaginary, or dissipative, part of the dy- 
namical spin susceptibility as a function of wave number 
and real frequency is 



x"{k,n)^xik) 



Dk^n 



+ {Dk^f 



(2.78) 



with D the spin diffusion coefficient and x{k) the static 
spin susceptibility, x" {k,^) is an analytic function of 
frequency. Upon Fourier transformation, this yields an 
exponential decay of correlations in real time, x{k, t) oc 
exp(— ZJfc^i). Now consider the dynamic structure fac- 
tor, S{k,Vl), and the symmet rized fluctuatio n function, 
</7(fe, il), which are defined by ijForsteil Il975j) 



S{k,n) 

^ik,n) 



— dx 



dt e-'''^-^+*^"(M(a;,i) Af(O,0)), 

(2.79a) 

OO 

dt e~*'"-"'+'"*(Af(a;, t) Af (0, 0) 

-OO 

+M{O,0)M{x,t)), (2.79b) 



2 ' '^'^ 



where M is the magnetization operator. (We neglect 
the vector nature of M for simplicitv-) The flu ctuation- 
dissipation th eorem ijC alien and Weltonl Il95ll see also 
iForsteil Il97l relates these two correlation functions to 



X via 

Sik,n) ^ ^_l^pn X"ik,n), (2.80a) 

ipik,n) = coth{pn/2) x"ik,fl). (2.80b) 

For any nonzero temperature, these correlations are again 
analytic functions of the frequency, and their Fourier 
transforms therefore decay exponentially in time. How- 
ever, in the zero-temperature limit, (3 —^ oo (and for 
fc 7^ 0), both of these correlation functions are not ana- 
lytic at 17 = 0; namely, they are proportional to 17 0(17). 
This nonanalyt icity lea ds to a power-law decay in real 
time (see, e.g.. lLighthillll95^ . 



>5'(A;, t oo) oc ip{k, t oo) oc 1/t^ 



(2.81) 



Physically, the nonanalyticities in the Eqs. H2.80|l at T = 
reflect the absence of excitations in the ground state. 
This is contained in the detailed balance relation 



S[k,-Vl) 



-fjn 



s{k,n), 



(2.82) 
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which imphes that S{k, f2) at T = must vanish for 
< 0. Note that these long-time tails are present at 
all wave numbers, not just in the long- wavelength limit. 
More generally, the above considerations imply that S 
and ip will display LTTs for any zero-temperature system. 
For low but nonzero temperature there will be a pre- 
asymptotic LTT followed by an exponential asymptotic 
decay. 

For later reference we also consider the spin suscepti- 
bility at imaginary Matsubara frequencies iri„ = i2TrTn, 

(2.83) 

which is a nonanalytic function of Qn- Accordingly, 
the imaginary time correlation function x(^: ''") — 
T ^„ e-*""^x(fc, i^n) at T = has a LTT for large imag- 
inary time T, 

x{k,T ^ oo) (x1/t^. (2.84) 

More generally, any imaginary time correlation function 
has a long imaginary time tail if the corresponding causal 
function is nonanalytic at zero frequency. 

III. INFLUENCE OF GENERIC SCALE INVARIANCE ON 
CLASSICAL CRITICAL BEHAVIOR 

We now combine the concept of GSI with critical phe- 
nomena. Given that critical phenomena are driven by 
soft modes, one expects that soft modes connected to 
GSI, provided they couple to the order-parameter fluctu- 
ations, will influence the critical behavior. This is indeed 
the case in both classical and quantum systems. In this 
section we will discuss two classical examples. The first 
one deals with the nematic-smectic-A transition in liq- 
uid crystals, where the generic soft m odes are the Gold- 
stone modes o f the nematic phase. Jch^ "dl. Il978t 
iHaloerin et a/.l . li974|l This turns out to be closely related 
to the normal metal-superconductor transition, where 
the generic soft modes are virtual photons. The second 
one considers phase separation in a binary liquid sub- 
ject to shear, and the generic soft modes are due to th e 
nonequilibrium situation (Onuki and Kawasaki. ri979|) . 
Another classical example, which we will not dicuss here, 
involves the ferromagnetic transition in compressibl e 
magnets l)Aharonvlll976l: iBergman and Haloerinl . [l976fl . 
Quantum phase transitions will be discussed in Sec. IIVI 
below. 



A. The nematic-smectlc-A transition In liquid crystals 

Nematic liquid crystals consist of rod-shaped 
molecules. In the nematic phase, there is direc- 
tional order, i.e., the rod axes tend to be parallel to one 
another, while the centers of gravity of the molecules 



do n ot show long-range order ijPeGennes and Prostl 
Il993|) . In the smectic-A phase, the molecules are 
additionally arranged in layers. Within each layer, the 
rod axes are on average aligned perpendicular to the 
layer, but the molecules still form a two-dimensional 
liquid. In either phase, there are Goldstone modes 
associated with the broken rotational invariance that is 
characteristic of the directional order. It turns out that 
these Goldstone modes couple to the order parameter 
for the smectic-A order, and have an important influence 
on the cr itical behavior. Let u s now discuss this effect, 
following IHalperin et gj] lll974l) . Ichen et all l)l978D . and 
IDeGennes and Prostl l)l993|) . 

1. Action, and analogy with superconductors 

The layers of the smectic-A phase are described as a 
density modulation in, say, the z-direction, 

p{x) = po + pi{x) cos {q z + ip{x)) . (3.1) 

In the smectic-A phase, the fluctuating amplitude pi{x) 
has a nonzero mean value, q is the wave number associ- 
ated with the smectic layer spacing, and is a phase. It 
is convenient to combine pi with to form a complex or- 
der parameter (see fPcGcnnes and Prostl (^19931 and ref- 
erences therein) 

i;{x) = pi{x) e"^^''\ (3.2) 

One piece of the action for the smectic-A order then takes 
the form 

^aM = Jdx {romx)f + c\VtP{x)\^ + uomx)\^} . 

(3.3a) 

This is the LGW functional for a complex scalar field, or, 
equivalently, the 0{2) version of Eq. (|2.3|l .'^^ rg and uq 
are bare parameters that will be renormalized later. This 
action is incomplete. One must add to it, (1) a term de- 
scribing the fluctuations of the nematic order parameter, 
and (2) a term describing the coupling between the two. 
The nematic rods are described by a unit vector (called 
the "director") n, and the director fluctuations are given 
by A{x) — n{x) — z. To linear order in the fluctuations, 
n? = \ enforces A = (ylj,, Ay, 0). The Gaussian contri- 
bution of A to t he action takes the form of t hree gradient- 
squared terms ijDeGennes and Prostlll99^ . viz., 

5n[A] ^ j dx {a'i [V • Af + if2 [i • (V X A{x))f 

+ K:,[dA{x)/dzfY (3-3b) 



The gradient terms parallel and perpendicular to the layers have 
different coefficients, but the action can be made isotropic by an 
appropriate real-space scale change. 
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where Ki, K2, and are related to elastic constants. 
The coupling between ip and A can be shown to be ad- 
equately represented by replacing the VV' in Eq. (j3.3a|l 
by 

V%Ij{x) ^ [V - iqA{x)]il:{x). (3.3c) 

We see that the coupling between the order parameter 
and the director fluctuations takes the same form as the 
one between the order parameter and the vector potential 
in Eq. I|2.15a|) . It is important to realize, however, that 
the action given by Eqs. H3.3|l . in contrast to Eqs. H2.15|l . 
is not invariant under local gauge transformations due to 
the additional gradient terms in Eq. Ij3.3bp . This has pro- 
found consequences for the ordered phase; for instance, 
it leads, via a Landau-Peierls instabili ty, to the a bsenc e 
of true long-range order in a smectic l|Lubenskvl 
see also footnote 1421 below *) . However, it has been shown 
that a simplified action obtained by neglecting Ki and 
puttin g K-2 = = I/ Sttij, has the same critical prop- 
erties l)Lubenskvl Il983|) . Although the full action has 
properties that are in general very different from those 
of the gauge theory given in Eq. H2.15a|l . the former thus 
reduces to the latter for the purposes of determining the 
critical behavior. Keeping in mind that A lies in the 
xy-plane, an appropriate effective action therefore is 

^na[^,A] = Jdx {roWx)\^ + c \[W ^ iqA{x)]^{xf 
+ uomx)\^ + ^[VxA{x)f}. (3.4) 

As we have discussed in Sec. III. A. 21 this action formally 
is gauge invariant. We also recognize it as the ac- 
tion of a spin-singlet superconductor, with tp the super- 
conducting order parameter, A the vector potential, q 
the Cooper p air charge, and jj , the magnetic permeabil- 
ity (see, e.g.. iDeGennesl [l989l) . The nematic-smectic-A 
transition can thus be mapped onto a superconductor- 
normal metal transition, even though the generic soft 
modes in the superconducting case, viz., the photons, 
have a very different origin than the director fluctuations: 
The latter are the Goldstone modes of a broken symme- 
try, while the former are the result of gauge invariance, 
as was discussed in Sec . HTl This remarkable analogy 
l)Halperin and Lubenskvl Il974|) is a good example of the 
versatility of effective field theories. This becomes even 



more apparent upon the realization that the action H3.4(l 
is closely related to the scalar electrodyn amics problem 
studied bv lColeman and Weinberg! (|^23), which demon- 
strated a mechanism for the spontaneous generation of 
mass in particle physics. 

We make on final remark to put this section into con- 
text. In the disordered phase we have generic long-ranged 
correlations, or GSI, from the gauge field fluctations, cf. 
Eq. H2.17a(l . In the following subsections we will see how 
these generic long-ranged correlations modify the critical 
behavior compared to the one that would result from Sa, 
Eq. (|3.3a|) . alone. In particular, see Eq. H3.7|) below, and 
the discussion following it. 

2. Gaussian approximation 

The simplest way to deal with the action ^na is to 
treat it in Gaussian approximation, that is, to neglect 
the nonlinear coupling between A and ip. The ip'^ term 
needs to be kept for stability reasons, of course. In this 
approximation the model reduces to an XY-model, which 
makes the liquid crystal analogous to a superfluid rather 
than to a superconductor, and it pre dicts a continuou s 
transition in the XY- universality class l)DeGennesl . ll972j) . 
As we know from Sec. III. A. 21 and will see again below, 
this Gaussian approximation misses some crucial aspects 
of the full model. 

3. Renormalized mean-field theory 

The action H3.4|l is characterized by two length scales. ^"^ 
In a superconducting language, they are, the coherence 
length 

e = y/^l (3.5a) 
and the London penetration depth 

A = fc-i = Vl/87r/^cg2(|^(a;)|2). (3.5b) 

Here r is the renormalized distance from the critical 
point, see below. The ratio k = A/^ = l/k\^ is the 
Landau-Ginzburg parameter. Notice that k is indepen- 
dent of r, since (j-i/'l) esc y^jrj. k < 1 characterizes type-I 
superconductors, where the order parameter coherence 
length is larger than the penetration depth. In a mag- 
netic field, these materials form a Meissner phase that 
completely expels the magnetic flux.'*'* In the extreme 



It is important to realize, however, that the only physical gauge 
is still the one where = 0. "Gauge transformations" that 
take one to, e.g., the Coulomb gauge used in Sec. 111. A. 21 lead 
to a transformed order parameter which no longer represents 
a physical observable. For instance, the transformed order pa- 
rameter exhibits true long-range order in the smectic phase, 
whereas, as mentioned in the text, the orig inal one does not 
ennes and Pros 1 1l993t iLubenskvl Il983l) . Another way to 
say this is that the phase of the order parameter in a smectic is 
observable, while in a superconductor, it is not. 



The full action for the liquid crystal problem contains many 
more length scales, and hence more cases than the type-I and 
type-II of the superconductor need to be distinguished. These 
have not been fully classified. 

To avoid misunderstandings, we stress that the conclusions 
drawn below are valid in the absence of an external magnetic 
field, as the vector potential A can describe spontaneous electro- 
magnetic fluctuations. 
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type-I limit, k <C 1, order-parameter fluctuations are neg- 
ligible, ij{x) Ri {i'ix)) = -0, and the part of the action 
that depends on the vector potential takes the form 



1 

87r/i 



dx 



kl A^{x) + (V X A{x)y 



We see that in the normal conducting phase, the trans- 
verse photons are soft, as they should be according to 
Sec. III.A.2.al (these are the generic soft modes), while in 
the superconducting phase they acquire a mass propor- 
tional to k\ (X ^p, cf. Sec. III.A.2.bl Choosing a gauge 
fixing term, Eq. H2.16|l . with rj — 1^^ one finds for the 
A-propagator 



{Aa{x) Ap{y)) = Safj S{x - y) 



V2 



(3.6) 



Since the vector potential no longer couples to any fluc- 
tuating fields, it can be integrated out to obtain an action 
e ntirely i n terms of the superconducting order parameter 
fHalperin et adll974 - 



Sc«m = In J D[A] e 



= V (rol^P + uolVI") - 2 Tr In (V^ + kl) . (3.7) 

Here Sna is given by Eq. I|3.4|l . V is the system volume, 
and we have neglected a constant contribution to the ac- 
tion. The above procedure is exact to the extent that 
fluctuations of ip can be neglected; otherwise, it is an 
approximation. It produces an action entirely in terms 
of ij}, but, since one has integrated out a soft mode, the 
price one pays is that this action is nonlocal.'*^ The ef- 
fect of the nonlocality is most conveniently studied by 
expanding the equation of state in powers of the order 
parameter if), and the free energy density / — T S^fi/V 
can be obtained by integrating the result or der by order. 
The details have heen given bv IChen~ei' ;iDri978). The 
result for the leading terms in three-dimensions is 

//r = r|V^|2-t;3|^|3+y|^|4_/,^ (d = 3). (3.8a) 

Here we have added a field h conjugate to the order 
parameter, r and u are given by and uq with ad- 
ditive corrections proportional to fiq^ and (fiq^)^, re- 
spectively. These changes do not affect the behavior 



Notice that the renormalized mean-field theory neglects all fluc- 
tuations of the field in contrast to the Gaussian theory dis- 
cussed in Sec. III.A.2.bl In particular, no particular gauge needs 
to be chosen in order to eliminate the Higgs field. One should 
also keep in mind that the full action for the liquid crystal, which 
is not gauge invariant due to the additional gradient terms in Eq. 
I3.3bi . has a different soft-mode spectrum in the ordered phase 
than the superconductor. 
*® More precisely, the mode that has been integrated out has a 
small mass (for small \ip\) that is given by the order parameter. 



of the theory. However, > is a positive coupling 
constant proportional to that drives the transi- 

tion first order. Notice that the new term generated by 
the generic soft modes is nonanalytic in the order pa- 
rameter; such a term can never appear in a local Lan- 
dau expansion. This provides a rather extreme exam- 
ple of generic soft modes infiuencing critical behavior: 
The critical point is destroyed, and a first-order transi- 
tion takes place instead. This phenomenon is known as 
a fluctuation-induced first-order transition in condensed 
matter physics, and as the Coleman- Weinberg mecha- 
nism in high-energy physics. Interestingly, it is not the 
soft mode per se that produces the first-order transition; 
it is the fact that a nonzero order parameter gives the 
soft mode a mass. This is an examp le of a more general 
principle that has been explored bv lBelitz et al\ l|2002'l . 
For later reference we also quote the free energy density 
in four dimensions l)Chen et a/.l Il978(l , although this case 
is not of physical relevance for liquid crystals, 

//T = rlV-P + v^\i\j\^ In + wlV'l^ -h^ = 4), 

(3.8b) 

with U4 > 0. 

Note that the free energy functionals given by Eqs. 
(|3.8|) are nonanalytic in ■(/;, and that this nonanalytic- 
ity has nothing to do with the phase transition. Indeed, 
they denote the proper free energy for the order param- 
eter caused by a conjugate field far from the transition. 
Physically, the nonanalyticity reflects GSI in the disor- 
dered phase, and it is directly related to a nonanalytic 
wave number dependence of the order parameter suscep- 
tibility x^. This becomes obvious if one remembers that 
in any mean- field theory, r — \ ~ Fxi/. , with F the appro- 
priate interaction coupling constant, and that V' scales 
like a wave number far from the phase transition, as can 
be seen from Eqs. Ij3.5bl I3.7|l . In a scaling sense, Eq. 
(|3.8a|l thus corresponds to a wave number dependence of 
Xv- given by 



X^(fc)/X^(0) = l-H6'^'^|fc| 



(3.9a) 



with 1)^^ > a positive coefficient. Notice that this is 
analogous to Eqs. (|2.36bl l2.35d|) . although Xi/j is not ob- 
servable for the superconductor, and related to an ob- 
servable susceptibility in a complicated way for the liq- 
uid crystal. In Sec. IIVI we will discuss examples where 
the magnetic susceptibility, which is directly measurable, 
plays an analogous role. In real space, this behavior cor- 
responds to 



Xs(r > 0, \x\ ^ cx)) cx l/|a;|^ 



(3.9b) 



This manifestation of GSI is ultimately responsible for 
the failure of the Gaussian theory to correctly describe 
the phase transition. 



Notice that the fluctuations in question are not the order- 
parameter fluctuations, but rather the generic soft modes. 
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The conchision that the nematic~smectic-A transition 
and the BCS superconducting transition are of first or- 
der is inevitable if order-parameter fluctuations are neg- 
ligible. Quantitatively, the effect turns out to be un- 
observably small in superconductors, as the first-order 
nature of the transition is predicted to manifest itself 
only within about one /iK from the transition tempera- 
ture. In liquid crystals, on the other hand , the n umbers 
are much more favorable l)Halperin et all Il974|l . After 
a long time of confusion, careful experiments have in- 
deed confirmed the first-or der nature of the transition in 
a var i ety of liquid crystals llAnisimov et 1199(1 iLdidi^. 
l200lHYethirai et oilEiol . 



4. Effect of order-parameter fluctuations 

Interesting and technically difficult questions arise 
when order-parameter fluctuations are taken into ac- 
count. First-order transitions are often thought of as 
unaffected by order-parameter fluctuations, almost so by 
definition, since they preempt a fluctuation-driven crit- 
ical point. However, if the ffi'st-order transition takes 
place within the critical region of an unrealized second- 
order transition, then the critical fluctuations associated 
with the latter can destabilize the mechanism that drives 
the transition ffist order. The resulting transition can 
then be continuous, but described by a fixed point differ- 
ent from the one that destabilizes the fluctuation-induced 
first-order transition. An explicit example for such a 
"fl uctuation-induced second -order transition" was given 
by iFucito and Parisl l)l98l|) . How this works technically 
becomes clear once one realizes that, within a RG treat- 
ment, a fluctuation-induced first-order transition comes 
about by the quartic coupling constant u flowing to neg- 
ative values. Critical fluctuations can weaken this ten- 
dency, and lead to a positive fixed point value of u after 
all. We will discuss an explicit example of this happening 
in Sec. IIV.AI below. 

For the superconducting/nematic-smectic-A transi- 
tion problem, a RG analy s is to one-loop order has been 
performed bv lChen et al\ ^jB)- To first order in e = 
4 — d, these authors found no critical fixed point and con- 
cluded that the transition is always of first order, both 
for superconductors and for liquid crystals, and for both 
type-I and type-H mat e rials. However, later work by 
iDasgupta and HaloerinI 1)1981*). prompted by mounting 
experimental evidence that the transition is of second 
order in type-II materials, concluded that, as one enters 
into the type-II region, the first-order transition becomes 
weaker and weaker until the transition reverts to sec- 
ond order in the so-called inverted XY-universality class. 
This has been confirmed by a sizea ble body of numeri- 
cal and analytical evidence (see, e.g.. lHerbut et g/Tl200ll 
and references therein). Why the perturbative RG does 
not show this fixed point is not quite clear. 



B. Critical beliavior in classical fluids 

Our second classical example deals with the criti- 
cal behavior of a classical fluid, both in equilibrium 
ijHohenberg and Halperin. 1977: Kawasak i, 1976), and in 
a nonequilibrium si tuation created by applyin g a con- 
stant rate of shear (jOnuki and Kawas akf. '1979^. Let us 
consider a fluid conflned between two parallel plates sub- 
ject to a constant rate s of shear, see Fig. ^2 The average 
local velocity is given by 

v{x) ^ sye^, (3.10) 

with Ex a unit vector in x-direction. This externally in- 
duced velocity has to be added to the Navier-Stokes equa- 
tions H2.21|l . Close to the liquid-gas critical point, these 
equations can be replaced by the following set of simpler 
equations l)Kawasa id 119701 lOnuki and Kawasakil[T979,) . 

(3.11a) 
(3.11b) 

Here the order parameter %p represents the specific en- 
tropy, u is the deviation of the local velocity from its 
average and A and rj are the bare thermal conductiv- 
ity and the bare shear viscosity as in Eqs. H2.21|l .^^ v± 
denotes the transverse part of a vector v. S" is a classical 
action, or free energy functional, given by 

Slip]^ J dx [r^A^ + c(V^)2 -l-uV'^] , (3.11c) 

and 9 and are Langevin forces that obey 

{e{x,t)0{x',t')) = -2XV^6{x-x')6{t-t'), 

(3. lid) 

{Co.ix,t)Cp{x',t')) =. -25^pr^V^5{x-x')5{t~t'). 

(3. He) 




FIG. 11 Velocity field in a fluid subject to a steady, plane 
Couette flow. 



Here we discuss a pure fluid for simplicity. For experimental 
purposes, binary fluids are preferred for technical reasons. This 
changes the interpretation of the various quantities in Eqs. 13. Hi . 
but the physics remains the same. 
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In the a bsence of shear, s = 0, Eqs. 13.11(1 describe 
Model H of lHohenberg and HalperinI ||1977.1 . It is a time- 
dependent Ginzburg-Landau theory for the conserved or- 
der parameter coupled to the conserved auxiliary field 
u. Since the shear modes are soft, this coupling influences 
the critical behavior. We will now discuss this influence, 
first in equilibrium, and then for a system with s ^ 0. 



1. Critical dynamics in equilibrium fluids 

Let us return to the discussion of long-time tails in 
Sec. III.BI There, we considered corrections to the bare 
kinematic viscosity v in the GSI region far from critical- 
ity, and mentioned that analogous results hold for other 
transport coefficients. Of particular interest for the crit- 
ical dynamics near the liquid-gas critical point is a con- 
tribution to the correction to the thermal conductivity, 
SX, from the coupling of the transverse-velocity fluctua- 
tions to the entropy fluctuations described by Eqs. ((3.11() 
for s = 0. Before performing the wave number integral, 
this part icular contribution, which we deno te by <5A_ls, is 
(plohcnb erg and HalperiiA Il977t iKawasakl ^76) 

^ P i 

(3.12) 

Here Dt = X/pCp is the thermal diffusivity in terms of 
Cp, the specific heat at constant pressure, x is the order 
parameter susceptibility for the phase transition, where 
we have anticipated that near the critical point we will 
need the momentum-dependent %, and p± = p ± k/2. 
Setting k = and carrying out the momentum integral 
leads to a t^'^^^ LTT. The correction to the thermal con- 
ductivity is obtained by integrating SX{t) over all times, 
as in the case of the kinematic viscosity.*^ 

By examining Eq. H3.12|l one easily identifies a mecha- 
nism by which the LTT effects can become even stronger. 
Consider a system with long-range static correlations, for 
instance due to Goldstone modes, or due to the vicinity 
of a continuous phase transition. In either case, some 
susceptibilities, e.g. the x i^i EQ- l|3.12fl . become long- 
ranged, amplifying the LTT effect. Before we discuss the 
realization of this scenario in the vicinity of a phase tran- 
sition we make one last point concerning the LTTs. So 
far we have stressed the leading tails that decay as t^'^^^, 
but there are numerous subleading LTTs as well. Most of 



These results imply that for d < 2, conventional hydrodynamics 
does not exist. Indeed, it is now known that for these dimensions 
the hydrodynamic equations are nonlocal in s pace and time. For 
a discussion of this topic, see, Forster et al. ( 1973) ^.nd references 
therein. In smectic liquid crystal phases this effect is stronger, 
since some susceptibilities behave, for certain directions in wave 
vector space, as 1/k*, amplifying the LTT effect. This causes 
a breakdown of local hydrodynamics for all dimensions d < 5 
jMazenko et adll983l) . 



them are uninteresting, but one becomes important near 
the critical point, via the mechanism discussed in the last 
paragraph. According to Eq. 1(3.12(1 . a central quantity 
for determining the critical contribution to A is the shear 
viscosity rj (which enters v, see the definition of v after 
Eq. 1(2.24(1 .). It turns out that the contribution to ry that 
is dominant near the critical point is a subleading LTT 
away from criticality. It involves a coupling of two heat 
or entropy modes and is given by ijKawasakl Il976j) 

^,(M) = ^i:^:(p)x(>'-p)(j5)-j^)' 

x(kx-p)'e-^-[p'+(''-P)')l* , (3.13) 

with A a constant. Away from the critical point, this 
LTT decays as t~^'^/^~^^\ so in fact it is a next-to-ncxt 
leading LTT. Nevertheless, it is the dominant mode- 
coupling contribution to rj near the critical point because 
of the two factors of x in the numerator of Eq. ((3.13() . 

Near continuous phase transitions, fluctuations grow 
and ultimately diverge at the critical point. For the 
liquid-gas critical point the order parameter is the dif- 
ference between the density and the critical density, and 
the divergent fluctuations are the density fluctuations as 
described by the density susceptibility. The susceptibility 
X in Eq. ((3.13() is proportional to this divergent suscepti- 
bility. In the Ornstein-Zernike approximation it is given 
by Eqs. 1(1.2(1 . Away from the critical point, r ^ 0, x 
decays exponentially in real space, see Eq. ((1.2bl) . 

Carrying out the time integral, the leading singular 
contribution to the static, wave number-dependent t her- 
mal conductivity is l|Hohenberg and HalperinI Il977() 

Using Eq. ((1.2a(l in this equation we see that the homoge- 
neous thermal conductivity is infinite at the critical point 
for all d < 4, diverging as |r|~('*~'^)/^. This is result of 
the amplification of the LTT by the critical fluctuations. 

By the same mechanism, we see that Eq. ((1.2a() leads 
to a logarithmically singular contribution to Si], if we take 
into account that Eq. ((3.141) implies that at the critical 
point, DT(k) - Ikl''-^^ 

The above conclusions result from one-loop calcu- 
lations that use the Ornstein-Zernike susceptibility. 
To go beyond these approximations it is necessary 
to, (1) use the correct scaling form for the suscep- 
tibility, and, (2) use either a self- c onsist ent mode- 
coupling theorv jKadanoff and Swiftl. Il968l) or a RG 
appro ach ijForster et all Il977t iHohenberg and HalperinI 
11977ft to improve on the one-loop approximation. The 
result is that the thermal conductivity diverges as 

Aoc|r|-0-" , (3.15a) 

and that the thermal diffusivity vanishes as 

Dt oc |r|0-67 . (3.15b) 
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FIG. 12 The thermal difFusivity Dt = X/pcp of carbon diox- 
ide in the critical region as a function of density at various 
temperatures (Tc = 304.12 K). The symbols indicate exper- 
imental data for Dt measured directly, and for X/pCp de- 
duced from thermal-conductivity data. The solid curves rep- 
resent values calculated from themode-coupling theory. After 
iLuettmer-Strathmann et al\ lllfl 



This is in good agreement with experimental results, as 
shown in Fig. 1121 which compares experimental and theo- 
retical results for the thermal diffusivity of carbon dioxide 
in the critical region. 

All of the above results are obtained by considering the 
effects on the dynamics (i.e, on transport coefficients) 
of the static critical behavior, as expressed by the sus- 
ceptibility X, which has been determined independently. 
There is no feedback of the critical dynamics on the stat- 
ics. This is an example of the missing coupling between 
the statics and the dynamics in classical equilibrium sys- 
tems that we alluded to in the Introduction. As we will 
see, this changes in a nonequilibrium situation. 



2. Critical dynamics in a fluid under shear 

Let us now consider Eqs. H3.11|l in the presence of 
shear, s ^ 0. In this situation the entire fluid can still be 
at its critical temperature and density. One therefore still 
expects a sharp phase transition (in an infinitely large 
system), albeit a nonequilibrium one. This is in contrast 
to driving the system out of equilibrium by means of, 
e.g., a temperature gradient. Inspecting Eq. H3.11a|l . we 
make the following observations. (1) Shear will have a 
tendency to make the order parameter susceptibility less 
soft. This is obvious since the operator yd/dx in Eq. 
(|3.11a|l . which scales like a mass by naive power count- 
ing, enters additively to the diffusive operator V^. In 
addition, the susceptibility will become anisotropic. (2) 
1/s sets a new time scale in the problem, which needs 
to be compared to the relaxation time r. One therefore 



expects equilibrium critical behavior in the weak-shear 
region sr < 1, and a crossover to different behavior in 
the strong-shear region sr > 1. Alternatively, one can 
define a characteristic wave number kg by T{ks) = 1/s. 
The strong-shear region is then given by k^^ > 1, with f 
the order parameter correlation length. 

It follows from point (1) above that the upper critical 
dimension cannot be greater than the one in the equi- 
librium case, which is equal to 4. At least for d > 4, 
the nonlinearities in the stochastic equations (|3.11l) must 
therefore be irrelevant. In the linearized equations, and 
u decouple, and the equation for the former becomes 



dtp 
'dt 



(3.16) 



The static order parameter susceptibility. 



1 f dn f , 

- T y 2^7 



(V'fe(i)V'-fc(O)), (3.17) 



can be obtained from Eq. H3.16|l Fourier transforming and 
using Eq. (|3.11d() . This yields the following differential 
equation for x-ip^ 



Xk^r + k^)--.sk, 



d 
dk,. 



(3.18) 



In the strong-shear region, t he solution of this equation i s 
adequately represented by l)Onuki and Kawasakil Il979j) 

X4,{k) = l/{r + const, x fcf/^|/C:,p/^ + fe^), (3.19a) 

with the constant of 0(1). This is the generalization 
of the Ornstein-Zcrnike susceptibility, Eq. H1.2a|l . to a 
sheared system. (Remember that x in Eq. (|3.14|) is 
proportional to Xtp-) expected, x^(fc) is strongly 
anisotropic, and less soft than for s — 0. The nonequilib- 
rium situation induces long-ranged static order parame- 
ter correlations even away from criticality that manifest 
themselves in the nonanalytic wave number dependence 
of Xii ■ In real sp ace in d = 3, one finds powe r-law corre- 
lations at r ^ ijOnuki and KawasakilllQT^ . 



X^(r > 0, \x\ ^ cxd) a 1/ 



.|7/5 



(3.19b) 



for |a;^| <^ k~^ ^ \x\, where x = x± + xex- For \x±\ » 
on the other hand, x^p decays exponentially. That is, 
exhibits GSI (and extreme anis otropy). At criticality, 

r = , Xip is of even longer range ijOnuki and Kawasakil 

Il979|) . If we use Eq. pri9a|l for r = in Eq. we 

find 



5A(fe = 0) oc / dpp'^-3^7«/5|p,|-2/5 
Jo 

cx fc,"«/^ / dp (3.20) 

In contrast to the equilibrium situation, where 5X{k — 0) 
diverged for all d < 4, we see that in the presence of shear 
the divergence occurs only for d < 12/5. 
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The above result suggests that the critical behavior 
of the thermal conductivity A is mean- field like, i.e., 
given by the time-dependent Ginzburg-Landau theory, 
for d > 12/5, and in particular in d = 3. The physical 
reason is that the long-range order parameter correla- 
tions stabilize the mean- field critical behavior. It is not 
obvious that this conclusion is correct, though, as it has 
been derived by using the approximation (|3.19a|l . which 
neglected the nonlinearities in the equations of motion. 
lOnuki and Kawasakil l)l979j) have performed a RG analy- 
sis which shows that the simple argument given above is 
indeed correct. In the presence of shear, the upper crit- 
ical dimensionality is d+ = 12/5, and for d > d+ there 
is a new simple critical fixed point where and 6X are 
given by Eqs. (|3.19a|l and (|3.2()|l . respectively. The RG 
treatment shows that the flow equations for Xip and the 
transport coefficients are coupled. This means that, con- 
trary to the equilibrium case, it is not possible to solve 
for the static critical behavior indep endently of the dy- 
namics. lOnuki and Kawasakil l)l979|) also considered the 
equation of state, which they found to be of mean-field 
form for d > 12/5. In particular, the critical exponent /3 
has its mean- field value f3 = 1/2. 

Another theoretical prediction regards the suppression 
of the critical temperature by the shear. Consider the 
inverse correlation length, as a function of r and s. 
Since s is an inverse time, it is expected to scale with the 
dynamical critical exponent z. Dynamical scaling then 
predicts the relation 

= r'^ri(l,s/r''^), (3.21) 

with b an arbitrary scale factor, and v and z the cor- 
relation length exponent and the dynamical critical ex- 
ponent, respectively, at the equilibrium transition. For 
s 7^ 0, the inverse correlation length thus vanishes not 
at r = 0, but rather at r = const, x s^/"^^. The shear 
dependence of the critical temperature, for small shear, 
is therefore 

Tc(s) = Tc(0) - const, x s^l''\ (3.22) 

A t the equilibrium transition, v » 0.63, and z w 3 in d = 
3 l|Hohenberg and Halr)erinLll977(l . Many of these results 
have been confirmed experi mentally in light-scattering 
exper iments on binary fiuids l|Bevsens an d Gbadamassi, 
For the Tc-suppression, the experimental results 
are shown in Fig. ^] together with the theoretical pre- 
diction. 

The classical nonequilibrium phase transition de- 
scribed above is very much analogous to the equilib- 
rium quantum phase transition in a ferromagnet with 
quenched disorder, which we will discuss in Sec. IIV.BI 
In that case, too, the order parameter susceptibility re- 
flects long-ranged correlations, although they are due to 
generic soft modes rather than a nonequilibrium situa- 
tion. This in turn leads to the stabilization of a simple 
Gaussian critical fixed point, much like in the nonequi- 
librium classical fiuid. An important difference is that in 
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FIG. 13 Critical temperature versus shear in the binary 
fluid cyclohexane-anihne. Symbols represent the experimen- 
tal data, and the solid line corresponds to the theo retical ex- 
ponent Ijvz K, 0.53. After ISevsens and Gbadamassi. (.1980>) . 



the quantum ferromagnetic case, the equation of state is 
not mean-field like, see Eq. (j4.24b|l below. The reason 
for this difference is not known. 



IV. INFLUENCE OF GENERIC SCALE INVARIANCE ON 
QUANTUM CRITICAL BEHAVIOR 

We now turn to quantum systems. As we have dis- 
cussed before, the only qualitative difference compared 
to the classical case is the coupling between the statics 
and the dynamics in quantum mechanics, which leads to, 
e.g., long-time tails affecting the critical behavior of ther- 
modynamic quantities even in equilibrium. Otherwise, 
the phenomena are qualitatively very much analogous. 
However, there is an important quantitative difference. 
In any itinerant electron system, there are generic soft 
modes, viz., the particle- hole excitations that can be un- 
derstood as Goldstone modes as explained in Sec. lII.B"2l 
which are soft only at zero temperature. Our first ex- 
ample for which these particle-hole excitations play the 
role of the generic soft modes is the case of the ferromag- 
netic transition in metals at T = 0. We will treat clean 
and disordered systems separately, since they turn out 
to behave quite differently, and are qualitatively remark- 
ably similar to the classical nematic-smectic-A transition 
discussed in Sec. IIII.AI and to the critical dynamics of a 
classical fiuid under shear. Sec. IIII.B.2I respectively. An- 
other example is the transition from a normal metal to 
a BCS superconductor T = 0, which we cover in Sec. 
IIV.CI The case of the quantum antiferromagnetic tran- 
sition, which is very different and at this point rather 
incompletely understood, is discussed in Sec. IIV.DI 



A. Quantum Ferromagnetic Transition in Clean Systems 

Let us consider a clean itinerant quantum Heisenberg 
ferromagnet as our first example of a quantum phase 
transition that is influenced by generic soft modes. In 



33 



Sec. lIV.Bl we will discuss the effects of quenched disorder 
and also give some additional technical details. 

Our exposition does not follow historical lines. The 
first detailed theory of the quantum ferrom agnet i c tran - 
sition was given in an influential paper bv yHert j l)l976|) . 
who derived an effective action from a microscopic model 
and concluded that the critical behavior is mean-field like 
in all dimensions d > 1. As we will show, this conclusion 
was the result of an approximation that ignores all but 
the most basic aspects of the generic soft modes. 



1. Soft-mode action for itinerant quantum ferromagnets 

The order parameter is the fluctuating magnetization 
field M{x) (with x = {x,t) the space-imaginary-time 
coordinate as in Sec. III.B.2|I whose expectation value is 
proportional to the magnetization m. The generic soft 
modes in this system are the diffusive particle-hole exci- 
tations. They were parameterized in Sec. lII.B?2l in terms 
of the soft sector q of the matrix field Q, see Eqs. (|2.46l 
12.471) . with the corresponding propagator given by Eq. 
(|2.48c|l . The coupling between the two is provided by 
the fact that the magnetization couples linearly to the 
spin density, which can be expressed in terms of the Q. 
The action will thus consist of a part that depends only 
on the magnetization, a part that depends only on the 
generic soft modes, and a coupling between the two, 

A[M,Q]=AM + Aq + AM,q (4.1) 

Am is a static, local, LGW functional for the magnetiza- 
tion fluctuations. It can be chosen static because, as we 
will see in Sec. IIV. A.2l below. the low-frequency dynami- 
cal part will be generated by the coupling to the ballistic 
modes. We thus write 

Am[M] = J dx [M{x) (ro - cV^) M{x) + uq M^{x)] . 

(4.2) 

Aq must yield the ballistic propagator of the soft modes, 
Eqs. (jngjl . The Gaussian part of the fcrmionic action 
will therefore have the form 

Af^^^fdxdy J2 E^9i2WrS34(^-y);'Z34(y)- 

1,2,3,4 r,i 

(4.3a) 

The particle-particle degrees of freedom are irrelevant for 
this problem, and we can therefore restrict the sum over 
the index r to r = 0, 3. The vertex function T^^^ is most 
easily written in momentum space, 

"'^ilsiik) = Si3S2iT^if\k) + 61-3,2-46,0 27:TGK, 
+ 61-3,2-i (1 - 6m) 2TrTGKf (4.3b) 

with 

rf2°'(fe) - |fe| + GTO_2 . (4.3c) 



The inverse of F^^^") yields the noninteracting propaga- 
tor, Eqs. (|2.48al IT^Sc|) . and the inverse of E^^^, its gen- 
eralization to interacting systems. A's is the spin-singlet 
interaction amplitude defined after Eq. H2.48c|l . and Kt 
is a 'residual' spin-triplet interaction.^" 

AM,q originates from a term Am-q that couples M 
and Q. Such a term must be present since in the presence 
of a magnetization the fcrmionic spin density will couple 
linearly to it. Using Eq. H2.44|l to express the spin density 
in terms of Q, we thus obtain 

r ^ 

Am-q = 2ciVT / dxY,J2^'^n{x) 

X (-l)''/^^tr[(Tr® s,)Q™^„+„(a;)] , 

r=0,3 m 

(4.4a) 

with a model-dependent coefficient ci. Defining a sym- 
metrized magnetization field by 

bi2{x) ^ (r^ ® s,) lb^^{x) , (4.4b) 
with components 

n 

(4.4c) 

allows to rewrite Eq. 14.4a|l in a more compact form, 

Am-q = ciVr j dx tr {b{x) Q{x)) . (4.4d) 

Using Eq. H2.46|l in Eq. H4.4a|l or H4.4d|) . and integrating 
out the massive P-fluctuations, obviously leads to a series 
of terms coupling Af and q, M and g^, etc. We thus 
obtain AM,q in form of a series 

AM,q = AM-q + AM-q2 + ■■■ (4.5a) 



One might argue that the spin-triplet degrees of freedom are 
included in M, so also including them in A.q constitutes double 
counting. To see that this is not true, imagine deriving the effec- 
tive action from a microscopic model, e.g., Eqs. 12,3712.391 . This 
introduces M by decoupling the spin-triplet interaction term in 
Eq, I2.39a t by means of a Hubbard-Stratonovich transformation 
I Hubbard, 19591: IStratonovichl 1195^ . which leaves spin-triplet 
degrees of freedom in the noninteracting part of the action So, 
(2) 

and hence in Aq , 

A weaker version of the same objection is that Kt should be 
zero, since it represents the interaction that has been decoupled. 
However, this argument ignores the fact that a spin-triplet inter- 
action will be generated from the spin-singlet one under renor- 
malization, even if there is none in the bare action. The only 
restriction on Kt is therefore that it must not be so large as 
to induce a ferromagnetic instability in the electronic "reference 
system" described by Aq in the absence of a coupling to M. 
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The first term in this series is obtained by just replacing 
Q hy q and , respectively, in Eq. (|4.4dp . 



In the context of the full action, A4 gives the Gaussian 
M propagator. 



12 •' r i=l 

(4.5b) 

The next term in this expansion has an overall structure 

Am-q^ oc C2VT J dx tr (^b{x) q{x) q\x)'^ , (4.5c) 

with C2 another positive con stant. The detailed struc- 
ture {Kirkoatrick and Bclit3, l2002j) can be obtained by 
integrating out the P-fluctuations in tree approxima- 
tion, in analogy to the deri vation of the nonlinear a 
model in the disordered ca se l)McKane and Ston3 . ll98lt 
IWegner and Schafeil Kml . Terms of higher order in q 
in this expansion will turn out to be irrelevant for deter- 
mining the behavior at the quantum phase transition. 

Let us pause here and compare our soft-mode action 
with the one for the classical liquid-crystal transition in 
Sec. nil. A. II The order-parameter part of the action, 
Eqs. (|3.3al) and (|4.2|l . respectively, is a (/(''-theory in either 
case, and the generic soft modes are described by a Gaus- 
sian action, Eqs. (|3.3bp and H4.3a|l . respectively. Finally, 
in either case there is a direct coupling between the order 
parameter and the generic soft modes. In the case of the 
liquid crystal, this coupling is between the square of the 
order parameter and the square of the soft-mode field, 
while in the case of the magnet the order parameter cou- 
ples linearly to all powers of the soft-mode field, but this 
does not have major physical consequences, as we will 
see. Apart from this, the main difference is that in the 
quantum case the fields depend on time or frequency in 
addition to position or wave vector, which leads to a more 
complicated detailed structure of the terms in the action, 
the coupling terms in particular. Given these structural 
similarities, it is natural to analyze the quantum ferro- 
magnetic transition in analogy to the liquid crystal one, 
Sees. IIII.A.21 and IIII.A.3I 



2. Gaussian approximation 

A strict Gaussian approximation neglects the term 
^M-gS, Eq. (|4.5c|l . One can then integra te out t he in- 
trinsic soft modes q to obtain an action ijHertzl. Il976|) 



Ak[M] ^^M{k)M-^{k)M{-k) + 0{M^), (4.6) 

A; 

where k = (fe, f2„) is a momentum-frequency four-vector, 
d — 4:Gci/tt, and A4^^ is the inverse of the paramagnon 
propagator 



(M;(fe) M^{p)) = <5„,-™ S,, - Mnik). (4.7b) 



M and the fermionic propagator T), Eq. (|2.48c|l . with the 
latter suitably generalized to allow for the interaction am- 
plitudes Kg and Kt, are the Gaussian propagators of the 
coupled field theory.^^ The |r2„|/|fc| structure of the dy- 
namical piece of Ai is characteristic of the itinerant elec- 
tron degrees of freedom that couple to the magnetic ones; 
it also manifests itself, e.g., in the fr equency dependence 
of th e Lindhard function (see, e.g., iPines and NoziereS. 
Il989j) . Technically, it results from the inverse of the ver- 
tex T^'^\ Eq. Ij4.3b|l . which gets multiplied by a frequency 
due to frequency restrictions inherent in the definition of 
the q. This illustrates how the coupling to the generic 
soft modes generates the dynamics of the order parame- 
ter, see the remark above Eq. (|4.2|) . 

An is Hertz's action, which predicts a continuous tran- 
sition with mean-field critical behavior for all dimensions 
d > 1. From our experience with the analogous classical 
transitions in Sec. IIII.Al above we suspect that the Gaus- 
sian approximation yields qualitatively incorrect results. 
Indeed, internal inconsistencies of H ertz's results have 
been d iscussed bv lSachdevI l)l994() and lDzero and Gorkovl 
(j2P0v?). More explicitly, the instability of Hertz's fixed 
point can be shown formally by means of arguments anal- 
ogous to those we will present in Sec. IIV.B.4.bl bclow for 
the case of disordered magnets. In what follows, we show 
that the transition can be either of first order, or of sec- 
ond order with non-mean-field critical behavior. 



3. Renormalized mean-field theory 

As a first step to improve upon the Gaussian approx- 
imation, it is natural to construct a renormalized mean- 
field theory in analogy to Sec. IIII.A.3I The relevant 
length scales are, the magnetic coherence length 



(4.8a) 



and a length 

A = l/GH^C2y/T{Ml{x)) = 1/Ac2 m, (4.8b) 

with m is the magnetization. A — G^J ttH^ /8 depends 
on microscopic length and energy scales only and en- 
sures that TO is dimensionally an inverse volume. In a 
magnetic phase, A determines the mass in the transverse 
spin-triplet g-propagator, just like the London penetra- 
tion depth determines the mass of the transverse photon 



Mn{k) 



1 



ro + ck^ + d\nn\/ i\k\ + GH\nr^ 



(4.7a) 



There also are mixed propagators {bq), but they do not enter 
any diagrams that are important for the phase transition. 
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in a superconductor, or the mass of the director fluctua- 
tions in a smectic-A phase, see Eq. (|3.6ll . The Ginzburg- 
Landau parameter, 



K — X/£_ cx \fulci\fc, 



(4.8c) 



is again independent of to or |r|. 

Let us now neglect the fluctuation of the magnetic or- 
der parameter, i.e., we put 



(4.9) 



Ml^ix) « 8a Sno m A/G VT. 



In the hmit A <C ^ this approximation becomes exact; in 
general, its validity will need to be investigated, q can 
then be integrated out, in exact analogy to the treatment 
of A in Sec. nil. A. 31 The result for the free energy density 
/, in a magnetic field h, and with /o = /(to = 0), is 

/ = /o + 7*0 TO^ + Mo to"' — hm 

-^^T^ln7V(fc,n„;m). (4.10a) 



fe<A 



A is an ultraviolet momentum cutoff, and 
iV(fe,r2„;TO) = IGclC^K^m^nl 
+ {\k\ + GHilnf \\k\ + G{H + kt)n„ 



(4.10b) 



Minimizing / with respect to the magnetization gives the 
equation of state. 

The integral in Eg . 14.10a|) has been analyzed by 
iKirkpatrick and Belitzl (2002.) . In d = 3, the equation 
of state was found to take the form 



h — 2rm + 4:V to ln(TO + T ) 
+ m-^ 4uo + 2v 



,2 _| rp2 

TO 

TO2+r2 



(4.11a) 



For the free energy density at T — h ^ this implies 
l|Belitz adll99d[l 



f = fo + r m + W3 TO In TO + um 
in d = 3, and 

f = fa + r — m'^^^ + u 



(4.11b) 



in generic dimensions. In these equations, /, to, and T 
are measured in terms of suitable microscopic quantities 
such that r, w^, and u are all dimensionless. r and u 
are given by tq and uq, respectively, plus additive renor- 
malizations from the soft modes, in analogy to the clas- 
sical example of renormalized mean-field theory in Sec. 
IIILA.3I Vd > IS quadratic in cf, so in strongly corre- 
lated systems Vd is larger than in weakly correlated ones. 
A comparison with Eq. (|3.8b|) shows that the free energy 
for the quantum ferromagnet in d = 3 is precisely analo- 
gous to that of the classical liquid crystal in d = 4 in the 
same approximation. 



This free energy functional is nonanalytic in to, in anal- 
ogy to Eqs. (|3.8|l in our classical example, and the same 
discussion applies. Namely, the nonanalyticity reflects 
GSI in the paramagnetic phase, and it is directly related 
to the nonanalytic wave number dependence of the spin 
susceptibility, Eq. (|2.36b|) . In real space, the latter cor- 
responds to 



Xs{r > 0, \x\ ^ oo) cx l/\x\ 



2d~l 



(4.12) 



This manifestation of GSI is the ultimate physical reason 
behind the failure of Hertz theory to correctly describe 
the quantum critical behavior. 

Note that the derivation of the above mean-field re- 
sult for the free energy, Eqs. (|4.11bl B.llc|l . implicity as- 
sumes that the nonanalytic wave number dependence of 
Xs given by Eq. Ij2.36bp or Eq. H4.12|l . which is exact in 
the paramagnetic phase, continues to hold in the critical 
region. This is a nontrivial assumption. In Sec. IIV.A.4I 
below we will discuss how critical fluctuations can mod- 
ify this result. This po int has recently been addressed by 
IChubukov enj] i2QQ?t) . 

The phase diagram predicted by these equations ind ~ 
3 is shown, in a more general context, in the first panel 
of Fig. |2ni below. There is a tricritical point at 



T = Tu 



-u/2v 



(4.13a) 



At T = 0, there is a first-order phase transition at r = ri, 
with the magnetization changing discontinuously from 
zero to a value toi. One finds 



TOl 



-(l+u/v)/2 



(4.13b) 



In d — 2, there is no finite-temperature magnetic phase 
transition. However, at zero temperature there is a QPT, 
which is predicted by the Eqs. (|4.10|) to be discontinuous. 
In d > 3 the nonanalytic terms produced by the soft 
modes are subleading, and the transition is described by 
ordinary mean-field theory. The generalized mean-field 
theory thus suggests an upper critical dimension d^ = 3. 
As we will see in the next subsection, a more sophisti- 
cated analysis confirms this result. 



(4.11c) 4. Effect of order-parameter fluctuations 



The renormalized mean-field theory predicts that the 
quantum ferromagnetic transition in clean systems is al- 
ways of first order. As we have discussed above, this 
conclusion is certainly valid in the limit A ^ ^. In gen- 
eral, however, the order-parameter fiuctuations need to 
be taken into account. This can be done by a systematic 
RG analysis of the action, Eq. (|4.1|) . In contrast to the 
liquid crystal case, it turns out that a one-loop calcu- 
lation predicts, under certain conditions, a critical fixed 
point that corresponds to a second-order transitio n. 

The fixed point found bv lKirkpatrick and Belit j (^2002^ 
has the property that G is marginal. To one-loop order. 
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in d = 3, and for r = 0, the coupling constants u, H, C2, 
and c obey the flow equations^^ 



du 
dlnb 
dH 

dlnb 
dc2 
dlnb 



-2u~Au cl/H, 



Ah/c 

-C2, 



dc 
dlnb 



(4.14a) 
-AJH, (4.14b) 
(4.14c) 



with b the RG length rescaling factor, and Ah, and 
Ac positive constants. In a purely perturbative treat- 
ment, C2 and H on the right-hand side of Eq. H4.14a|l are 
constants, and u inevitably becomes negative at large 
scales (6 ^ oo). A negative u is usually interpreted as 
signalizing a first-order transition;^'^ this is the RG ver- 
sion of the conclusion that the transition is always of 
first order. However, this perturbative argument is not 
consistent since, at the same level of the analysis, Eq. 
(|4.14b|) predicts a singular specific-heat coefficient H, 
which couples back into the coefficients of the Landau 
theory. This feedback is taken in to account by solving 
the one-loop equations H4.14|l self-consistently, whereby 
the one-loop correction to the w-fiow equation decreases 
with increasing scale, due a combination of C2 being ir- 
relevant and H increasing under renormalization. The 
conclusion that u flows to negative values is therefore no 
longer ine vitable. Indeed, a solution of the above flow 
equations ijKirkpatrick and Belitd . l2002(l shows that 



u{b — > oo) 



{Au/Ac){c2,Qf Co] r 



(4.15) 



with Mo = u{b ~ 1), etc., the bare coupling constants, 
K the Ginzburg-Landau parameter from Eq. I|4.8c|) . and 
Kc c>c Au/Ac- We see that the RG analysis pre- 
dicts an asymptotically negative value of u, and hence 
a fluctuation-induced first-order transition, for small val- 
ues of the Ginzburg-Landau parameter, k < Kc, in agree- 
ment with the renormalized mean-field theory, and in 
analogy with the liquid-crystal case. For k > Kc, how- 
ever, u stays positive and one finds a critical fixed point, 
in contrast to the liquid-crystal case. As mentioned in 
Sec. nil. A. 41 above, the mechanism for this "fluctuation- 
induced second- order transition" i s very similar to the 
one discussed by iFucito and Parisil l)l98lh for a classical 



At this level of detail of the discussion it is not obvious that 
there are two time scales in the problem, one related to the order- 
parameter fluctuations, and the other to the fermionic degrees 
of freedom. As a consequence, C2, which couples the two, does 
not have a unique scale dimension and can be either marginal or 
irrelevant, depending on the context. The flow equation 14.14cl 
applies to its irrelevant incarnation. We will discuss this point 
in more detail in Sec. IIV.B.4I below. 

While this interpretation of u flowing to negative values is not 
always correct, in the current case it is bolstered by the renor- 
malized mean-field theory. 
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FIG. 14 Phase d iagram of MnSi as measured by 
iPfleiderer et al\ lll997t) . The insets show the be havior of the 
susceptibility close to the transition point. From fVoita et gp 
L1999.V 



system. The critical behavior at this transition will be 
summarized in the next section. 

There is experimental evidence that the quantum fer- 
romagnetic transition in clean systems is of flrst order 
in some materials, but continuous in o thers, consistent 
with the theoretical picture given above. IPfleiderer et al\ 
(jfOO?) have reported a continuous transition in MnSi 
at moderate hydrostatic pressures corresponding to rel- 
atively high values of Tc, while the transition becomes 
first order if Tc is driven to low values by higher values 
of the pressure, see Figs . ITU and IT^ ^'^ Qualitatively the 
same magnetic phase diagram has been observed in UGe2 
l|Huxlev et adl200lD . see Fig.[ll In ZrZuj, on the other 
hand, the magnetic transition is of se cond order down 
to the lowest temperatures observed (jPfleiderer et alV 
The latter two materials are of particular inter- 
est since they also have a phase of coexistent ferromag- 
netism and superconductivity. The superconductivity is 
very vulnerable to disorder, which testifies to the fact 
that the samples in question are very clean. Another ma- 
terial where a continuous transitio n is observed down t o 
very low temperatures is Ni^^Pdi-^, l)Nicklas et QiJ .ll999). 
While the substitutional nature of this system inevitably 
introduces some disorder, the quantum phase transition 
occurs at low Ni concentration [x w 0.02), and disorder 



Strictly speaking, MnSi is not a ferromagnet, but rather a 
weak helimagnet. Local remnants of the helical order have been 
speculated to be responsible for the n on-Fermi-liquid proper- 
ties observed in the paramag netic phase ijDoiron-Levraud et aZl 
l2003t IPfleiderer et alV l2004li . The influence of the helical or- 
der on the quantum crit ical behavior has been studied by 
IVoita and Sknepneld J2001D . 
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FIG. 15 Magnetic susceptibility (in SI units) versus temper- 
ature for MnSi at pressures corresponding to the data points 
in Fig. 1141 Pressure values for the curves from right to left are, 
1.80, 3.80, 6.90, 8.60, 10.15, 11.25, 12.15, 13.45, 13.90, 14.45, 
15.20, 15.70, and 16.10 kbar. The change from a sharply 
peaked susceptibility consistent with a second order transition 
to a discontinuous one as expected for a first-order transition 
is apparent. From lPfleiderer et al\ ljl997t l. 



5. Critical behavior at the continuous transition 

a. Critical behavior in d — 3 A solution of the RG flow 
equat ions for the case k > Kc (jKirkpatr ick and Belitzl 
120021) yields a wave number dependent coefficient c in 
the paramagnon propagator, Eq. Ij4.7b|l . 

c(fe ^ 0) cx (lnl/|fe|)-i/26 ^ (4 ;Lg) 

where the exponent is determined by the ratio Ah /Ac = 
27. Expressing the logarithmic corrections to power-law 
scaling in terms of scale-dependent exponents, the crit- 
ical exponent rj, which describes the wave number de- 
pendence of the paramagnon propagator at criticality, is 
therefore given by 

77 = ^ lnlnVln& . (4.17a) 

The parameters r and d in the paramagnon propagator 
are not renormalized. The correlation length exponent 
1^, the susceptibility exponent 7, and the dynamical ex- 
ponent z can therefore be directly read off Eqs. (j4.7b|l . 
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FIG. 16 Magnetization versus hydrostatic pressure at T = 
2.3K for UGe2. The intial first-order transition at p ~ 15kbar 
is followed by a metamagneti c transition at p « 12kbar. After 
iPfieiderer and HuxlevI J2002D . 



is believed to play a minor role at the transition. 

While the soft-mode theory is consistent with these 
observations, it is not the only possible explanation. 
Another obvious possibility are band-structure effects, 
which can change the sign of the coefficient u in an 
ordinary Landau theory, Eq. (|4.10a|l without the last 
term, and thus lead to first-order transitions in some 
materials, and second-order transitions in others. In- 
deed, band structure effects have been proposed to be 
the source of the first-order transition, the superconduc- 
tivity, and an observed m etamagnetic transiti on wi thin 
the ferromagnet i c pha se l)Pfleiderer and HuxlevI l2002t 
ISandeman eit a/.U2003(l . 



i/=l/(2-r;) , z^3-r] , 7=! . (4.17b) 

The order parameter exponents (3 and 6 can be obtained 
from scaling arguments for the free energy. The result is 

P=l/2 , (5 = 3 . (4.17c) 

Finally, one can generalize the definition of the specific 
heat exponent a familiar from thermal phase transitions 
by defining Cv oc T~" at criticality. One finds 

a = -1 + (Inln6/ln6 - ?7)/z . (4.17d) 

The result for i] is valid to leading logarithmic accuracy; 
the values of 7, /3, and S, as well as the relations between 
T] and v, z, and a, respectively, are exact. 

The theory thus predicts the critical behavior in d = 3 
to be mean- field like with log arithmic corr ections. Al- 
though the fixed point found bv lHertd l)l976fl is unstable, 
it is so only marginally. Hertz's re sults, and the ir exten- 
sion to nonzero temperatures bv lMilli3 (^2^, should 
therefore apply apart from corrections that are too small 
to be detectable with current experimental accuracies. 

Although the critical behavior has not been probed 
directly experimentally, a combination of various expo- 
nents determines the shape of the phase diagram at low 
temperatures. To see this, consider a homogeneity law 
for the magnetization, which can be obtained by differ- 
entiating Eq. p. 10(1 twice with respect to h, and putting 
/i = 0: 

Tn{r,T,u) = &-''-^+2^^'•TO(r6^/'',T6^^i6["l). (4.18a) 

Here we have included the coefficient u of the quartic 
term in the free energy, even though its scale dimension 
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FIG. 17 Phase diagram of Nia;Pdi_a;. Equation is 
obeye d over an x-ra ngc of two decades close to the transition. 
From lNicklas et al\ IHQQal . 
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FIG. 18 Phase diagram of MnSi. These are the same data as 
in Fig. 1141 scaled to display the relation given in Eq. 14.191 . 
The tricritical point separating second and first order transi- 
tions coincides with the poin t where the scaling breaks down. 
After iPfleiderer et al\ ^ 



[w] < is negative, and u is thus irrelevant. The reason is 
that u is a danger ous irrelevant variable w ith respect to 
the magnetization l)Fisheil. Il983 iMaL Il976fl . which influ- 
ences the critical behavior of m. To see this, consider the 
equation of state, Eq. (|4.11a|) . In the absence of the soft- 
mode corrections, one has to (x \J —r ju. m thus depends 
on r/u, and the homogeneity law can be written 

m{r/u,T) = 6-''-"+22"'m((r/ii)6i/''"["l,r6^) 

= {r/uf TO (l, T (r/w)-'^^/(i-'^["l)) , (4.18b) 

where /3 — v{d+ z — 2yh) /{l — v[u]). The relation between 
Tc and r is now obtained from the requirement that to 
must have a zero for T = Tc. This yiel ds r — const, x 
j.(i-i,iu])/^z ilMilliaL Il993t ISachdevl \l99^ . [u] =A-d-z, 
as can be seen from Eq. (|4.2|l by power counting, and 



neglecting the lo garithmic cor rections to scaling we find 

(4.19) 



3 garithmic cor r 
from Eqs. (|CT7|l l|Millisl Il993|) 



This agrees well with the phase diagram measured in 
Nij^Pdi-a;, and also with the portion of the phase di- 
agram in MnSi where the transition is continuous, see 
Figs.[T7|[THl 



b. Critical behavior in d ^ 3 In d > 3, the RG analysis 
shows that Hertz's fixed point is stable, and all exponents 
have their mean-field values. This confirms the sugges- 
tion of the renormalized mean- field theory that the upper 
critical dimension is d+ = 3. 

In d < 3, the critical behavior can be studied by means 
of an expansion in e = 3 — d. KirkDatrick and Bc lit3 
(12002') have found a fixed point where G, if, c, and ci are 
marginal. C2 can again be either marginal or irrelevant. 



depending on the context. A one-loop calculation of the 
exponent 77 and the specific heat exponent a yields 

77 = -e/26 , a = -d/{3^r]). (4.20) 

i^, z, 7, (3, and 6 are still given by Eqs. (|4.17bl ICTfcll . 



B. Quantum Ferromagnetic Transition in Disordered 
Systems 

If one adds quenched disorder to the problem con- 
sidered in the previous subsection, the action changes 
relatively little. The nature of the generic soft modes 
changes; they are now diffusive rather than ballistic in 
nature, but this change by itself will clearly only change 
the upper critical dimensionality. Apart from this, the 
disorder needs to be averaged over, e.g. by means of the 
replica trick. It turns out that the disorder average leads 
to some terms in perturbation theory having a sign that 
is opposite from the corresponding one in the clean case. 
For instance, the one-loop correction to the quartic cou- 
pling constant u is positive, and as a result the phase 
transition is always of second orde r. Treating t he generic 
soft modes in tree approximation. Ulertj l)l976() had con- 
cluded that the critical behavior is mean-field like for all 
d > 0. This turns out to not be true, for reasons anal- 
ogous to those that invalidate Hertz theory in the clean 
case. It turns out, however, that the critical behavior 
can still be determined exactly for all dimensions d > 2, 
although it is not mean-field like. Historically, this ex- 
tension of Hertz's theory for the disordered case was de- 
veloped earlier than the corresponding treatment of the 
clean case (Bclitz et al, 2001a b; Kirkpatrick and Bclit^ 
119961) . Here we explain the structure of this theory and 
summarize its results. 
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1. Soft-mode action 



The action as given in Sec . II V. A . 1 1 remains vahd, with 
two modifications. First, all fields carry replica indices 
that need to be summed over. Specifically, M carries 
one replica index, while Q and q carry two replica in- 
dices, as explained in connection with Eq. H2.44|l . The 
symmetrized magnetization field b is diagonal in its two 
replica indices. Second, the fermionic vertex function 
r(^) changes; Eqs. H4.3bp and H4.3cll get replaced by 



*ri2,34(fe) — (513 (524 



x27tTG d,oK, + il-S,o)Kt . (4.21a) 



with replica indices ai, and 

.(2,0) 



1^12' (^) — fe + GHQjij^^n^ 



(4.21b) 



All other terms in the action remain unchanged, except 
for the addition of replica indices. 

It is obvious from simple physical consideration that 
the quenched disorder must lead to additional terms in 
the action. For instance, consider the bare distance from 
the critical point, r in Eq. H4.2|l . a random function of 
space, and integrate out this 'random mass' with respect 
to some distribution function. This will generate terms of 
higher order in M, starting at 0{M'^), that have a differ- 
ent imaginary-time structure than the u M'^ term in Eq. 
(|4.2|) . Such random-mass terms do indeed get generated 
by the basic disorder term in an underlying microscopic 
action, Eq. H2.37c|) or H2.42|l . and thus need to be added 
to the effective action. However, they turn out to be ir- 
relevant for the critical behavior in all dimensions except 
d — 4:, and we therefore neglect them. 



2. Gaussian approximation 

Keeping only terms that are bilinear in the fields yields 
the Gaussian approximation, as in the clean case. The 
paramagnon propagator now reads 



Mnik) 



1 



ro + ck^ + d\nn\/ {k^ + GH\n„\) 



(4.22) 



Formally this is the same as Eq. H4.7a|l with |fc| replaced 
by fe^, but the coefficient G has a different physical in- 
terpretation as was explained in connection with Eqs. 
(|2.48(l . The basic fermionic propagator is given by Eq. 
(I2.48b|) . As in the clean case, there is a mixed propagator 
that turns out not to be important for the determination 
of the critical behavior. 



3. Renormalized mean-field theory 

A renormalized mean-field theory can be constructed 
in exact analogy to the treatment of the clean case in 



Sec. IIV.A.3I ("S essions and Belitil200^ . The free en ergy 
in this approximation is again given by Eq. H4.1Ua|l . but 
the quantity N in Eq. (|4.10bp is replaced by 



fc2 + G{H + kt)n. 



320Fc2G2m2 



(fe2 + GHVLn) + 32V^C2G2to2 

(4.23) 

An analysis of the integral in Eq. (|4.1Ua|l with this ex- 
pression for N yields a free energy at T = of the form 



f = fo + rm 



Ad+2)/2 



+ um^-h, (4.24a) 



or, equivalcntly, an equation of state 

h — r m + w m'^^^ + u m'^. 



(4.24b) 



In Eqs. (|4.24() . w > is a positive coefficient that is 
proportional to the disorder, and r and u again represent 
additive renormalizations of vq and uq, respectively. 

Regarding the nonanalytic dependence of / on m, the 
same comments apply as in the clean case. Sec. IIV.A.3I 
and an assumption analogous to the one mentioned after 
Eq. H4.12|l has been made. As far as the GSI is concerned, 
the only difference is that the decay of the spin suscepti- 
bility in real space far from the transition is given, instead 
of Eq. g32l, by 



Xs{r > 0, |a;| 00) cx l/\x\ 



2d-2 



(4.25) 



which is the Fourier transform of Eq. (j2.35d|l . Also, 
the magnetization now scales as a wave number squared, 
which explains the exponent of the nonanalyticity in Eq. 

Despite these similarities, there is a crucial sign dif- 
ference between the nonanalytic term in the clean case, 
Eqs. (|4.11bl I4.11c(l . and in the one in Eq. (|4.24(l . Al- 
though this is not of much significance from a GSI point 
of view, it leads to the quantum phase transition to be 
of first order in the clean case (at least at the level of 
the renormalized mean-field theory), while it is of sec- 
ond order in the disordered case. Physically, this can be 
understood by means of arguments very similar to those 
that explained the sign difference between Eqs. (|2.35dp 
and Ij2.36b|l . respectively. In the clean case, the nonan- 
alyticity is produced by fluctuation effects, represented 
by the generic soft modes, that weaken the tendency to- 
wards ferromagnetism. Accordingly, the constant contri- 
bution from the integral over InA^ in Eq. (|4.1Ua|) . which 
changes rg to r, is positive. The appearance of a nonzero 
magnetization gives the soft modes a mass, and hence 
weakens these fluctuations, so the leading m-dependent 
correction to the constant contribution is negative. In 
the disordered case, on the other hand, the quenched 
disorder leads to diffusive motion of the electrons, which 
is much slower than the ballistic dynamics in the clean 
case. This increases the effective interaction strength be- 
tween the electrons, which in turn favors ferromagnetism. 
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FIG. 19 Phase diagram in the G-r plane resuhing from Eq. 
ll^:^ at T = for = 1, ?; = 0.5, A = 0.5, Bt = Bo ^ 
1, showing a second-order transition (dashed hne), a first- 
order transition (solid line) ^ a critical end poi nt (CEP), and 
a critical point (CP). From lBelitz et al\ lll999t) . 



Consequently, r is smaller than tq. A nonzero magnetiza- 
tion again weakens this effect, and therefore the leading 
nonanalytic m-dependence of the free energy is positive. 

A more general form of the free energy in renormahzed 
mean-field approximation, which is valid for T > and 
allows the zero-disorder limit to be taken, has been con- 
sidered bv lBehtz et all l)l999fl . The free energy density 
discussed by these authors takes the form 

/ - fo + rm' + iG/eF)Am* [m^ + (StT)^] 

+v In [m^ + {T + BcGf] +um^ + Oim^). 

(4.26) 

Here Bt and Bq are dimensionless parameters that mea- 
sure the relative strengths of the temperature and the 
disorder dependence, respectively, in the two nonanalytic 
terms, and A is a measure of how strongly correlated the 
system is. They are all expected to be of order unity. 
For G — the resulting equation of state reduces to Eq. 
(|4.11a|l . For sufhcicntly large disorder G, the logarith- 
mic term is unimportant, and at T = one recovers 
Eqs. H4.24fl . This form of the free energy thus interpo- 
lates between the clean and disordered cases. It displays 
a rich phenomenology, as shown in Figs. E| [201 Notice 
that, for a disorder larger than a threshold value, they 
predict a second magnetic transition inside the ferromag- 
netic phase that is always of first order. 

The critical exponents f3 = 2/{d — 2) and S = d/2 (for 
2 < d < 6) can be read off Eq. H4.24b|l . but a determi- 
nation of the remaining exponents requires the consid- 
eration of order-parameter fluctuations. We will discuss 
this in the next subsection, where we will also see that 
the critical behavior predicted by the renormahzed mean- 
field theory is very close to the exact one. 



4. Effect of order-parameter fluctuations 

The effects of order-parameter fluctuations on the dis- 
ordered ferromagnetic quantu m phase tran sitio n have 
been studied at various levels. iHertd il97(t) and lMilhsl 



l|1993|) considered order-parameter fluctuations, but inte- 
grated out the generic soft modes in a tree approximation 
that n eglects all fermionic loops. iKirkpatrick and Belit j 
(|l99ff) kept fermionic loops, but still integrated out the 
fermions, which led to a nonlocal field theory in terms of 
the order paramete r only, which they ana lyzed by means 
of power counting. iBelitz et all |2001a"b') kept all of the 
soft modes explicitly, and on equal footing, in analogy to 
the theory for clean systems discussed in Sec. IIV.AI and 
to the classical theories covered in Sec. Illll This theory 
contains the previous ones as approximations, and we 
will sketch it in the remainder of this section. 



a. Coupled Geld theory, and power counting In contrast to 
the clean case, Sec. IIV.B.II for disordered magnets the 
fermionic part of the soft-mode action is known in closed 
form, namely, the nonlinear a model, Eqs. (|2.61l [TH^ . 
with no spin-triplet interaction. The order-parameter 
part of the action is again given by Eq. (|4.2|l with a 
replicated M field as explained in Sec. IIV.B.II and the 
coupling between the two is given by Eq. H4.4a|l . with Q 
replaced by Q and replica indices added as appropriate. 
This action can be written down on general principles, or 
it can be derived from the complete nonlinear a model, 
Eqs. (12.6111. by per forming a Hubbard-Stratonovich de- 
couphng (^Hubbardlil959,: Stratonovich. .1957.) of the spin- 
triplet interaction and adding an term.^^ Since the 
nonlinear a model has been derived from a microscopic 
action, this also provides a derivation of the final effec- 
tive action from a microscopic starting point, if that is 
desired. 

We now expand this effective action in powers of M 
and q. In a schematic notation that omits everything not 
necessary for power counting, we have 



Acs[M,q] =AM + Ag + A, 



M,q- 



(4.27a) 



with 



Am = 



- J dx M [r + ca-2 + cdl] M 

+0{diM^,M^), (4.27b) 

= J dx{d^qf+H j dxnq^ + KsT J dx q^ 

--^ j dxdlq^ + Hi j dx^q^ 

+0{Tq\dlq\^q''), (4.27c) 

AM,q = VTci j dxMq + \/Tc2 j dx M q^ 

+0{VTMq^) . (4.27d) 



The quartic and higher terms in M originate from massive 
modes in the underlying miscroscopic action, which have been 
dropped in deriving the nonUnear cr model. 
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FIG. 20 Phase diagrams in the T-r plane resulting from Eq. llOSt iov u = Ba ^ 1, v = Bt = A ^ 0.5. TC P and TCEP 
denot e a tricritical point and a tricritical endpoint, respectively. All other symbols are the same as in Fig.EHl From lBelitz' etall 



Here the fields are understood to be functions of position, 
frequencies, replica labels. Only quantities that carry a 
scale dimension are shown. Accordingly, frequency and 
replica sums have been omitted, but appropriate powers 
of the temperature are shown. All of the terms except 
the second one in Am can be obtained by combining Eqs. 
(|0 I13al lO ESI and Eqs. g^U- The bare value of 
the coefficient Cd-2 is zero, but we will see that such a 
term is generated under renormalization. 

For a power-counting analysis of this effective action, 
we assign scale dimensions [L] — —1 and [T] = [fl] = —z 
to lengths and temperatures or frequencies, respectively. 
It is important to note that z is unlikely to have a unique 
value: Since the critical paramagnon propagator will in 
general have a frequency-wave number relation that is 
different from that of the fermionic q propagator, we ex- 
pect at least two different time scales in the problem, 
hence two different values of z. This needs to be taken 
in to account in the power counting analysis. 




b. Her tz's fixed point The fixed point identified bv lHertj 
l)l976j) is recovered from this formalism by the coefficients 
c and ci to be dimensionless. The field q is expected to 
be diffusive, and we therefore choose its scale dimension 
to be [q] — {d— 2)/2, in accord with Eqs. H2.63|l . From 
Am we then obtain [M] = (d - 2)/2 and [t] = 2. r 
is thus relevant, as expected, and the correlation length 
exponent is = l/[r] = 1/2. Cd-2 is also relevant, of 
course, but since its bare value is zero we ignore it for 
now. The first term in AM,q produces the dynamical part 
of the paramagnon, as was demonstrated in Sec. lIV.iOl 
The scale dimension of the y/T prefactor thus yields the 
paramagnon or critical dynamical exponent z^ — 4. The 
frequency and the temperature in the second and third 
term, respectively, in Aq carry the fermionic time scale 
Zq = 2 that is consistent with diffusion. G, H, and 
are then all marginal, and all higher order terms in Aq 
are irrelevant. For the scale dimension of C2 one finds 



FIG. 21 One-loop renormalization of the M'^ vertex. Dashed 
lines denote M fields, solid lines denote q propagators, and 
the vertices carry a factor of C2 each. 



Due to the existence of two time scales mentioned above, 
one needs to distinguish between two different variants of 
C2, which differ with respect to the value of z that enters 
their scale dimension. With z = Zc = 4, C2 = is irrel- 
evant for all d > 2. However, with z = = 2, C2 = is 
relevant with respect to the putative fixed point for d < A. 
This instability of Hertz's fixed point is indeed realized. 
Consider, for instance, the renormalization of the M"^ 
vertex by means of the diagram shown in Fig. 1211 Since 
it is a pure fermion loop, the factors of ^/T which come 
with the vertices, and which are absorbed into the loop 
integral over the frequency, carry indeed the fermionic 
time scale, and the relevant coupling constant is cj. We 
note that this is the same perturbative contribution that 
leads to the nonanalytic wave number dependence of the 
spin susceptibility deep in the paramagnetic phase, Eq. 
(|2.35d|l . The term with coupling constant Cd-2 is thus 
generated from the C2 coupling. It is convenient to ex- 
plicitly add this term to the bare action, as we have done 
in Eq. Ij4.27bp . although the physics it represents is al- 
ready contained in the C2 coupling term. Indeed, with 
respect to Hertz's fixed point, (cj)^ and Cd~2 have the 
same scale dimension, [cd-2] = 2[cJ] = 4 — c?. 

One thus finds that the generic soft modes render 
Hertz's fixed point unstable, and this first becomes ap- 
parent at one-loop order. This instability, although 
not its source, can also be deduced from the fact that 
the mean-field value for the correlation length exponent 
V = 1/2, violates, for all d < 4, the rigorous bound pro- 
vided by the Harris criterion ijChaves et al ]. ll98(il:lHarri^ 
Il974f) for systems with quenched disorder, v > 2/d. 



~-{d+z-Q)/2. 



(4.28) 



42 



c. A marginally unstable fixed point The behavior of the 
renormahzed mean-field theory is recovered from the 
RG if one (1) takes the Cd-2 term into account, and 
(2) drops the requirement that C2 be marginal. We 
thus require only that ci be marginal, which implies 
[M] = I + [d - z) /2, and that [q] = {d-2) /2 and Zq ^ 2 
as before. At this point it is useful to consider an anal- 
ogy between the present problem and that of a classical 
fluid under shear discussed in Sec. IIII.R21 In the latter 
case, long-ranged correlations between order-parameter 
fluctuations stabilized mean-field critical behavior in di- 
mensions lower than the usual upper critical dimension. 
In the current problem, the Cd-2d^~^ term in the action 
represents long-range order parameter correlations that 
decay as |a;|~^'^'^~^^. By analogy with the classical fluid 
example, it is therefore natural to expect the stabilization 
of a Gaussian fixed point where Cd-2 is marginal in the 
critical paramagnon. This implies [M] = 1, and hence a 
critical time scale and an exponent rj characterized by 

z^ = A-7j = d. (4.29a) 

r is relevant with [r] = d — 2, which corresponds to a 
correlation length exponent 

jy = l/(d-2). (4.29b) 

In contrast to the situation at Hertz's fixed point, this 
respects the Harris criterion. 

Equation H4.29a|l implies that the magnetic suscepti- 
bility at criticality decays in real space according to 

Xs(r = 0,|a;|->c3o)cxl/|a;|2. (4.30) 

In the relevant dimensionality range, d > 2, these correla- 
tions are of longer range than the power-law correlations 
in the paramagnetic phase, Eq. H4.25|l . We see that, as in 
the case of the classical fluid, Sec. IIII.Rl2l the long-ranged 
correlations or GSI in the disordered phase influence the 
critical behavior, and lead to correlations with an even 
longer range at the critical point. 

The results given above hold for 2 < d < 4. For d > 4, 
77 = and V = 1/2 have their mean-field values, and 
Zc — 4. The exponents P and 6 can be ob tained either 
from s caling arguments for the free energy ijBelitz et all 
I2OOIIJ) . or from repeating th e above counting argument s 
in the ferromagnetic phase ijSessions and BelitzL l2003|) . 
By either method one finds the same result as from Eq. 
llt24hll . viz., 

(3^2/{d-2) , S^d/2. (4.31) 

These relations hold for 2 < c? < 6, while for d > 6 one 
has the mean-field values (3 — 1/2, 5 = 3. A detailed 
analysis reveals that Hertz's fixed point is actually stable 
for all d > 4, but for 4 < d < 6 the coefficient w in 
Eq. I|4.24b|l acts as a dangerous irrelevant variable with 
respect to the magnetization, which explains why P and 
6 lock into their mean-field values only for d > 6. 



What remains to be done is to investigate the stability 
of the Gaussian fixed point with the critical behavior de- 
scribed above. It turns out that it is marginally unstable, 
and that the exact critical behavior is given by the power 
laws given above with additional logarithmic corrections 
to scaling. We give these results next. 



d. Exact critical behavior The results of the previous 
subsection characterize a Gaussi an fixed point that was 
descri bed by the theory given by iKirkpatrick and Belitj 
|)1996|) . To check its stability, one needs to consider the 
scale dimensions of the remaining terms in the action. C2 
has a scale dimension 

[C2] - 1 - z/2 (4.32) 

with respect to the Gaussian fixed point, which implies 
that C2 (i.e., C2 with z = Zq) is marginal. This was to be 
expected, since Cd-2^ which describes the same physics 
as C2, is also marginal. However, it implies that the non- 
Gaussian C2 term must be kept as part of the fixed-point 
action. All other terms are nominally irrelevant. How- 
ever, it turns out that, for subtle reasons again related to 
the existence of multiple time scales, the terms of order 
in Aq can also be effectively marginal and must be 
kept as part of the fixed point action. For all d > 2, the 
latter then consists of the terms shown explicitly in Eqs. 

Since the fixed point action is not Gaussian, the ex- 
act critical behavior cannot be obtained simply by power 
counting. However, it turns out that the problem can 
still be solved exactly, by me ans of an infinite re sum- 
mation of perturbation theory ijBelitz et 0^112001^ . For 
2 < d < 4, the result consists of logarithmic corrections 
to the exponents derived in the previous subsection. For 
instance, the correlation length depends on r via 

^ cx e-™"''*-^('"''^(i/''»', (4.33) 

to leading logarithmic accuracy. Notice that the correc- 
tion to the power law varies more slowly than any power 
of r, but faster than any power of Inr. This behavior 
is conveniently expressed in terms of exponents that are 
scale dependent via a function g whose asymptotic be- 
havior for large arguments is 

g{x » 1) « [21n(d/2)/7r]-l/" g[ln(c(<J).)]V21n(d/2)_ 

(4.34) 



A more conventional approach would be to not include the Cd—2 
term in the action, but rather have its effects included in a renor- 
malization of c while working in d = 4 — e. This can be done, 
but the loop expansion in this case does not lead to a control- 
lable expansion in powers of e since C2 remains strictly marginal 
order by order. N evertheless, such a proce dure yields interesting 
technical insig hts iRoUbuhler et a;.ll20ol) . 
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The dimensionality dependent coefficient c{d) goes to 
zero for d ^ A, and to a constant for d ^ 2. In the 
case of ly, Eq. (|4.33|l corresponds to 



l/i' = d - 2 + \n g{\nb)/ Inb. 
For the other exponents one obtains 

Zc = 4: - T] = d + lng{lnb)/ \nb, 

S = ~a/2d^ zj2, 

(3 = 2iy , 7=1, 

Zq = 2 + In.g(ln6)/ln6, 



(4.35a) 



(4.35b) 
(4.35c) 
(4.35d) 
(4.35e) 



where the specific heat exponent a is defined as in the 
context of Eq. Ij4.17d|l . 

The critical behavior of various transport coefficients 
and relaxation rates, as well as the tunneling; density of 
state s, have also been determined (feelitz et al) . l2001bl 
l200f]l) . In particular, these authors showed that the quasi- 
particle prop erties at criticality are those of a marginal 
Fermi liquid f Varma et al\ .\l98^). Here we just quote the 
result for the electrical conductivity. In d = 3 at critical- 
ity, one finds a nonanalytic temperature dependence in 
addition to a noncritical background term. 



cr(T ^ 0) = cr(T = 0) + const. X 



Tg{- HeF/T) 



-,1/3 



(4.36) 

with g{x) from Eq. H4.34|l . 

Experimentally, there are few systematic studies of 
the infiucncc of quenched disorder on the ferromagnetic 
quantum phase transition. URu2-2;Rea:Si2 shows a rich 
phase diagram in the T-x plane, with phases display- 
ing ferromagnetic, antiferromagnetic, and superconduct- 
ing order in addition to a paramagnetic non-Fermi-liquid 
("strange metal") phase llBaueill2002HBauer etai\ . \200A 
iDalichaouch et all Il989t ISteward. l200l|) . see Fig. 1^ 
Near x = 0.3 there is a quantum phase transition from 
the ferromagnetic phase to the strang e-meta l phase that 
has been st u died i n some detail by iBaued l)2002() and 
iBauer et all l)2003(l . These authors found a value of 
S = 1.5, in agreement with Eq. H4.35c() (ignoring the log- 
arithmic corrections to scaling). The specific heat coef- 
ficient for T ^ is found to diverge logarithmically, or 
as a very small power of T, which is, within the exper- 
imental accuracy, consistent with Eq. (|4.35c|l . However, 
the logarithmic divergence of the specific heat coefficient 
is observed to persist away from the quantum critical 
point in both the paramagnetic and the ferromagnetic 
phases. Similar non-Fermi liquid properties have been 
found in the paramagnetic phase of other materials that 
display a quantum ferromagnetic transition, e.g., in MnSi 
|Pfleiderer et ai, 2001a). We will come back to these ob- 
se rvations in Se es. IIV.UI and IVI 

iDiTusa et al\ |2003l|l have investigated Fei_a;Coa;S2, 
which displays a T = insulator-metal transition at a 
very small value of a; (a; < 0.001), followed by a quantum 
paramagnet-to-ferromagnet transition at a; « 0.032. The 
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FIG. 22 Phase diagram of URuReSi, showing ferromagnetic 
(FM), antiferromagnetic (AFM), superconducting (SC), and 
non-Fermi liquid (NFL) phases. The different symbols refer 
to dif ferent metho ds for determining the phase boundaries. 
After iBaueJ ll2002D . 



strength of the quenched disorder is substantial in this 
system, with values of k^i ranging from 2 to 15 for x 
between 0.001 and 0.17. The authors have reported scal- 
ing of the conductivity near the ferromagnetic transition 
that is consistent with Eq. H4.36|l . 



5. Effects of rare regions on the phase transition for Ising 
system 

We now turn to rare regions and nonperturbative dis- 
order effects. In contrast to the previous subsections, 
which considered Heisenberg ferromagnets, we will now 
discuss a special effect for magnetic order parameters 
with an Ising symmetry. In such systems, the rare re- 
gion effects are much stronger than the quantum Grif- 
fiths phenomena discussed in Sec. III. B. 61 because of the 
coupling between the order parameter and the generic 
soft modes. To see this, consider a particular rare region 
devoid of impurities. Within such a region, the appropri- 
ate paramagnon propagator is given by Eq. H4.7a|l . and 
the coupling between the magnetization and the conduc- 
tion electrons is reflected by the dependence on |r2„|.^^ In 
imaginary time space, this corresponds to a long-time tail 
given by 1/r^ , see Sec. lII.B^Tl At T = 0, every rare region 
thus maps onto a clean, classical ((i-l-l)-dimensional Ising 
model that is finite in d dimensions and infinite in one 
dimension. The interactions between the spins are short- 
ranged in the former, and long-ranged, proportional to 
1/t^, in the latter. 

This long-range interaction can have drastic conse- 
quences, namely, the sharp quantum phase transition can 



^'^ This is sometimes referred to as Landau damping, in analogy 
to the fate of the plasmon when it enters the particle-hole con- 
tinuum. In the absence of such a coupling, e.g., in the pure spin 
model given by Eq. i2'75l . one would have a propagating mode 
with a frequency dependence on H^. 
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be destroyed by smearing. This can be understood as 
follows. A one-dimensional Ising model with a 1/rMn- 
teraction is known to possess an ordered phase l)Cardvl 
ll98lHThouleslll969D . A rare region can therefore de- 
velop a static order parameter independently from the 
rest of the s ystem. This result i s consistent with the one 
obtained bv lMillis et all ()200 2atl. who directly calculated 
the tunneling rate of a Griffiths island in an itinerant 
system and found it to vanish for a sufficiently large is- 
land. Suppose a nonzero average order parameter first 
appears on a rare region, rather than in the bulk of the 
system. Since the order is truly static, the ordered region 
will effectively act like a small permanent magnet imbed- 
ded in the system, and it will be energetically favorable 
for subsequent rare regions to align their order parameter 
with the first one (assuming that the effective interaction 
between the regions is also ferromagnetic). As a result, 
a finite magnetization appears as soon as a finite vol- 
ume of rare regions start to order. This magnetization is 
exponentially small, and the correlation length does not 
diverge at this point. However, the ordered rare regions 
provide an effective magnetic field seen by the bulk of 
the system; therefore, the l atter cannot u ndergo a sharp 
quantum phase transition ijVoital . l2003a(l . Notice that 
the number of ordered regions changes continuously with 
the control parameter, so the smeared transition is not 
simply given by the behavior of electrons in a fixed mag- 
netic field. Rather, order develops in different parts of 
the system at different values of the control parameter, 
and the order parameter close to the smeared transition 
is very inhomogeneous in space. The same mechanism 
also destroys classical phase transitions in systems with 
planar defects (Sknconck and Voita, 2003: Voita, 2003b). 

The above arguments are valid at T = 0, where the sys- 
tem is infinite in imaginary time direction. The behavior 
at finite temperatures is less well established. Since the 
rare regions are far apart, their interaction is very small. 
Therefore, the relative alignment between the order pa- 
rameter on different rare regions vanishes already at a 
temperature that is exponentially small in the density of 
the rare regions, i.e., double-exponentially small in the 
disorder strength. Above this temperature, the rare re- 
gions act as independent classical moments. Whether or 
not there is a second crossover at even higher tempera- 
tures to quantum Griffiths behavior similar to that in un- 
damped systems (as discussed in Sec. Ill.fj.6|) is not fully 
understood; it appears t o be a question of the micro- 
scopic parameter v alues ijCastro Neto and Jonesl l2000t 
iMillis et oiJ . l2002aD . 

In systems with continuous order parameter symme- 
try the rare region effects are weaker. Specifically, the 
quantum phase transition will remain sharp because the 
rare regions cannot develop static order. This can be 
seen by mapping each rare region onto a classical one- 
dimensional XY or Heisenberg model with inter- 
action. In contrast to the corresponding Ising model, 
these models do not have an ordered phase l|Brund . f200lt 
lKosterlit3. 1 1 976fl . The quantum Griffiths behavior in the 
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FIG. 23 Superconducting versus normal-state resistivity 
for three superconductin g materials. The data are taken from 
iRoweLL and DvnesI (llQSOTl : the lines are fits to a theory for the 
Tc degradation that is not important for the present discus- 
sion. From .Belitz and Kirkpatrick (.1994') . 



vicinity of the transition has not yet been worked out in 
detail, but it is likely to be of a more conventional type. 



C. Metal-Superconductor Transition 

Sufficiently strong quenched nonmagnetic disor- 
der^* decreases the critical temperature Tc for the 
superconductor-metal transition in bulk conventional su- 
perconductors,^^ see Fig. ESI At a critical value of 
the disorder, Tc vanishes, and there is a quantum phase 
transition from a normal-metal phase to a supercon- 
ducting phase where the concepts of additional soft 
modes and GSI effects play a crucial role. Since the 
generic soft modes in question, particle-hole excitations 
again, are massive at nonzero temperature, this quan- 
tum phase transition requires a theoretical description 
that is qualitatively different from BCS theory and the 
theories that describe fluctuation corrections to it. Such 
a theory, whose main predictions should not be hard to 
check experimentally, is presented below. An alterna- 
tiv e, and physically different, theory has been proposed 
bv 'Vishve shwara et al\ l)200C() . These authors have pro- 
posed that the normal-metal-superconductor transition 
triggered by disorder leads generically to a gapless su- 
perconducting state, and they assert that the resulting 
gapless quasiparticle excitation spectrum will have conse- 



Weak disorder can actually enhance Tc , sometimes substantially 
so. Aluminum, and many other low-Tc superconductors, are ex- 
amples for this effect. 

Also very interesting are the analogous effects in thin super- 
conducting films, where nonmagnetic disorder leads to a transi- 
tion from a superconducting phase to an insulating one, with or 
without an intermediate metallic phase. This topic is outside the 
scope of this review. 
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quences for the critical behavior. It remains to be worked 
out to what extent this assertion is correct, or to what 
extent the gaplessness is a property of the stable fixed 
point that describes the superconducting phase, rather 
than the critical fixed point. We will further comment 
on this in Sec. IV.B.2.bl below. 

One can construct a theory for the metal- 
superconductor transition that is structurally very simi- 
lar to the one for the disordered f erromagnetic transition 
discussed in the previous section ijZhou and Kirkpatrickl 
120041) . so we will keep the discussion brief. There are, 
however, a few crucial differences. The most important 
one is that the superconducting order parameter field 
^{x), whose expectation value is closely related to the 
anomalous Green function and the superconducting gap 
function, is related to the bilinear product of fermion 
fields at equal and opposite frequencies, 



(4.37) 



which in turn can be expressed in terms of the matrix 
elements q of the matrix Q, see Eqs. (|2.44l IT^^ . Since 
the q are soft modes, see Eq. (I2.48f) . this means that the 
coupling of the order parameter to the generic soft modes 
that is even stronger than in the ferromagnetic case, and 
therefore the effects of GSI are even more dramatic. 
We write the action as. 



= Am, + Aq + Am,, 



(4.38) 



Am, is analogous to Am, Eq. H4.2|l . It is a static, local, 
LGW functional for a three-dimensional XY-model. The 
order parameter field one can take to be complex val- 
ued, or as having two components '^^^{x), with r = 1,2. 
Structurally it is identical to Eq. H4.2() with Ml^{x) re- 
placed by \E'^(a;). Since we are dealing with a disordered 
system, ^ also carries a replica label. The fermionic part 
of the action, Aq, is identical to the one in the previous 
section, except now the particle-particle or Cooperon de- 
grees of freedom, r = 1,2, need to be taken into account. 
In particular, the Gaussian part of Aq is given by Eq. 
^ with r given by Eqs. (|4.21l) for r = 0, 3 and by 



2ttTGKc (4.39) 

for r = 1,2. Here r(^'°' is given by Eq. (j4.21b|l . and is 
a repulsive Cooper channel interaction that is generated 
by the particle-hole channel electron-electron interactions 
Kg and Kt in analogy to the mechanism that generated 
Kt in Eq. (j4.3b|l . It is also disorder dependent. As we will 
see below, Kc is responsible for driving to zero. The 
bare attractive Cooper channel interaction that is due to, 
e.g., phonon exchange, is denoted by — — |-fl^c| < 0. 
The coupling between 4" and q originates from a term 
Am,-q that can be written, in analogy to Eq. (|4.4d|l . 



The functional form of this term becomes plausible if 
one realizes that the order parameter field \1/ acts like an 
external field that couples to the particle-particle number 
density. To bilinear order this yields a contribution to 
AM,^q given by 

AM,-q ^ -SciVTY^ dxY^ r(3i2{x)°q^^{x). 
12 r=l,2 

(4.40b) 

Here ci oc li^d^/^, and 

rf3i2{x) = 5aia2 X^^^.^i+^^^n^'l^c). (4.40c) 
n 

This structure is in exact analogy to Eqs. (|4.41 14.51) for 
the magnetic case, only the factor of i in Eq. H4.40a|l . 
the sign in Eq. (|4.40b|) . and the frequency structure in 
Eq. H4.40c() reflect the particle-particle channel. Again in 
analogy with the magnetic case, one can further expand 
A^/^q in powers of q. The next term in this expansion has 
an overall structure 

AM,-q2 cx C2VT J dx tr {f3{x) q{x) q''ix)). (4.41) 

It is now straightforward to determine the Gaussian 
propagators. In particular, the order parameter correla- 
tion function reads 



with 

1 



AA„(fc) = 



Here 



ro-C7(fc,r!„) ■ 

4 ln(»o/(jJfc^ + !»„!)) 
1 + {KjH)ln{no/iDk^ + \^n\)) 



(4.42a) 



(4.42b) 



(4.42c) 



where f2o is an ultraviolet frequency cutoff and ro is the 
coefficient of the Vt^ in A^,. For asymptotically small 
wave numbers and frequencies the critical propagator is 
given by 



AA„(k) 



H -1 



const. 



in{nQ/{Dk^ + \n„\)) 



(4.42d) 



Am,-q = ci iVT / dx tr (/3(a;) Q{x)) . (4.40a) 



with const. > and r = rg — c\H/Kc the Gaussian 
distance from the quantum critical point. 

The Gaussian theory represented by Eqs. (|4.42() has 
several interesting properties. First, for given tq > 0, 
H, and Kc there is a quantum phase transition at a 
critical strength of oc \Kc\ which yields r = 0. Sec- 
ond, for fixed other parameters in the superconducting 
phase, an increase of K^, which is an increasing func- 
tion of disorder, will drive the system into the normal- 
metal phase. This is consistent with other theories that 
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have identified the Coulomb pseudopotential as an im- 
portant source of Tr-degra dation in disordered supercon- 
du ctors llFkLkelsteinl.ll987t [KirkDatr ick and BeUtl Il992t 
see iBehtz and Kirkpatrick, 1994 for a discussion of ear- 
her theories of this effect). Third, as in the ferromagnetic 
case there are long-ranged order parameter correlations 
in the disordered phase away from criticality. Equation 
Ij4.42d|) implies for the order parameter susceptibility 

> 0, |a;| ^ oo) cx l/|a;|'* In \x\. (4.43) 

A comparison with Eq. H4.25|l shows that the correlations 
are of even longer range than in the ferromagnetic case 
for all d > 2. This reflects the strong coupling of the su- 
perconducting order parameter to the generic soft modes 
mentioned in connection with Eq. (|4.37|) . 

Most of the critical behavior predicted by the 
Gaussian theory can s i mply be read off Eqs. (|4.42|l 
ijKirkpatrick and Belit j . Il997(l . The correlation length 
depends exponentially on r, rather than as a power law, 

For the correlation length exponent this implies 1^ = 00. 
The exponent 7 has its mean-field value, 

7=1, (4.44b) 

and the exponents rj and z for the order parameter (5*) 
and the fermionic {q) degrees of freedom, respectively, 
are 

?7* = 2 , 7]q = , z^^, ^ Zq ^2. (4.44c) 

The order parameter susceptibility at criticality decays 
in real space as 

X>]/(r = 0, |a;| -> 00) cx ln|a;|/|a;|'*. (4.45) 

As expected, the critical order-parameter fluctuations are 
of even longer range than those reflecting GSI in the dis- 
ordered phase, Eq. H4.43|l . 

These results also follow from a tree-level RG analysis 
of the field theory. The exponents rj are related to the 
scale dimensions of the fields via 

[q{x] = -{d-2 + r,,)/2, (4.46a) 
[^{x)] = -(rf-2 + r;*)/2, (4.46b) 

As in the magnetic case, Sec. IIV.B.4.c1 there is a critical 
fixed point where ci is marginal, and the fermions are 
diffusive, with exponents given by Eq. 14.44|l . However, 
in contrast to the magnetic case, the coupling constant C2 
of the term ,4^_q2 is RG irrelevant, and so are all higher 
terms in the expansion in powers of q. We therefore con- 
clude that the Gaussian critical behavior is exact. The 
most obvious technical reason for this surprising result 
is the fact that the time scales for the order-parameter 
fiuctuations and the fermions, respectively, are the same. 



which renders inoperative the mechanism that led the 
possibility of C2 being marginal in Sec. llV.lj.4l Physi- 
cally, the very long range of the order-parameter fluctu- 
ations reflecting the GSI stabilizes the Gaussian critical 
behavior. This is in agreement with the fact that long- 
ranged order parameter correlations in classical system s 
stabilize mean-fleld critical behavior (i Fisher et a/].ll972(l . 

The renormalized mean-fleld theory, the equation of 
state, and the critical exponent can be discussed in 
analogy to the magnetic case i Kirkpatrick and Belit j . 
1997; Zhou and Kirkoatrick, 200^. 

The theory of the metal-superconductor transition 
presented above is assuming conventional s-wave spin- 
singlet superconductors. It is interesting to ask how 
the phase transition scenario is modifled by exotic (i.e., 
nonzero angular momentum, ^ > 0) pairing. First of 
all, in contrast to conventional superconductivity, exotic 
superconductivity is rapidly destroyed by non magnetic 
disorder because Anderson's theorem ijAnderson, 195^ 
does not hold. This ha s been observed, e.g., in ZrZn2 
(iPfleiderer et a/.U2001b(l which is believed to be a p-wave 
spin-triplet superconductor, i = 1. Recently, the re- 
sulting quantum phase transition between a dirty metal 
and an exotic superconductor has been studied within 
a La ndau-Ginzburg- Wilson approach ( Skneonck et all 
It turns out that the nonzero order parameter 
angular momentum suppresses the effects of GSI on this 
transition. The nonanalytic wave number dependence in 
the Gaussian order parameter propagator Mn{k) takes 
the form fe^^lnfc. Thus, compared to the s-wave case, 
Eq. Il4.42dll , the nonanalytic term is suppressed by a fac- 
tor of fc2^. This can be understood as follows: In the 
presence of nonmagnetic quenched disorder, the domi- 
nant electronic soft modes are those that involve fluctu- 
ations of the number density, spin density, or anomalous 
density in the zero angular momentum channel, while 
the corresponding densities in higher angular momen- 
tum channels are not soft. Since the different angular 
momentum modes are orthogonal at zero wave number, 
the coupling between a flnite angular momentum order 
parameter and the zero angular momentum soft modes 
must involve powers of the wave number |fc|. 

In conclusion, for £ > 1, the GSI induced nonanalytic 
term is subleading compared to the conventional fe^ term, 
and it will not influence the critical behavior. For 1=1, 
the nonanalytic term is marginally relevant. The ulti- 
mate fate of the transition has not yet been worked out, 
mainly because of (unrelated) complications stemming 
from the disorder fluctuations similar to that in disor- 
dered itinerant antiferromagnets. 



D. Quantum Antiferromagnetic Transition 

Let us now consider a quantum antiferromagnetic tran- 
sition, in analogy to the ferromagnetic one discussed in 
Sec. lIII.B.ll A crucial difference between these two cases 
is that in the latter, both the order parameter fleld M 
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(whose average is the magnetization) and the generic soft 
modes are soft at zero wave number, while in the former, 
the order parameter field N (whose average is the stag- 
gered magnetization) is soft at a nonzero wave number 
As a result, the hydrodynamic wave number k in 
the dynamical piece of the ferromagnetic Gauss ian ac- 
tion, Eq. (|4.6|l . gets replaced by p, and one has (|Hertj . 
[1976) 

Au[N] = J2 ^(k) {r + ck'^ + d laj) Nik) + 0(Ar4). 

k 

(4.47) 

Here k is the wave vector measured from the reference 
wave vector p. The missing inverse wave number in the 
frequency term reflects the much weaker coupling, com- 
pared to the ferromagnetic case, of the particle-hole ex- 
citations to the order parameter field. Because of this 
weaker coupling, one expects the Gaussian approxima- 
tion to be much better here than in the ferromagnetic 
case. Indeed, it is easy to show that the effects that 
needed to be taken into account for the ferromagnet are 
irrelevant with respect to the mean-field transition de- 
scribed by Eq. H4.47|l . One thus expects a continuous 
transition with mean-field critical behavior in all dimen- 
sions d > 2. The finite- temperature pro perties of thi s 
theory have been worked out in detail bv iMillisI (^2^. 
It is important to remember that Eq. (|4.47|l assumes that 
the only relevant soft modes at the quantum antiferro- 
magnetic transition are the order-parameter fluctuations. 

Experimental observations arc not in agreement with 
this expectation, probably because most metallic ma- 
terials that display easily accessible quantum antifer- 
romagnetic transitions are far from being simple met- 
als. They fall into the class of heavy-fermion ma- 
terials; an overview has been given by Coleman et al. 
l)200l|) . The best-studied system is CeCug-xAux, which 
shows a quantum phase transition to an antiferro- 
magi ietic state at a criti c al gol d concentration Xc ~ 
0.1 ijvon Lohnevsen et all Il994j) . There are experi- 
mental indications for the quantum critical fluctua- 
tions being two-dimension a l or quasi-two-dimensional 
in nature l|Stockert et oil Il99i^ . A detailed phe- 
nomenological analysi s of neutron scattering experiments 
l|Schr5der all\200(l(\ has shown that the magnetic sus- 
ceptibility is well described by the form 

X{k, n) = a [{-in + bT)" + 6i(fc)"]~^ . (4.48) 



This holds for a generic shape of the Fermi surface. Special 
geometric features of the Fermi surface can cause particle-hole 
excitations to be soft at the same wave vector as the order pa- 
ramet er even in an antiferromagnet, see. lAbanov and Chubukovl 
120001) and the discussion below. 

An interesting counterexample is the case of Cri—x^x, which is 
a common transition metal, yet displays properties, both at its 
antiferromagnetic quantum critical point and away from it, that 
is very similar to the exotic behavior sho wn by the heavy-fe rmion 
materials and high-Tc superconductors JYeh et all. l2002l) . We 
will come back to this in Sec. |V] 
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FIG. 24 The magnetic structure factor in CeCus.gAuo.i as 
measured by neutron scattering at various frequencies and 
temperatures, as a function of Q/T. The inset shows a mea- 
sure of the scatter of the scaling plot for various values of the 
exponent a. After fSchrodor et al... (,200Q') . 



Here a and b are constants, a — 0.75 ± 0.05, and 9{k) is 
the wave vector dependent Weiss temperature. This type 
of behavior is often referred to as "local quantum criti- 
cality" , as it is believed (although not universally so, see 
below) to point to atomic-scale physics as the cause of the 
antiferromagnetism, as opposed to the Fermi-liquid or 
spin-density- wave description that underlies Eq. H4.47|l . 
The possibility of s uch a loca l quan tum critical point was 
first pointed out bv ISi et al\ ()l999(l . 

Perhaps the most obvious violation of the expected 
mean-field behavior is the occurrence of the exponent 
a ^ 1. Other discrepancies between the observed behav- 
ior and the one expected fr o m He rtz t heory ha v e bee n 
discussed bv lColeman et al\ (|2001[1 and ISi et all l|200lD . 
A prominent one is that frequency and temperature are 
observed to scale in the same way, see Fig. |21 This ob- 
servation is usually referred to as "fi/T scaling", and it 
is expected to hold at a quantum critical point that ex- 
hibits hyperscaling (Sachdcv and Yc, 1992). By contrast, 
in Hertz theory scales as T^/^ (Mlllis, 1993). Another 
one is the behavior of the specific heat coefficient, which is 
observed to diverge logarithmically, while in Hertz theory 
it remains finite a nd shows a square-roo t cusp singular- 
ity llMilliaLll99.'l) . ISchroder et all l)2000|) have proposed 
the following physical picture to explain this violation. 
In CeCuAu, or any heavy-fermion material, the highly 
localized f-orbitals of the lanthanide (in this case, Ce) 



Naively, one might expect Q/T scaling to simply follow from 
the fact that the scale dimension of T is given by the dynamical 
critical exponent z, and hence the same as that of Q. In general, 
this is not true due to the presence of, (1) multipl e temperatur e 
scales, and (2) dangerous irrelevant variables. See. lMillia 

(1223); 

ISachdevI il997D . 
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FIG. 25 Schematic phase diagram for CeCug-ajAu^. Shown 
are the antiferromagnetic (AFM) phase, the local moments 
(LM) region, and the Fermi liquid (FL) region. The Neel 
temperature Tn (solid lines) represents a true phase transi- 
tion, while the Kondo temperature Tk (dashed lines) denotes 
a crossover. The three scenarios shown correspond to the 
three different physic al situations discussed in the text. After 
ISchroder ei"airil2000fl . 



or actinide provide local magnetic moments, while the 
more extended s-, p-, and d-orbitals provide a Fermi sur- 
face. An increase in Au doping increases the hybridiza- 
tion between the localized and extended orbitals, which 
increases the coupling between the local moments via the 
conduction electrons, and hence the Neel temperature 
Tn. If the hybridization becomes too strong, however, 
the f-electrons are incorporated into the Fermi surface, 
and the antiferromagnetism disappears. As a function 
of doping, Tn thus goes through a maximum and disap- 
pears at a critical doping concentration, see Fig. If 
this were all that happens. Hertz theory should apply. 
However, there is a second temperature scale besides Tn, 
namely, the Kondo temperature Tk. Only for T < Tk do 
the local moments of the f-orbitals become screened by 
the condu ction electron s, and the hybridization becomes 
effective. ISchroder et a l. (2000) have proposed that Tk 
vanishes at the same critical concentration as Tn, see 
Fig. I25r b). In this picture, there is no heavy-electron 
Fermi surface at the quantum critical point, and Hertz 
theory is inapplicable. Notice that, according to this pic- 
ture, there is nothing technically wrong with Hertz the- 
ory; the reason for its failure is that the model does not 
contain crucial slow degrees of freedom. Alternatively, 
one can imagine a scenario where Tk remains nonzero 
into the antiferromagnetic phase, see Fig. TiSf a). but the 
experiments show that this is not the case in CeCuAu. 
Finally, it is conceivable, at least in principle, that the 
local-moment region might extend to zero temperature in 
an entire region of the phase diagram, as shown in Fig. 
I25r c). although there is currently no explicit theory that 
can explain why the unscreened local moments in heavy- 
fermion systems would not order at sufficiently low tem- 
peratures. Still, it is interesting to note that there are 
cases where a non-Fcrmi liquid phase is observed adja- 
cent to an antiferromagnetic one, e.g., in URu2-a;Rea;Si2, 
see Fig. |21 

Theoretically, the quantum antiferromagnetic transi- 
tion is an open problem that is under very active con- 
sideration. A theoretical description of the above sce- 
nario involves the so-called Kondo lattice problem, i.e.. 



the interaction of many localized spins with each other 
and with a band of conduction electrons. Si and col- 
laborators have studied such a model within a dy nam- 
ical mean-field approach HI, l2003t ISi et all l200l|) . In 
this approximation, these authors find what they call 
a local quantum critical point, which has many prop- 
erties that are consistent with the observations. The ex- 
ponen t a is nonuniversal. More recently, Scnthil et al\ 
l)2004j) have proposed that the magnetic state is an un- 
conventional spin-densit y wave where s pin-charge sepa- 
ration has taken place. IColeman et al\ (|2000) have de- 
veloped a supersymmetric representation of spin oper- 
ators that allows for a treatment of both magnetism 
and the Kondo effect in the context of a large- expan- 
sion. These authors have speculated that, if applied to 
the Kondo lattice problem, this formalism can give rise 
to two different types of fixed points: A weak-coupling 
fixed point of Hertz type, and a non-Fermi-liquid fixed 
point that displays spin-charge separation. There is, 
however, no consensus that the existence of a Fermi 
surface prec l udes an explanation of the observations. 
iRosch et al\ l)l997(l have proposed an explanation in 
terms of three-dimensional conduction electrons coupling 
to two-dimensional ferromagneti c fluctuations. In a two- 
dimen sional spin-fermion model, Eb anov and Chubukovl 
(2000) have found that special nesting properties of the 
Fermi surface lead to an exponent a ~ 0.8, but their 
results do not display il/T scaling. A breakdown of 
the LGW expansion due to nesting has also been dis- 
cussed bv jLe rchcr a nd Wheatlev (2000), for a revie w, see 
lAbanov et airij2003|) . fSachdev ai^dMoriMjil l)2002|) have 
found that similar "nonlocal" forms of the dynamic sus- 
ceptibility are obtained in two-dimensional models where 
an order parameter couples to long- wavelength deforma- 
tions of a Fermi surface. Notice that all of these the- 
ories involve some coupling of soft modes, generic or 
otherwise, to the order parameter, although there cur- 
rently is no consensus about their nature and origin. In 
the order of theoretical ideas listed above, they are, (1) 
Slow local moment fluctuations that serve as "Fermi sur- 
face shredders", or (2) other fluctuations that destroy 
the Fermi surface, ^'^ (3) ferromagnetic fluctuations, (4) 
particle-hole fluctuations across a Fermi surface with spe- 
cial geometric features, and (5) volume and shape defor- 
mation of the Fermi surface. We finally mention that an- 
tiferromagnetic quantum criticality, in particular in two- 
dimensions, has received much attention in connection 
with high-Tc superconductivity, a topic that is beyond 
the scope of this article. 



These soft modes are 'generic' in the sense of our definition only 
in the scenario depicted in panel (c) of Fig. 1251 
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V. DISCUSSION AND CONCLUSION 

In this section we summarize the main theoretical ideas 
we have presented, as weh as the experimental situation. 
We also discuss a number of open problems in the field 
of quantum phase transitions, and suggest a number of 
new experiments to address some of these. 



A. Summary of Review 

Early work on quantum phase transitions had sug- 
gested that most of them were conceptually quite simple, 
since they are related to corresponding classical phase 
transitions in a higher dimension. This led to the con- 
clusion that quantum critical behavior generically would 
be mean-field like. For a number of reasons, these conclu- 
sions have turned out to be not valid in general. In this 
review we have discussed one of the mechanisms that in- 
validate the mapping of a quantum phase transition onto 
a simple classical one in higher dimensions. The central 
idea behind this mechanism is as follows. The soft-mode 
spectrum of many-body systems is in general different at 
T = from the one at T > 0, since there are soft modes 
at T = that develop a mass at nonzero temperature. 
These soft modes lead to long-ranged correlations in en- 
tire regions of the phase diagram, a phenomenon known 
as generic scale invariance. Physically, both the generic 
soft modes and the critical order-parameter fluctuations 
are equally important in the long- wavelength limit, and if 
the coupling between them is sufficiently strong, the for- 
mer will influence the leading critical behavior. This can 
happen for classical phase transitions as well, but it is less 
common since at T > there are fewer soft modes than 
at T = 0. These observation led to the general paradigm 
that effects related to generic scale invariance are of fun- 
damental importance for the theory of generic quantum 
critical points. As a result, quantum phase transitions 
are typically related to classical phase transitions in the 
presence of generic scale invariance, rather than to simple 
ones. 

In Sec. ^ we have reviewed the concept of generic 
scale invariance in both classical and quantum systems. 
Although in both cases examples are plentiful, in the 
quantum case this is especially so because of additional 
Goldstone modes that exist at zero temperature. We 
have distinguished between direct and indirect generic 
scale invariance effects; the former being immediate con- 
sequences of Goldstone's theorem, conservation laws, or 
gauge symmetries, while the latter arise from the former 
via mode-mode coupling effects. Classical examples we 
have discussed include: Goldstone modes in Heisenberg 
ferromagnets or analogous systems, which lead to long- 
ranged susceptibilities everywhere in the magnetically or- 
dered phase; local gauge invariance in superconductors or 
liquid crystals, in which context we have stressed inter- 
esting analogies between statistical mechanics and parti- 
cle physics; long-time tails in time correlation functions 



which determine the transport coefficients in a classical 
fluid in equilibrium; and long-ranged spatial correlations 
in a classical fluid in a nonequilibrium steady state. Ex- 
amples of generic scale invariance in the quantum case 
all involved interacting electron systems, namely, weak- 
localization effects in disordered materials and the anal- 
ogous effects in clean ones; and nonequilibrium effects 
analogous to those in classical fluids. Throughout this 
discussion we have stressed the coupling between the 
statics and the dynamics in quantum systems, which 
leads to long-ranged spatial correlations in quantum sys- 
tems even in equilibrium. 



In Sec. IIIII we discussed some classical phase tran- 
sitions where generic scale invariance plays a central 
role. Our flrst example was the nematic-smectic-A tran- 
sition in liquid crystals, which maps onto the classical 
superconductor-normal metal transition. In this case the 
relevant generic soft modes are the director fluctuations, 
which are Goldstone modes due to a broken rotational 
symmetry. They can drive the transition first order, even 
though in their absence one expects a continuous transi- 
tion. The other classical phase transition discussed was 
the critical point in a classical fluid. Here the generic 
soft modes are due to conservation laws. In equilibrium, 
they influence the critical dynamics only. However, in a 
fluid subject to shear, they also couple to the static criti- 
cal behavior. In this case, long-ranged static correlations 
due to the generic soft modes stabilize mean-fleld critical 
behavior below the equilibrium upper critical dimension. 



Section IIVI was devoted to four examples of quan- 
tum phase transitions. For the flrst three, namely, the 
quantum ferromagnetic transition in clean and disor- 
dered itinerant electrons systems, respectively, and the 
normal metal-superconductor transition at T = 0, the 
generic soft modes are the particle-hole excitations that 
cause the weak-localization effects and their clean coun- 
terparts. The clean ferromagnetic case turned out to be 
analogous to the classical nematic-smectic-A transition 
in that the generic soft modes can lead to a fluctuation- 
induced flrst-order transition. The disordered quantum 
ferromagnetic transition is analogous to the classical fluid 
under shear, in the sense that static long-ranged corre- 
lations due to the generic soft modes stabilize a simple 
critical fixed point, and the critical behavior can be de- 
termined exactly even in d = 3. In the case of the su- 
perconducting transition this effect is even stronger, and 
the critical fixed point is Gaussian. Our fourth exam- 
ple dealt with the quantum antifcrromagnetic transition, 
where the critical behavior is currently not understood. 
There are many indications that, in the materials stud- 
ied so far, generic soft modes related to local magnetic 
moments exist and couple to the critical order-parameter 
fluctuations, but a detailed theory of this effect remains 
to be worked out. 
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B. Open Problems, and Suggested Experiments 

We conclude by listing a number of open questions, 
both theoretical and experimental, the answers to which 
would help shed light on some of the problems we have 
discussed. Our remarks are necessarily incomplete and 
speculative. 



1. Generic scale invariance, and soft modes 

There are strong indications that the list of mecha- 
nisms for generic scale invariance in Sec.^is incomplete, 
even for relatively simple systems. Here we briefly men- 
tion two examples. 



a. Generic non-Fermi-liquid behavior As mentioned in 
Sec. lIV.B.4l behavior not consistent with Landau's Fermi 
liquid theory has been observed in a variety of ma- 
terials far from any quantum critical point. While 
such exotic behavior may not come as too much of a 
surprise in strongly correlated systems involving rare 
earths or actinides, the well documented case of MnSi 
ijPfleiderer et all l2001a ^ makes it likely that some more 
basic understanding is lacking. The observations clearly 
show that long-ranged correlations exist in a large region 
of parameter space that are not caused by any quan- 
tum phase transition. For instance, the resistivity has 
a T'^/^ asymptotic temperature dependence up to a fac- 
tor of two in parameter space away from the quantum 
ferromagnetic transition. This generic long-time tail be- 
havior cannot be explained in an obvious way by any ex- 
isting theory. Apart from being intrinsically interesting, 
this also raises the question whether there are unknown 
generic soft modes that will be important for a complete 
understanding of the quantum phase transition. One ob- 
vious candidate are slow fluctuations of local magnetic 
moments, which have been invoked for many mysteri- 
ous effects, from hard-to-understand aspects of metal- 
in sulator transitions (for a rev iew see, e.g.. Sec. IX. B 
of lBelitz and Kirkpatrickl Il994|) to the critical behavior 
at the quantum antiferromagnetic transition (see. Sec. 
I1V.D|I . So far, however, there is no detailed theory that 
is capable of describing these soft modes. 



b. Soft modes due to nested Fermi surfaces 
EE anov and Chubukol l|2000(l have discussed nest- 
ing properties of Fermi surfaces that are important for 
quantum antiferromagnetic transitions. In the language 
of this review, the effect considered by these authors 
is likely a change of the properties of the generically 
soft particle-hole excitations. It would be useful to 
explicitly confirm this conjecture, and to develop a 
general classification of the effects of Fermi surface 
geometries on generic soft modes. 



2. Aspects of the quantum ferromagnetic transitions 

a. Disordered ferromagnets A detailed experimental 
study of the critical behavior of both the thermodynamic 
and transport properties near the quantum ferromagnetic 
transition in a disordered itinerant electron system would 
be very int eresting. Some recent results on Fei_j,Coj;S2 
l)DiTusa et ai. 2003) appear to be consistent with the 
log-log-normal corrections to power-law scaling discussed 
in Sec. lIV.Bl A more precise determination of the critical 
behavior in this or other systems would be of great help 
in confirming or refuting the theoretical ideas. 

Also of interest would be a systematic study of the 
destruction of the first order ferromagnetic transition by 
nonmagnetic disorder that is predicted by the theory dis- 
cussed in Sec. IIV.BI For instance, the theory predicts 
that the tricritical point observed in MnSi and UGe2 
will move to lower and lower temperature with increasing 
strength of quenched disorder, turn into a tricritical end 
point, then a critical end point, and finally the first-order 
transition should disappear, see Fig. 1201 An experimen- 
tal check of this prediction would be very valuable. For 
instance, it would give an indication of whether or not 
an understanding of the non-Fermi-liquid nature of the 
paramagnetic phase, mentioned in Sec. IV.B.II is impor- 
tant for describing the quantum phase transition. 



b. Itinerant ferromagnets with magnetic impurities An in- 
teresting quantum phase transition to study experimen- 
tally would be the ferromagnetic transition in an itin- 
erant electron system with dilute magnetic impurities, 
which would fundamentally change the soft-mode struc- 
ture in systems with or without additional nonmagnetic 
disorder. Such a study would provide an important 
check of the soft-mode paradigm, especially if it were 
possible to start with no magnetic impurities and study 
the crossover induced by their gradual introduction. In 
the absence of complications due to Kondoesque effects, 
one would expect a local LGW theory for the order- 
parameter fluctuations to be valid, since the magnetic 
impurities cut off the soft modes that strongly couple to 
the order-parameter ffuctuations. Without nonmagnetic 
disorder, this would be Hertz's original theory, while in 
the nonmagnetically disordered case, rare-region effects 
would likely play an important role. 

It is likely, however, that Kondo screening effects, and 
interactions between the impurity sites, will lead to com- 
plications not unlike those believed to be responsible 
for the observations in quantum antiferromagnets. Sec. 
IIV.DI which are not understood. Generally, the physics 
of local moments, and their influence on the properties 
of conduction electrons, is one of the most important un- 
solved problems in condensed matter physics. 



c. Clean ferromagnets The finite-temperature behavior 
of the observables near the continuous quantum ferro- 
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magnetic transition in clean itinerant electron systems 
has not been worked out. Given the existing RG descrip- 
tion that was briefly discussed in Sec. IIV.A.4I and the 
existing theory of finite-temperature effects near Hertz's 
fixed point UMillis. 1993,) . this should be a relatively 
straightforward problem. 



3. Metal-superconductor transition 

a. Experimental study of critical behavior It would be in- 
teresting to experimentally study the quantum phase 
transition between a disordered metal and a disordered 
conventional bulk superconductor. The theory reviewed 
in Sec. lIV.Cl predicts that the quantum phase transition 
is governed by a Gaussian fixed point. If this is correct, 
there will be no crossover to mean-field behavior, and 
the asymptotic critical behavior will be observable, in 
contrast to the situation at the thermal transition. Since 
the predicted Gaussian critical behavior is radically dif- 
ferent from mean-field critical behavior, this should be 
easy to observe. 



b. Quantum phase transition to a gapless superconductor 
The metal-superconductor transition described in Sec. 
llV.CI assumed that the superconducting states was a con- 
ventional gapped state. It would be interesting to con- 
struct an analogous theory for the case of a gapless su- 
perconductor , especial l y in t he light of the suggestion 
bv fVishvcsh wara et all l)200(i) that the superconducting 
ground state in a disordered system will generically be 
gapless. On general grounds it is likely that the critical 
behavior is be the same in both cases if the transition is 
approached from the metallic side. However, the gapless 
superconducting state will have additional soft modes, 
and these modes may influence some of the critical be- 
havior, such as the exponents governing the equation of 
state, when the transition is approached from the super- 
conducting side. 



4. Aspects of the quantum antiferromagnetic transitions 

a. Phenomenological theory for the quantum antiferromag- 
netic problem As discussed in Sec. IIV.DI the quan- 
tum antiferromagnetic transition in real systems remains 
rather incompletely understood. In most materials stud- 
ied so far, a basic ingredient of the problem seems to be 
itinerant electrons coupled to local moments. The latter 
are not screened since, at or near the quantum critical 
point, the Kondo temperature vanishes due to critical 
fluctuations effects. In the absence of a more fundamen- 
tal approach to this very complicated problem, it would 
be very interesting to construct a purely phenomenolog- 
ical theory of these coupled fluctuations, including the 
soft modes associated with the slow temporal decay of 
the local moments. 



b. Simple quantum antiferromagnets As we mentioned in 
Sec. IIV.UI most of the work on the quantum antiferro- 
magnetic transition has focused on heavy-fermions com- 
pounds and high-Tc superconductors, partly due to the 
general interest in these systems. These materials are far 
from being simple metals, and their behavior is heavily 
influenced by local-moment physics, among other com- 
plications. It would be interesting to flnd and study an 
itinerant electron system without local moments that has 
a quantum antiferromagnetic transition, provided such 
materials do indeed exist. The complicated behavior 
of Cri-xVx, see footnote |^ is discouraging in this re- 
spect, although at least some properties of this material 
have been exp lained as due t o nesting properties of the 
Fermi surface iB azaliv e t all l2004HNorman et a/.l. [20031 
iPepin and Norman. 2004j) . In the absence of disorder, 
the quantum phase transition in such a systems should 
be described by Hertz's theory. The introduction of non- 
magnetic disorder would likely lead to a quantum phase 
transition where statistically rare events are important. 



5. Nonequilibrium quantum phase transitions 

In Sec, nil. B. 21 we discussed a classical nonequilibrium 
phase transition where generic scale invariance plays a 
role. The general topic of classical phase transitions in 
driven s ystems has received a substa ntial amount of at- 
tention llSchmittmann and Zi ilHil. In Sec. lII.ij.5l we 
discussed how nonequilibrium situations at zero tempera- 
ture can lead to correlations that are of even longer range 
than in equilibrium. Experimentally, however, the prob- 
lem of nonequilibrium quantum phase transitions has re- 
ceived little attention so far. In quantum Hall systems, 
which we have not discussed in this review , elect ric-field 
effects have been studied (see, ISondhi et aA\l99(t\ . but in 
this case the electric field acts as a relevant operator with 
respect to the transition. Experiments that show an ac- 
tual nonequilibrium quantum phase transition would be 
of interest. 



6. Other quantum phase transitions where GSI might play a 
role 

In addition to the examples covered in this review, 
one can think of quantum phase transitions that have 
not been investigated so far either theoretically or ex- 
perimentally, for which GSI effects are likely to play a 
central role. For example, the order parameter for the 
spin-triplet analog of the isotropic-to-nematic phase tran- 
sition that has been proposed to o ccur in quantum Hall 
systems bv lOeanesvan et oil l)200l|) is expected to couple 
to generic soft modes. For this transition, in contrast to 
its spin-singlet analog, a simple LGW theory is therefore 
unlikely to be valid. 
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